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Introduction 


Throughout this module you will use computer software called Maxima to 
investigate mathematical concepts and perform calculations that may be 
too laborious or difficult to do by hand. Maxima is a computer algebra 
system (often abbreviated to CAS); that is, a computer program that 
can perform algebraic manipulations, such as expanding brackets, 
factorising expressions and solving equations, as well as performing 
numerical calculations and plotting graphs. It is typical of the type of 
computer software used by mathematicians today. 


This Guide explains how to use Maxima for some of the mathematics in 

MST124. It is designed to be read in conjunction with the main study 

texts and is not meant to be read all at once because its later sections use 
mathematical concepts that you may not have yet met. As you work EE 
through the main study texts in Books A-D, you will come across sections EAS 
or activities marked with the icon shown here. Some of these ask you to 

study part of this Guide. Other activities ask you to solve a problem using 

Maxima commands that you have previously learned. While the 

mathematical answers to such activities are included within the 

appropriate book, the computer methods needed to solve them are given in 

the section ‘Computer methods for CAS activities in Books A-D’ towards 

the end of this Guide. 


The ‘Maxima accessibility guide’ (Section 14) is primarily aimed at those 
who may have difficulty using Maxima because of a disability, and provides 
some advice and suggestions. 


The final section of this Guide contains a summary of all the Maxima 
commands introduced. This should be useful for reference, both as you 
work through the module and afterwards. 


Note that this Guide does not attempt to describe all the features of 
Maxima. If you need more information, please consult the extensive 
documentation provided on the Maxima website at 
maxima.sourceforge.net or use Maxima’s help system, which is 
described in Subsection 2.10. 


Maxima’s origins lie in a system called Macsyma, a well-respected 
computer algebra system first developed in the late 1960s at the 
Massachusetts Institute of Technology and funded by the United 
States Department of Energy. From 1982 until his sudden death in 
2001, William (Bill) Schelter, a professor of mathematics at The 
University of Texas at Austin, maintained a version of Macsyma 
which we now know as Maxima. In 1998 he obtained permission from 
the Department of Energy for the system to be made freely available 
to all. Maxima is now dedicated to his memory. William Schelter (1947-2001) 
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“Rather than learning how to solve that, 
shouldn't | be learning how to operate 
software that can solve that?” 


A word of warning 


Using a computer algebra system is not a substitute for learning the 
mathematical techniques taught in this module! You need to become 
proficient at using the methods taught so that you can apply them in the 
examination and in your future studies, and also to help you develop your 
mathematical understanding and intuition. 


Maxima is not intelligent — it is just a tool, similar to your calculator. You 
need to use your own mathematical skills to guide it and interpret its 
results critically. 


Conventions 


Here are the various computing terms used in this Guide, and how they 
should be interpreted on Microsoft Windows and Apple Mac computers. 
The Guide concentrates on using Maxima on a Microsoft Windows 
computer, although Maxima is available for many different types of 
computer. 


Computer terminology 


Term Microsoft Windows Apple Mac 


Click 


Right click 


Select 


Ctr1-Q 
Etc. 


Click the left mouse/touchpad button Click the mouse/touchpad button 


Click the right mouse/touchpad button Hold the key while clicking 


the mouse/touchpad button 


Either click to select an option, Either click to select an option, 


or select text by dragging the cursor over 
the text while pressing the left mouse button 


Hold down while pressing (Q) 


or select text by dragging the cursor over 
the text while pressing the mouse button 


Hold down while pressing (Q) 


Etc. Etc. — but there are exceptions to this 


pattern that we’ll mention as they arise 


On some computer keyboards the key might be labelled i 
or simply , and the key might be labelled : 

In this Guide, these keys are referred to as (_ Enter ) and Shift _}, 
respectively. 

Sometimes | Enter ) and Return } are two separate keys. For such keyboards 
references in this Guide to the key should be interpreted as 


references to the key. 
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1 Installing, testing and configuring 
Maxima 


1.1 Introducing and installing Maxima 


Maxima can be used on many different types of computer, and there are 
several different intermediary systems, known as interfaces or 
front-ends, that provide a means of communication between you, the 
user, and the Maxima system itself. This is illustrated in Figure 1. The 
interfaces enable you to enter commands that are then processed by 
Maxima, and they display the results obtained. 


Interfaces 


wxMaxima 


XMaxima 


Maxima 


Command line 


Figure 1 Different interfaces for Maxima 


All the commands mentioned in this Guide can be used with any interface 
to Maxima, unless indicated otherwise. 


The recommended interface is wxMaxima, which is shown in Figure 2. 
(The letters ‘wx’ in the name of this interface refer to the software used to 
create it.) If you need to use an interface other than wxMaxima, for 
example if you are not using a Windows computer and you are unable to 
install wxMaxima, then you should try to achieve all the outcomes shown 
in this Guide using your chosen interface. 


If you use a screenreader, then the basic command-line interface shown 
in Figure 3 may be the easiest interface for you to use. 


If you think you may have difficulty using the wxMaxima interface, 
because of a disability for example, then you should read Subsection 14.1 
now. This contains details of how wxMaxima might be configured to meet 
your needs, and should help you decide whether to use wxMaxima or the 
command-line interface. 
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File Edit View Cell Maxima Equations Algebra Calculus Simplify Plot Numeric Help 


Dy 3 = | 
S@ASelixHOrlaleoes| |© 
(311) expand ( (2*x+1) * (S-3*x%2) ); 
(%01) -6 x9 -3 x?+10 x+5 
Welcome to wxMaxima Ready for user input 


Figure 2 The wxMaxima interface 


E command line Maxima 


Maxima 5.30.0 http://maxima.sourceforge.net 

using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING. 
Dedicated to the memory of William Schelter. 

The function bug report() provides bug reporting information. 
(Sil) 


E 


an) > 


Figure 3 The command-line interface to Maxima 


Before you can use Maxima, you need to install it on your computer. You 


will do this in the next activity. 


Computer activity 1 Installing Maxima 


(a) Go to the OU Maxima website at learn1.open.ac.uk/site/maxima 
and follow the instructions there to download and install Maxima on 
your computer. (You can also find a link to the OU Maxima website 


on the MST 124 website.) 


1 Installing, testing and configuring Maxima 


(b) Create a folder on your computer in which to store your Maxima work 
for this module. You might like to put this on your computer desktop, 
or in your usual documents area. 


If you are studying within a prison or other closed institution, you can 
install Maxima from the Offline resources disc. 


1.2 Starting, testing and configuring 
wxMaxima 


If you are using the command-line interface you should study 

Subsection 14.2 instead of this subsection. If you are using a different 
interface, you should make sure that you can achieve similar outcomes with 
that interface. 


Computer activity 2 Starting wxMaxima 


Start wxMaxima on your computer, and keep it open to work with as you 
study this section. 


If you cannot remember how to do this on your computer, see the 
OU Maxima website: learnl.open.ac.uk/site/maxima. 


The wxMaxima interface to Maxima, illustrated in Figure 4, consists 
mainly of a large white working area (called the worksheet) where you 
can type Maxima commands. These commands instruct Maxima to either 
perform calculations or change system settings. 


When the software first starts, a secondary window showing the ‘Tip of the 
Day’ is also displayed. You can close it by clicking the ‘Close’ button. If 
you want to prevent this window from appearing the next time you start 
wxMaxima, click the box in its bottom left-hand corner, so that it is no 
longer ticked, before clicking ‘Close’. (You can always reinstate these 
messages at any time by selecting Show tips from the Help menu.) 


Along the top of the main window (or along the top of the screen in the 
case of an Apple Mac) is a list of menu names, including the File menu 
which allows you to save or print your work and re-open previously saved 
sessions (there is more on this later). You can also close the software from 
the File menu by selecting Exit. (On an Apple Mac you can do this by 
selecting Quit wxMaxima from the wxMaxima menu.) 


In Figure 4, a row of icons known as the toolbar is shown below the row 
of menu names. This is displayed only if Main Toolbar is selected in the 
View menu. The icons provide shortcuts to commonly used functions. 


The View menu also allows you to display palettes of commonly used 
commands. You can use these palettes as an alternative to typing some of 
the commands given in this Guide. 
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o laima 16.04.2 [ unsaved 


File Edit View Cell Maxima Equations Algebra Calculus Simplify Plot Numeric = 


W ERE APEL ER =|! |© 
| ~ |Table of Contents X 
Toolbar Menus 
Tip of the Da 
® Did you know... 
To start using wxMaxima right away, start 
typing your command. An input cell should 
appear. Then press Shift-Enter to evaluate your 
command. 
© Show tips at startup Next Tip Close 
Status line 
2 | 
Welcome to wxMaxima Ready for user input 


Figure 4 The initial wxMaxima window 


Also shown in Figure 4 is the Table of Contents. This shows the names 
of any sections or subsections that you may add to your worksheet (see 
Subsection 2.8). You may want to close the Table of Contents, by clicking 
the cross near its top right-hand corner. 


At the bottom right-hand corner of the wxMaxima window is a message 
that indicates the current status of the Maxima system. When Maxima is 
waiting for instructions from you this reads ‘Ready for user input’. Other 
possible messages include ‘Maxima is calculating’, or ‘Reading Maxima 
output’. 


The following activity shows you how to use wxMaxima to perform a 
simple calculation and checks that your system is working properly. Follow 
the steps using wxMaxima yourself as you read the activity. 


If you encounter unexpected problems when working through this activity, 
check the Frequently asked questions (FAQs) section on the OU Maxima 
website: learnl.open.ac.uk/site/maxima. 


Computer activity 3 Your first Maxima calculation 


Follow these steps to use Maxima to calculate 2 + 3 x $. 


(a) Click anywhere on the large white area of the wxMaxima window, to 
ensure that the text you type next is directed to the correct place. 


1 Installing, testing and configuring Maxima 


(b) Type the following exactly as it appears here: 
2+3%4/5; 


including the semicolon (;) at the end. 


Notice that as you type your command, it is shown in 
wxMaxima preceded by a red arrow: -->. 


This is the input prompt and shows where commands can be typed. 


To the left of the input prompt is the cell marker 


[ 


which will be described in more detail later. 


As you type, the colour of the text changes to highlight the different 
elements of your expression: numbers are coloured differently to the 
operations. 


(c) Press Shift-Enter, that is, hold down the key while pressing 


| 


(Remember: if, as on some Macs, your keyboard has separate 


and keys, then you should hold down the key while 
pressing .) 


If you make a mistake, try clicking elsewhere on the worksheet, below 
your error, and try again. 


Pressing Shift-Enter instructs wxMaxima to send the expression to 
Maxima to calculate. This is known as evaluating the expression. 


The input prompt (-->) is now replaced by (%i1). This is how 
Maxima labels lines. The i stands for ‘input’: this is input line 
number 1. Within each Maxima session line numbers start at 1 and 
increase for each new input line. The percentage symbol (%) is used 
by Maxima to indicate objects built into the system. There is more 
about this later. 


The status message in the bottom right-hand corner of the wxMaxima 
window first changes to ‘Maxima is calculating’, then ‘Reading 
Maxima output’, before the result of the calculation is displayed. The 
message then displays ‘Ready for user input’ when the system is ready 
to accept another command. 


The result is displayed under the input line, preceded by the label 
(%01). The o stands for ‘output’, so this label identifies the line as 
output line 1. It gives the output corresponding to input line 1. 
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22 
($01) = 
5 


= 2+3*4/5; 


Welcome to wxMaxima Ready for user input 


If, in Computer activity 3, Maxima behaved as described, then you have a 
working system to use throughout MST124. Before finishing though, you 
should make a couple of small changes to how the wxMaxima interface 
behaves. 

In Computer activity 3 you used Shift-Enter rather than to 
instruct Maxima to evaluate expressions, and the input prompt (-->) 
appeared only once you started typing. Both these behaviours can be 
disconcerting, but can be changed, as described in Computer activity 4. 


It is recommended that you configure wxMaxima as described in 
Computer activity 4, and in the remainder of this Guide it is assumed 
that you have done so. 


Computer activity 4 Configuring wxMaxima 


(a) From the Edit menu, select Configure (or on an Apple Mac 
computer, from the wxMaxima menu, select Preferences). 


Alternatively, click the following icon on the toolbar. 
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This opens the ‘wxMaxima configuration’ window. 


(b) Within the Configuration window, click the icon labelled ‘Worksheet’ 
to access the worksheet settings, then tick the ‘Enter evaluates cells’ 
box by clicking it. 


This changes the behaviour so that pressing alone evaluates 


commands. 


(c) In the same list of settings, tick the ‘Open a cell when Maxima 
expects input’ box by clicking it. 


This ensures that the input prompt is shown whenever Maxima is 


waiting for input. 


(d) In the same list of settings, ensure the ‘Show user-defined labels 
instead of (%oxx)’ box is not ticked. 


This ensures that the output line numbers are always shown. 


(e) Click on the icon labelled ‘Options’ in the Configuration window, then 
tick the ‘Save sidebar status’ box by clicking it. 


This ensures that your preferred configuration of the wxMaxima 


window, for example hiding the Table of Contents, is retained between 
Maxima sessions. 


(£) Click ‘OK’ (or close the window if using an Apple Mac). 


Closing wxMaxima 


When you have finished using Maxima, you can close Maxima and the 
wxMaxima interface by using one of the following methods. 


e Select Exit from the File menu (or Quit wxMaxima from the wxMaxima 
menu on an Apple Mac computer). 


e Press Ctrl-Q. (Here, Q stands for quit.) 


e On a Windows computer, click the usual small cross button at the top 
right-hand corner of the wxMaxima window. 


You may be asked if you want to save your changes to the worksheet 
before closing. At this stage, click ‘No’. You will see how to save your 
work in Subsection 2.9. 
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2 Getting started with Maxima 


In this section you will learn about the basics of the wxMaxima interface 
and how to use it to perform simple calculations with Maxima. 


If you encounter unexpected problems when working through this section, 
remember to check the Frequently asked questions (FAQs) section on the 
OU Maxima website. 


2.1 Calculating with Maxima 


In Computer activity 3 you were asked to enter a numerical expression into 
Maxima using the symbols + for addition, * for multiplication and / for 
division. Maxima evaluated it using the rules for the correct order of 
mathematical operations; that is, the BIDMAS rules. The way in which a 
calculation or a command has to be entered in Maxima is known as its 
syntax. The Maxima syntax for the basic mathematical operations and 
other functions is summarised below. 


In this Guide we will use ^ for powers. Many Maxima commands take the 
form of a command name, such as sqrt or abs, followed by a pair of round 
brackets containing one or more objects, such as numbers, that the 
command operates on. Such an object is called an argument of the 
command. In the summary tables in this Guide, such as the one below, 
arguments of commands are often represented by 


Basic Mathematical operations 


Operation or function Syntax Example 
Addition + 4+3; 


Subtraction - 4-3; 
Multiplication * 4*3; 
Division / 4/3; 
Brackets (and ) 2* (3+4); 
Powers, for example 23 a 2^3; 

Or ** 2**3; 


Square root, for example V5 | sqrt(_) | sqrt(5); 
Magnitude (absolute value) | abs(_) abs (-3) ; 


You should end all Maxima commands with either a semicolon (;) or a 
dollar sign ($). Ending a command with a dollar sign tells Maxima to 
perform the instruction, but not to display the answer. When using 
wxMaxima, if you forget to finish your command with a semicolon or 
dollar sign, then the system will add a semicolon for you. Other interfaces 
may not do this. In such systems, pressing without ending the line 
with a semicolon or dollar sign simply moves the editing cursor to the next 
line, waiting for you to finish entering a complete command. 


You can also enter several commands on one line, ending each one with a 
semicolon or dollar sign. 


2 Getting started with Maxima 


Maxima generally ignores all spaces that you type while entering a 
calculation or command, so you can use spaces to make a command easier 
to read. Once wxMaxima has been configured to use to evaluate 
commands (as described in Computer activity 4), you can then use 
Shift-Enter to enter line breaks, which are also ignored by Maxima. 


Warning: Do not omit multiplication signs! 


A common error when first using Maxima is to omit multiplication 
signs. For example, you might type 


e 2sqrt(2) rather than 2*sqrt (2), or 
e 2(3+4) rather than 2*(3+4). 


Doing this will result in an error. The displayed error message might 
include one of the following statements: 


e parser: incomplete number; missing exponent? 


e incorrect syntax: Syntax error 
incorrect syntax: Too many )’s 
incorrect syntax: Premature termination of input at ;. 


e … is not an infix operator. 


If you receive one of these error messages, first check for missing 
multiplication signs! 


Warning: Take care with question marks (?) 


Do not type question marks when asking Maxima to evaluate an 
expression. The symbol ? is a special symbol in Maxima which 
provides access to underlying systems. 


The one exception to this is when you are using the commands for 
obtaining help from the system, which you will see in Subsection 2.10. 


You can calculate roots other than square roots, such as cube roots, by 
using the index laws. 


Computer activity 5 Calculating cube roots 


(a) Remembering that V27 = 273, calculate /27 by entering the 
following line into Maxima. 


27* (1/3) ; 


Remember to press to evaluate the command. 
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(b) What happens if you enter 
27/1/3; 


instead? Can you explain this? 


During the previous activity you may have noticed that when you enter an 
opening bracket, the wxMaxima interface automatically adds a closing one 
for you. This behaviour can be turned off by unticking ‘Match parenthesis 
in text controls’ in the ‘Worksheet’ settings of the Configuration window. 
(See Computer activity 4 for how to open the Configuration window.) 


2.2 wxMaxima cell markers 


Corresponding matching input and output lines in wxMaxima form a cell 
identified with the marker shown in Figure 5. 


Figure 5 The wxMaxima cell marker 


The top part of the cell marker (between the upper two horizontal marks) 
corresponds to the input line(s), and the bottom part corresponds to the 
output line(s). The cell marker turns red while an input line is being 
edited, and is shown surrounded by a box while the cell is being evaluated, 
as illustrated in Figure 6. 


Figure 6 The cell marker during evaluation 


Computer activity 6 Using cell markers 


(a) Click on the triangle at the top of the cell marker. 


This collapses the cell to a single line. This can be useful if a 
calculation or result is very long and so makes the worksheet difficult 


to read. 


The triangle at the top of the cell marker becomes filled. 


(b) Click the filled triangle at the top of the cell marker. 


2 Getting started with Maxima 


This expands the cell again. 


(c) Click part of the cell marker other than the top triangle to highlight 


it, then press the keyboard key. 


This deletes the cell from the worksheet. 


2.3 Decimal numbers 


Maxima usually tries to calculate quantities exactly rather than 
approximating them by decimal numbers. It does this by manipulating the 
symbols in the expression according to the rules of mathematics. This is 
sometimes known as symbolic computation. 


If you wish to have a result displayed as a decimal number, you can 
instruct Maxima to do so with the float command. The name of the 
command float arises from the fact that the method used by computers 
to store decimal numbers internally using the binary digits 0 and 1 is 
called a floating-point representation. The float command instructs 
Maxima to convert a number stored symbolically as a mathematical 
expression to one stored using a floating-point representation. 


Computer activity 7 Obtaining exact and approximate 
answers 


In Maxima, calculate 1309/3451 in three different ways, by entering each 
of the following expressions in turn. 


(a) 1309/3451; (b) 1309.0/3451; = (c)_ float (1309/3451) ; 


Maxima returns an exact answer for part (a). 


Since part (b) involves a decimal number, 1309.0, Maxima returns a 
decimal approximation to the answer. 


In part (c), the command float instructs Maxima to return a decimal 
answer. 


Converting to decimal numbers 


Convert to a decimal number | float( ) | float(sqrt(2)); 
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Maxima usually displays decimal numbers to 16 significant figures (though 
it suppresses trailing zeros at the end of the decimal part of a number; for 
example 0.125 is displayed simply as 0.125 rather than as 
0.1250000000000000). This is also the precision to which Maxima 
performs decimal calculations. You will see in Subsection 2.7 how to 
change the number of significant figures displayed. 


You may also see decimal numbers displayed in the form 
6.223015277861141 107% 


which is wxMaxima’s way of displaying scientific notation. It means 
6.223 015 277 861 141 x 10'8°. (The command-line interface would display 
this as 6.223015277861141E+139.) You can also enter numbers given in 
scientific notation. For example, to enter 6.022 x 10°, you can type 


6.022e23 


Computer activity 8 Simplifying and approximating 
numerical expressions 


For each of the following expressions, use Maxima to simplify it and then 
find a decimal approximation for it. 


2/3 
ov wl) om 


The expression in part (c) doesn’t involve decimal numbers and so the 
result is evaluated exactly. Notice, however, that the middle 91 digits 
of the answer are not displayed because the system does not expect 
you to read them all! 


The decimal form of the result, 3.273390607896142 101°, is given in 
wxMaxima’s way of displaying scientific notation. It means 
3.273 390 607 896 142 x 10150, 


Maxima knows the values of standard mathematical constants such as the 
irrational numbers e and zr. It treats them symbolically unless you 
explicitly ask for them in floating-point form. 


The constants e and 7 


e he he”“2; 
T Lpi Qx*hpi; 


As with the input and output line labels, the ‘%’ symbol indicates the 
name of a quantity built into Maxima. 


When displaying mathematics in an output line, wxMaxima displays Api 
as T. 
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Computer activity 9 Working with built-in constants 


T 
Use Maxima to calculate the value of the expression 7 + 3 + e°, both 


exactly and as a decimal number. 


Notice that in the exact answer, hpi is displayed as 7, but %e is 


unchanged. 


2.4 Troubleshooting wxMaxima 


If you are using wxMaxima and it does not seem to be responding, try the 
following suggestions. 


(Suggestions for troubleshooting Mazima when using the command-line 
interface are given in Subsection 14.8.) 


Troubleshooting wxMaxima 


Cell Marker Status message Comments 


Boxed: (| Maxima is calculating Maxima is evaluating a command that is 
taking a long time to complete. 


Boxed: í Ready for user input Maxima is waiting for you to type something 
before continuing. 


In either case you might wish to abort the 
command in one of the following ways. 


x 
e Click 8 on the toolbar. 


e Select Interrupt from the Maxima menu. 
e Type Ctrl-G (on an Apple Mac hold 


down the key while pressing C) J 


Maxima might take some time to respond to 
this request. 


Maxima process terminated | Maxima (but not the wxMaxima interface) has 
stopped working. 


Select Restart Maxima from the Maxima menu. 


If you encounter other problems, check the Frequently asked questions 
(FAQs) section on the OU Maxima website: 
learn1.open.ac.uk/site/maxima. 
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2.5 Reusing and editing commands 


In Computer activity 8 you may have entered each expression twice: once 
to evaluate it exactly and once for a decimal answer. This is, however, not 
necessary. When you have entered and evaluated an expression, you can 
edit it and then evaluate it again, as shown in the following activity. 


(Details on reusing and editing commands when using the command-line 
interface are given in Subsection 14.4.) 


Computer activity 10 Editing a command 


(a) Calculate 2/3 + 527 by entering the following line into Maxima. 
2*sqrt (3)+5*sqrt (27) ; 
(b) Edit your entry in part (a) to calculate 2/3 — 5/27 as follows. 


(i) Click on the expression you entered in part (a), or move to it 
using the up and down keyboard arrow keys. 


As you do this the cell marker will turn red and you will see a small 
flashing vertical line (the editing cursor) appear in the expression. 


(hili) 2*sqrt (3) +B¥sqrt (27); 


(ii) Edit the expression by using the (backspace) or 


key to delete the + sign and enter a - sign instead. 


(iii) Press . 


Pressing re-evaluates the cell. When a cell is re-evaluated in 


this way, the input and output line numbers are updated as if a new 
cell had been added. 


Should you click in the space between two cells, a horizontal cursor will be 
shown, as in Figure 7. This indicates the position at which a new cell will 
be created, if you begin typing here. This is also shown as you move 
between cells using the up and down keyboard arrow keys. 


(hil) *pi+ňpi/3+%e^2; 
AT y2 
(hot) = + he 


(%i2) float (Api+kpi/3+he“2) ; 
(402) 11.57784630371704 


Figure 7 ‘The wxMaxima horizontal cursor, positioned between two cells 


2 Getting started with Maxima 


In wxMaxima you can also copy and paste expressions, or parts of 
expressions, as you might copy and paste when using other software; that 
is, using Copy and Paste from the Edit menu, or the Ctrl-C and Ctr1-V 
keyboard shortcuts. 


Maxima also permits the result of one calculation to be used in another 
calculation. The symbol % represents the result of the last command 
evaluated. Due to the updating of line numbers when a command is 
re-evaluated, this might not be the last command in the worksheet, but it 
will be the command with the highest line number. So, for example, if you 
enter sqrt (%); Maxima calculates the square root of the last evaluated 
output. You can also use the output (or input) of other lines, by typing 
the line label in a calculation. For example, if you enter %405°3; Maxima 
calculates the cube of the expression in output line 5. 


Computer activity 11 Using and editing previous 
calculations with wxMaxima 


(a) Enter 6+2; into a cell at the bottom of your worksheet and evaluate it. 


(b) On the following line type 3+%; and enter it. 


Maxima returns 11, the sum of 3 and the previous answer. 


(c) Edit the line where you typed 6+2;, change it to 6*2; and re-evaluate 


the cell by pressing . 


Notice that the cell number and output of the line you edited is 
updated, but the result of 3+%; does not change. This is because 


when the 3+%; cell was evaluated, the symbol % represented the result 
of the old version of the edited command. 


(d) From the Cell menu, select the Evaluate All Visible Cells 
option, or press Ctrl-R. 


This re-evaluates all visible cells in the worksheet, in the order in 
which they appear, updating all the line numbers as it does so. 


Notice that the result of 3+%; has now changed. At the time it was 
re-evaluated the symbol % represented the value of the last evaluation, 
which was the revised version of the previous line. 


(e) In a new cell at the bottom of your worksheet enter %o +5; 
with | replaced by the output line number corresponding to the 
result of 6*2 entered earlier. 
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This command adds 5 to the result of the line referenced. 


Re-evaluating a wxMaxima worksheet 


Selecting the Evaluate All Visible Cells option from the Cell 
menu re-evaluates all visible cells in the worksheet, in the order in 
which they appear. 


This can also be achieved by pressing Ctr1-R. 


By default, the re-evaluation of a worksheet will stop if an error is 
encountered. This behaviour can be changed so that evaluation continues 
past the error by unticking ‘Abort evaluation on error’ in the ‘Maxima’ 
settings of the Configuration window. 


2.6 Variables 


Maxima allows you to assign a value to a variable, and then use that 
variable in further calculations. 


For example, you can assign the value 23 to the variable a, and then use a 
in subsequent expressions or commands, such as a? — 3a + 2. You can also 
assign expressions, both numerical and (as you will see later) algebraic, 
and indeed any other Maxima object, to a variable, and use that variable 
in subsequent commands. 


To assign a value or expression to a variable you use a colon (:). For 
example, the command 


a:23; 


assigns the value 23 to the variable a. The colon can be read as ‘is 
assigned the value’. 


Warning 


You cannot use = to assign a value to a variable. This symbol is used 
to define equations in Maxima. 


Variable names can be any combination of letters and numbers that begins 
with a letter. For example, Maxima will accept any of the following 
variable names. 


a A solution solution2 x2b 


Some words are built-in to Maxima and cannot be used as variable names. 
If you try to assign such a variable, you will get an error such as: 
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assignment: cannot assign to ... 
-- an error. To debug this try: debugmode(true) ; 


If you encounter such an error, choose a different name for your variable. 


Variable names are case-sensitive; for example, the variable x is different 
from the variable X. If a variable name is the name of a Greek letter, then 
wxMaxima will display it in output lines as the proper Greek character; for 
example, the variable alpha will be displayed as a. 


Once a variable has been assigned a value, Maxima will remember this 
value for the rest of your session. This may sometimes surprise you. It is 
easy to forget that earlier in the day you assigned a value to, say, x; you 
may then be confused later on when a calculation involving x gives an 
unexpected result. You can, however, tell Maxima to ‘forget’ about a 
variable to which you have previously assigned a value. You can also ask it 
to list all the variables that you have assigned values to. These commands 
are demonstrated in the following activity. 


Computer activity 12 Working with variables 


(a) Assign the value of v8 to the variable a by entering 
a:sqrt (8); 


The value of a is displayed as output, in an equivalent form. 


Note that when assigning variables, wxMaxima precedes the output 
with either the output line number or the name of the variable being 


assigned a value. Which is shown is determined by whether or not 
‘Show user-defined labels instead of (%oxx)’ is ticked in the 
‘Worksheet’ settings of the Configuration window. In this Guide, the 
output line numbers are shown. 


(b) Enter a; 


This displays the current value of a. 


(c) Assign the value of av2 to b. 


Don’t forget to include a multiplication sign when inputting av/2. 


The value of b is simplified when displayed. 


(d) Edit the line where you assigned /8 to a, so that the value /7 
becomes assigned to a. 


(e) Enter b; to display the value of b. 
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Notice that the value of b has not changed. The variable b was 


defined when a had its original value. 


(f£) Enter values; to list the names of all the variables that are currently 
assigned values. 


The variable names are given within square brackets, separated by 


commas. This is how Maxima displays lists. 


(g) Enter kill(a); to remove the variable a from the system. 


You should obtain the output done. 


(h) Enter a; to display the value of a. 


The variable name is output, as a no longer has a value. 


(i) Enter b; to display the value of b. 


b is still defined. 


Notice that in Computer activity 12(d) and (e), when the value of a was 
changed, the value of b was not affected, since it was defined when a had 
its original value. If you want to update the value of b using a different 
value of a, then you should enter a new value of a then re-evaluate the line 
defining b. You can do this in wxMaxima by selecting the line and pressing 
, or re-evaluating the whole worksheet by selecting the option 
Evaluate All Visible Cells from the Cell menu. 


Working with variables 


Operation Command Example 


Assign a value to a variable variable: value | a:23; 
Display the value of a variable | variable a; 

List all user-assigned variables | values values; 
Remove an assigned variable kill( variable) | kill(a); 
Remove all assigned variables | kill(all) kill(all) ; 


Note: here the placeholder variable represents any variable name. 


You can practise using variables in the following activity. 
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Computer activity 13 Using variables 


(a) Define the variable a to have the value 42. 
(b) Define the variable b to be equal to a? — 3a + 4. 


(c) Define the variable c to be the decimal approximation of the square 
root of b. 


2.7 System variables 


There are some variables built into Maxima whose values affect the 
behaviour of the system. These are called system variables. You met 
one system variable in Computer activity 12: the variable values holds a 
list of the names of all the variables you have defined. 


Another system variable, fpprintprec, specifies the number of significant 
figures of decimal numbers that are displayed. The name fpprintprec is 
an abbreviation of ‘floating-point print precision’. 


To change the system behaviour so that, for example, only 4 significant 
figures are displayed, use the command fpprintprec:4; to assign the 
value 4 to the variable fpprintprec. 


You can set the variable fpprintprec to any value between 2 and 16. 
Also, setting it to 0 restores the default behaviour of displaying 16 
significant figures. The value of fpprintprec sets the requested number of 
significant figures, but Maxima does not always exactly meet the request. 


Note that setting fpprintprec only affects how decimal numbers are 


displayed: it does not change the accuracy with which calculations are 
performed. 


Computer activity 14 Displaying a specified number of 
significant figures 


Display the decimal approximation of m to 8 significant figures by doing 
the following. 


(a) Set the system variable fpprintprec to be 8. 
(b) Display the decimal approximation of r. 


You can reset the values of fpprintprec and all other system variables to 
their original values by using the reset () command, which has no 
arguments. Try entering reset () ; followed by float (%pi); and check 
that the value of m is displayed to 16 significant figures. Note that the 
reset command also resets your Maxima line numbers to start at 1 again. 
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Resetting system variables 


Operation Command | Example 


Reset all system variables 


Note that the list of variables displayed when you enter values; does not 
include any system variables. If one of your variables does not appear in 
the list displayed by values;, then you have probably used a variable 
name that is also the name of a system variable and hence changed the 
value of that system variable. In extreme circumstances this may change 
the behaviour of Maxima, which can be restored by resetting all system 
variables as described above. 


2.8 Annotating your wxMaxima worksheet 


It is important to be able to clearly explain and present your mathematical 
work, both so that someone else can understand it, and so you can 
understand what you have done if you return to your work after some 
time. To help you do this, wxMaxima allows you to enter text comments 
within your worksheet, as illustrated in Figure 8. 


mw Maxi 16.04 y = 


File Edit View Cell Maxima Equations Algebra Calculus Simplify Plot Numeric Help 


[SA|SEXHOE Else: K2 


Circles 


V The radius, r, of the circle is defined below. 


The circumference, c, of the circle is calculated below. 


--> c: float (2*$%$pi*r); 
18.84955592153876 


The area, a, of the circle is calculated below. 


Pass a: float (%pit*r“2); 


28.27433388230814 
v 
|Welcome to wxMaxima Ready for user input 


Figure 8 Adding text to a wxMaxima worksheet 


To enter text into your worksheet, first you need to create a suitable cell in 
which to type. To do this, position the cursor and then select the 
appropriate one of the following commands from the Cell menu. 


e Insert Text cell 

e Insert Title cell 

e Insert Section cell 

e Insert Subsection cell 


The new cell will be created under the cell in which the editing cursor 
appears, or between two existing cells at the position indicated by the 
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horizontal cursor, if shown. Remember, the horizontal cursor appears 
whenever you click between two cells, or move within the worksheet using 
the up and down keyboard arrow keys. 


As you type your text, to start a new line within the cell press : 
When you have finished typing your text, click elsewhere on the worksheet 
or use the up or down keyboard arrow keys to move out of the text cell. 

If your text cell is at the bottom of the worksheet, and you click or move 
below it, an input prompt may not be shown until you start typing again. 


Section and Subsection cells are automatically numbered. 


Each Title cell, Section cell and Subsection cell has a small square in its 
top left corner. Clicking this square hides (or reveals) the contents of the 
worksheet from that point onward, until the next cell of the same type. 


Whenever you add a Section or Subsection cell, a corresponding entry is 
added to the Table of Contents window. You can display this window by 
selecting Table of Contents from the View menu, or by pressing 
Alt-Shift-T, that is, by holding down the (Alt) and (_ Shift ) keys while 


pressing : 


Computer activity 15 Including text in a worksheet 


(a) Create the worksheet shown in Figure 8. Use a Title cell for the title 
‘Circles’ and Text cells for the remaining text. Remember to click 
elsewhere on the worksheet or move position with the up or down 
keyboard arrow keys when you have finished entering each piece of 
text. 


After evaluating a Maxima command, you will need to correctly 
position the horizontal cursor before inserting the next text cell. 


(b) Use your worksheet to find the circumference and area of a circle of 
radius 5. (Remember to re-evaluate the worksheet, for example, by 
pressing Ctrl-R, after changing the value of r.) 


2.9 Saving and printing your work 


(Details on saving and printing your work when using the command-line 
interface are given in Subsection 14.5.) 


After working with Maxima you will probably want to save your 
calculations so that you can reuse or read them at a later date. You can do 
this in wxMaxima using the Save or Save As options from the File menu. 


There are various different formats in which you can save your work. 


It is recommended that you save your wxMaxima work as The input 
without images with the .wxm file type. This saves all your input 
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but not any calculated output. You have to recalculate the output 
when you re-open the file, as demonstrated in Computer activity 16. 


The alternative file type available is the Whole document with the .wxmx 
file type. This saves all the input and the output. Note that when you load 
a worksheet previously saved as a .wxmx file, although both the saved 
input and output are displayed, the input has not been evaluated by 
Maxima, so you cannot immediately use any results of the worksheet in 
new calculations. To re-evaluate the worksheet, and so recalculate all the 
results it contains, select Evaluate All Visible Cells from the Cell 
menu, or press Ctr1-R. 


Input commands can also be exported to a Maxima batch file, with the 
„mac file type, using the Export option of the File menu. This is for more 
experienced users of Maxima. 


You can open a saved wxMaxima worksheet by using the Open option from 
the File menu, or by double-clicking on the saved file on your computer. 


Computer activity 16 Saving your work 


Save your current worksheet then reload it, by following these steps. 
(a) Select Save As from the File menu. 
(b) In the window that appears: 
(i) Choose the folder in which to save your work. You might like to 
save it in the folder that you created in Computer activity 1. 
(ii) Choose a name for the file in which the worksheet will be saved. 
(iii) Make sure the file is saved as a The input without images with 
the .wxm file extension. 
(On a Windows computer, make sure that The input without 
images (*.wxm) is selected from the ‘Save as type’ drop-down 
list.) 
(iv) Click ‘Save’. 
(c) Close wxMaxima, by selecting Exit from the File menu, or pressing 
Ctr1-Q, or clicking the small ‘x’ at the top right of the window. (On a 
Mac select Quit wxMaxima from the wxMaxima menu or hold down the 


key while pressing (Q) .) 


(d) Open your worksheet in wxMaxima again, by either 


(i) double-clicking on the file you just saved, or 
(ii) following the steps below: 
e Start wxMaxima again. 
e Select Open from the File menu. Maxima will ask if you want 
to save your just opened, blank worksheet. Click ‘No’. 
e Find the file containing your work and click ‘Open’. 
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(e) Re-evaluate the worksheet by selecting Evaluate All Visible 
Cells from the Cell menu or pressing Ctrl-R. 


The calculations in the worksheet are evaluated and results displayed. 


Printing is best achieved by first exporting the whole worksheet to a 
HTML document (a web page) by selecting Export from the File menu, 
then opening the HTML document using a web browser and printing from 
that. The worksheet can also be exported as source code for the ATEX 
mathematical typesetting system using the same menu items. (Currently, 
these exports are not perfect representations of the worksheet.) 


Alternatively, you can obtain an image of the entire wxMaxima window by 
taking a ‘screenshot’. To do this on a Microsoft Windows computer, first 
click on the wxMaxima window, then press Alt-PrintScreen (that is, 
hold down the keyboard key while pressing the PrintScreen key, 
which may be labelled , or something similar). This stores 
the image, which you can then paste into a suitable application. If you are 
using Windows Vista, Windows 7, 8 or 10, you might like to use the 
‘Snipping Tool’, available by typing ‘snip’ into the Start menu. To take a 
screenshot on an Apple Mac computer, press Cmd-Ctr1-Shift-4 (that is, 
hold down the Command) , (Ctrl) and ED keyboard keys while pressing the 
key). Next press the space bar, move the camera pointer over the 
wxMaxima window and click. This stores the image, which you can then 
paste into a suitable application. You may also like to use the Grab 
application from the Utilities folder. 


You can also paste individual output lines (or parts of output lines) from 
the worksheet into, for example, word-processed documents. Select the 
part of the output that you wish to paste, then from the menu obtained by 
clicking the right mouse button, select one of the following options. 


e Copy — This copies the Maxima syntax for the selection. It can be useful 
if you wish to use the Maxima syntax in another document, or use the 
copied output as input to another Maxima command. 


e Copy as LaTeX — This copies a ATEX representation of the selection 
that can be pasted into a source file of the ATEX mathematical 
typesetting system. 


e Copy as MathML — This copies the selection as MathML, which is an 
international standard for the representation of mathematics. The 
MathML can be pasted into, for example, versions of Word from 2007 
onwards, where it will be pasted as mathematics that can be edited 
using the ribbon maths editor. It can also be imported into LibreOffice 
documents, by, within LibreOffice, first creating a maths object using 
the Insert > Object > Formula menu options, then importing using 
Tools > Import MathML from Clipboard. 
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e Copy as Image — This copies the selection as an image that can be 
pasted into another document. 


Whole cells can also be cut and pasted using the above options, by 
selecting the cell marker before clicking the right mouse button. However, 
when doing this, the TFX source contains a significant amount of 
formatting information that you may not require, and the MathML option 
for a whole cell does not work in the version of wxMaxima current at the 
time of writing. 


2.10 Getting help 


There are several ways in which you can obtain help with using Maxima, 
beyond this Guide. In wxMaxima, you can access the complete Maxima 
manual by selecting Maxima Help from the Help menu. This opens the 
window shown in Figure 9. You can use the Contents, Index and Search 
tabs to look for the information that you want. 


@ Help: Maxima 5.38.1.5.gdf93b7b Manua = Oo 
er SAV SS A 
Contents Index Search a 


[Top] [Contents] [Index] [2] 


Maxima 5.38.1.5.gdf93b7b 
Manual 


(bookmarks) ~| ĦA A 


$ wxMaxima manual 
m- Maxima manual 


Maxima is a computer algebra system, implemented in Lisp. 


Maxima is derived from the Macsyma system, developed at 

MIT in the years 1968 through 1982 as part of Project MAC. 

MIT turned over a copy of the Macsyma source code to the 
Department of Energy in 1982; that version is now known as 

DOE Macsyma. A copy of DOE Macsyma was maintained by 
Professor William F. Schelter of the University of Texas from 

1982 until his death in 2001. In 1998, Schelter obtained 

permission from the Department of Energy to release the ¥ 

< > 


Done 


Figure 9 The wxMaxima help system 


Alternatively, Maxima also includes commands that you can use to obtain 
help. These are demonstrated in the following activity and summarised in 
the table that follows. 


Computer activity 17 Getting help! 


(a) Enter 
? float; 


to show the help for the float command. 


Note the space between the ? and float, which is needed. This is one 
of the few times in Maxima when spaces matter! 
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The ‘?’ is the Maxima syntax to request help on the command given. 


The help information for the command may include more detail than 
you need! 


After giving the help information, Maxima displays the output true. 
This means that the command you entered was successful. 


(b) Suppose that you cannot remember the Maxima command for square 
root, but you do remember that it was something like sqr. 


Enter 
?? sqr; 


to list all the help information titles that contain the letter sequence 


6 bs 


sqr’. 
A list of possible titles is displayed. 


Title number 1 is what you were looking for, so enter 1; to display 
the relevant information. 


If more than one title is displayed in response to a ?? command, then 
Maxima will not continue until you enter one of the following: 


the appropriate number, to display your chosen information 


several numbers, separated by spaces, to see help information 
on several topics 


all to display the information on all topics listed, or 


none to see no information. 


A box will be shown around the cell marker while Maxima is waiting 
for your response. 


Help commands 


Operation Command Example 


Get help on a command ? command ? float; 
or describe( command ) describe(float) ; 
Find help information whose 2? ?? flo; 
title contains the given text | or describe( , inexact) | describe(flo, inexact); 


Note 1: the space after the ? is needed. 


Note 2: the second pair of commands displays a list of the titles of all 
help information pages whose title contains the given text. To obtain 
help on a particular topic listed, type the number corresponding to 
the required title in the list, followed by a semicolon. 
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There is further help available in the Frequently asked questions (FAQs) 
section of the OU Maxima website: learn1.open.ac.uk/site/maxima. 


3 Algebra 


In Section 2 of this Guide you used Maxima to perform numerical 
calculations. Much of the power of Maxima, however, lies in its ability to 
perform algebraic operations. In this section you will learn how to do this. 


3.1 Algebraic manipulations 


You can enter algebraic expressions into Maxima in a similar way to 
numerical expressions, using the symbols +, -, *, / and ^, and the 
command sqrt(_) for square roots. You can also assign algebraic 
expressions to variables, using the assign (:) command, in the same way as 
numerical expressions. 


In the next activity you will meet some commands for manipulating 
algebraic expressions. 


Computer activity 18 Using algebraic expressions 


(a) Enter 
expand( (2+x)*(3*x+4)%3 ); 


to expand the expression (2 + x)(3x + 4)°. 


The command expand multiplies out the expression that is its 
argument. 


Spaces have been included simply to help with the reading of the 
command. 


(b) Enter 
factor( 2*x^3-7*x^2-10*x+24 ); 


to factorise the expression 21° — 7x? — 10a + 24. 


The command factor factorises, where possible, the expression that 


is its argument. 


3 
x? +1 

Defi to be th i : 
(c) Define p to be the expression “ae 


Remember to use brackets as necessary when typing this expression. 
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(d) Enter 
fullratsimp(p) ; 


to simplify the expression assigned to the variable p. 


The command fullratsimp simplifies, where possible, the expression 
that is its argument. It does this by using techniques for multiplying 
out brackets, collecting terms over a common denominator and 
cancelling out common factors. Simplifications requiring other 


techniques may not be handled by fullratsimp in the way you would 
expect. 


The name fullratsimp is short for ‘full rational simplify’. 


(e) Enter 
subst (3,x,p); 


to find the value of the expression p (defined above) when x = 3. 


The command name subst is short for ‘substitute’. 


The command subst (3,x,p); substitutes the value 3 for the variable 
x in the expression p. 


The commands introduced in the previous activity are summarised below. 


Algebraic manipulation commands 


Operation Command Example 


Expand brackets | expand(_) expand( (x+1)^2 ); 

Factorise factor ($) factor( 2*x+4*x%2 ); 

Simplify fullratsimp( ) fullratsimp( (2*x+4*x^2)/x ); 

Substitute subst (value, variable, | subst(4, x, x^2+1); 
expression ) which substitutes 4 for x in x^2+1 


It is hard for a computer to simplify algebraic expressions: often there is 
no single ‘simplest form’ for a particular expression, and the computer has 
no way of knowing which form you prefer. You can use the commands in 
the table above to simplify an expression in the way you want. 


You can access many of these commands from wxMaxima menus. For 
example, to factorise an expression you can first enter the expression in the 
worksheet, select it by highlighting the text using the mouse or clicking on 
the corresponding output line number, then select Factor Expression 
from the Simplify menu (or from the menu obtained by right clicking on 
a highlighted output line number). 
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Other options on these menus, such as Substitute... open a window for 
you to enter the command options. Options that work in this way are 
identified by ... at the end of the option name. 


3.2 Equations and their solutions 


Maxima can manipulate and solve equations. To type an equation in 
Maxima, you use an equals sign in the usual way. For example, you can 
type 3*x+2=4. 


You can assign equations to variables in a similar way to expressions. 
For example, by entering 


q: 3*x+2=4; 
you assign the equation 3x + 2 = 4 to the variable q. 


You can extract the left- and right-hand sides of an equation by using the 
lhs and rhs commands, as demonstrated in the following activity. 


Computer activity 19 Working with equations 


(a) Assign the equation x? + 2x = 4 to the variable r. 
(b) Extract the left-hand side of the equation by entering 


lhs(r); 
(c) Extract the right-hand side of the equation by entering 


rhs(r); 


Extracting the left- and right-hand sides of an equation 


Left-hand side of an equation lhs( equation) | lhs(4*x+1=2*x-2) ; 
Right-hand side of an equation | rhs (equation) | rhs(4*x+1=24x-2); 


Note: equation indicates that the argument of the command must be 
an equation, or a variable whose value is an equation. 


Maxima can find the exact solution(s) of equations, as demonstrated in the 
following activity. 


Computer activity 20 Solving equations 


(a) Solve the equation x? + 2x = 1, by entering 


solve( x42+2*x=1 ); 


Notice that Maxima displays the two solutions of the quadratic 
equation within square brackets and separated with a comma: 


[x=-/2-1, x=/2-1] 


This is a Maxima list. 


(b) Use Maxima to solve 4x? — 127 + 9 = 0. 


Notice that although there is only one solution of this quadratic 
equation, Maxima still gives it as a list, containing just one element. 


(c) Use Maxima to solve 2z? + z? — 5x + 2 = 0. 


This is a cubic equation, since the highest power of x is 3. Each cubic 
equation has at most three solutions. 


(d) Try to use Maxima to solve x? + 2x — 4 = 0. 


Here, Maxima returns (within a list) the equation, slightly rearranged, 
but unsolved. 


The solve command always tries to find an exact solution of the 
equation using algebraic manipulations. If it is unable to do this, then 
it returns the original equation, as here. 


In Section 6 of this Guide you will see how to find approximate 
solutions of such equations. 


(e) Use Maxima to solve x? — 4x + 5 = 0. 


This equation has no real solutions, since the discriminant, 
(—4)? — 4 x 1 x 5 = —4, is negative. 


However, Maxima finds solutions that are complex numbers. These 
involve the imaginary number, usually denoted by 7, whose square 
is —1. Maxima denotes this number by %i. 


Complex numbers were mentioned in Unit 1 and you will learn more 
about them in Unit 12. 


The solve command is summarised below. 


Solving an equation 


Operation Command Example 


Solve an equation, exactly | solve( equation) | solve(2*x^2-1=0); 
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You saw in Computer activity 20 that the output of the solve command is 
a list. In Maxima, a list is a collection of individual elements separated 
by commas and enclosed within square brackets. You can learn more 
about using lists in the following activity. 


Computer activity 21 Working with lists 


(a) Assign the list containing the first five prime numbers to the variable 
A by entering 
A: [2,3,5,7,11]; 
(b) Display the fourth element of the list by entering 
A[4]; 


The position number of an element in a list is known as its index. 


You can obtain individual elements of a list by typing the index of the 


element in square brackets after the list, or after the name of a 
variable to which the list is assigned. 


(c) (i) Assign the solutions of the equation z? — x — 1 = 0 to the 
variable solns. 
(ii) Display decimal approximations to the answers by entering the 
command 
float(solns) ; 
(iii) Display the exact value of the first solution by entering the 


command 


rhs(solns[1]); 


The command solns[1] refers to the first element of the list solns 


Bel 
which is x=-———. This is an equation. The value needed is the 


right-hand side of this, which is extracted using the rhs command. 


The operations introduced in the previous activity are summarised below. 


Lists 


Form a list 1A EER A: [a,b,c]; 
Extract an element of a list | list [index] | A[2]; 


Note: ... indicates the list of arguments could continue. 


3 Algebra 


As well as solving equations, you can also manipulate equations using 
Maxima. The next activity shows how you can do this. 


Computer activity 22 Manipulating equations 


(a) Assign the equation 3y + 4 = 3x? + 9x to the variable eq. 


(b) Enter the following command. 


eq-4; 


This subtracts 4 from both sides of the equation eq. 


(c) Now divide both sides of the rearranged equation by 3, by entering 
1/3; 


Remember that % represents the result of the last evaluated cell. 


The variable y is now alone on the left-hand side of the equation. 


We have rearranged the equation to make y the subject. 


You can also use the solve command to change the subject of an 
equation: you solve the equation for the variable that you want to be the 
subject, as shown in the next activity. 


Computer activity 23 Rearranging an equation 


2a +b + 3c 
a—b 
(b) Rearrange the equation to make a the subject by entering 


(a) Assign the equation c = to the variable eqn. 


solve(eqn, a); 


Here, the second argument of the solve command is the variable that 


you want to solve for; that is, the variable that you want to be the 
subject of the equation. 


This use of the solve command is summarised below. 
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Changing the subject of an equation 


Operation Example 


Solve an equation for a variable | solve( equation, variable) | solve(2*a*b-3*b=0, a); 


Another use of the solve command is to solve two or more simultaneous 
equations. To use the command in this way, you must enter the equations 
as a list, and enter the variables to be solved for as another list, as 
indicated below. 


Solving simultaneous equations 


Solve simultaneous equations | solve( list of equations, | solve([2*x+y=4, 3*x-2*y=-1], 
list of variables ) [x,y]); 


Computer activity 24 Solving simultaneous equations 


Use Maxima to solve the simultaneous equations 
3x? + 2y = 20 
4x — 3y = —4 
as described below. 


(a) Assign the first equation to the variable eq1 and the second equation 
to the variable eq2. 


Maxima might reorder the terms of the equations when displaying the 


output. 


(b) Solve the pair of equations (eq1 and eq2) for the two unknowns 
(x and y) by entering 


solve( [eqi, eq2], [x,y] ); 


Maxima’s output is [[x=2, y=4], [x=- 2 > =- $1] 


This tells you that there are two solutions to the equations: one 


solution is x = 2, y = 4 and the other is x = —26/9, y = —68/27. 
Notice that Maxima’s output here is a list of lists. 
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3.3 Plotting graphs 


You can plot graphs in wxMaxima by using the wxplot2d command. This 
command inserts your graph into your worksheet, as demonstrated in the 
activity below. Maxima plots the graph by calculating a large number of 
points on the curve and joining them up. The ‘2d’ in the name of this 
command comes from the fact that its graphs are two-dimensional plots. 
The ‘wx’ is because it is a wxMaxima command. 


If you are using another interface to Maxima, you should use the command 
plot2d instead. You use it in the same way as wxplot2d, but it will plot 
the graph in a different computer window, which you then need to close 
before you can create another plot. 


Computer activity 25 Plotting a graph 


Plot the graph of the equation y = 4x? — 8x + 2 for values of x such that 
—1 < x < 2, as described below. 


(a) Enter the command 


wxplot2d(4*x*2-8*x+2, [x,-1,2]); 


The first argument of the command is the expression to be plotted, 
and the second argument is a list that specifies the variable in the 
expression together with the least and greatest values that you want it 
to take in the graph. 


wxMaxima displays the graph next to a line label beginning %t, such 
as ({t1). This is how wxMaxima shows intermediate output, that is, 
output other than the final result of the command. 


The wxplot2d command has no final result, just intermediate output, 
which is the graph. 


(b) Change the vertical range of the graph to be —5 < y < 5 by including 
an additional argument [y,-5,5]. 


That is, change the command to 


wxplot2d(4*x*2-8*x+2, [x,-1,2], ly,-5,5]); 


Notice the warning: 


plot2d: some values were clipped. 


This means that some of the points on the graph with z-coordinates 
between —1 and 2 could not be displayed due to the range of values of 
y specified. 


3 Algebra 


39 


Computer algebra guide 


40 


(c) Change the colour of the curve to green, by including the additional 
argument [color, green] after the other arguments of the wxplot2d 
command, but before the closing round bracket. Remember 
arguments are separated with commas, so include a comma before the 
new argument. 


This argument is a list whose first element is a keyword indicating 
the property to be set, in this case color (with the American 


spelling!). The second element is the value of the property that you 
want to choose. 


You can plot multiple graphs in the same diagram by including a list of 
expressions as the first argument of the command, as shown in the 
following activity. 


Computer activity 26 Plotting multiple graphs 


Suppose that the displacements sı and s2 (in km) of two cars along a 
straight road from a given starting point are given by sı = 15¢? + 10t and 
s2 = 60t, respectively, where t is the time (in hours) and 0 < t < 2. 


(a) Plot graphs of the displacement of the cars for 0 < t < 2 by entering 


wxplot2d([15*t*%3+10+#t, 60*t], [t,0,2]); 


Notice that the first argument of the command wxplot2d is a list 
containing the two expressions to be plotted. In this case, the 
expressions involve the variable t, so the second argument specifies the 
range of values of t to be plotted. 


The curves are automatically plotted in different colours, and a legend 
showing which expression corresponds to which curve is included at 
the top right-hand corner of the graph. 


(b) Change the colours of the curves by adding the following additional 
argument to the command used above: [color, green, magenta]. 


That is, change the command to 


wxplot2d([15*t*%3+10*t, 60*t], [t,0,2], [color, green, magenta] ); 


Maxima uses the first colour listed for the graph of the first expression 


listed and so on. If you do not include enough colours, Maxima starts 
using the listed colours again. 


3 Algebra 


(c) Now change the labelling within the legend to be the variable names 
given to the expressions in the question, rather than the expressions 
themselves. Do this by including another argument 


[legend, "si", "s2"] 


in a similar way to the colour argument; that is, within the round 
brackets of the wxplot2d command and separated from the other 
arguments by a comma. 


In general, to change the legend text, you include an extra argument 
of wxplot2d: a list that contains the keyword legend, followed by the 


new pieces of text, each enclosed in double quote marks. 


(d) Finally, change the labelling of the horizontal and vertical axes, by 
adding two more arguments to the command: 


[xlabel, "time, t (h)"] 
and 


[ylabel, "displacement, s (km)"]. 


In general, to change the x-axis label, you include a further argument 
of wxplot2d: a list that contains the keyword xlabel, followed by the 


new label text enclosed in double quote marks. You can change the 
y-axis label in a similar way, using the keyword ylabel. 


Changing the graph size 


Finally, the size of the plot within your worksheet can be changed using 
the wxplot_size system variable. For example, entering 


wxplot_size: [800,600] ; 


changes the size of subsequently plotted graphs from the default setting of 
600 x 400 pixels to 800 x 600 pixels. 


Summary 


The different ways in which you can use the wxplot2d command are 
summarised below, along with some of the optional arguments. If you are 
using a different interface to Maxima, you can use the plot2d command in 
the same ways. 
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Plotting graphs 


Operation Command Example 


Graph plotting options 


Option Argument Example 
Vertical range ly. B. HD ly, 0, 5] 


Curve colour 
Legend text 


Turn legend off 
Horizontal axis label 
Vertical axis label 


Change the size of the graph 


42 


Plot a graph wxplot2d( expression , wxplot2d(x*2, [x,0,1]); 


horizontal range,...) which plots the graph of x? 


Plot several graphs | wxplot2d( list of expressions, | wxplot2d([x*2, 2*x], [x,0,1]); 
horizontal range ,...) which plots graphs of x? and 2x 


[color, patel (color, red] 
(legend, " ", ...] | (legend, "Car A"] 
(a label should be listed for each curve plotted) 
[legend, false] [legend, false] 
[xlabel, " "] [xlabel, "time, t (h)"] 
[ylabel, " "] [ylabel, "displacement, s (km)"] 


wxplot_size:[ , | wxplot_size: [800,600] ; 


for x between 0 and 1 


wxplot2d(x*2,[x,0,1],[y,0,2]); 
which plots the graph of x? 

for x between 0 and 1, and with 
vertical axis values between 0 and 2 


for x between 0 and 1 


Note: for a full list of the possible options, view the Maxima help 
obtained using the ?? plotting options; command. 


If you began studying Sections 2 and 3 of this Guide from Activity 31 of 
Unit 2, this completes your study of Unit 2. 


4 Functions 


You can define a function in Maxima using the := operator, as shown in 
the following activity. 


Computer activity 27 Defining functions 


(a) Define the function f in Maxima, with the rule f(x) = 2° — 3x? + 4, 
by entering the following command. (Use the characters : followed 
by = to type the := operator.) 


f (x) :=x^5-3*x^2+4; 
(b) Calculate the value of f(2), by entering f (2) ; 


(c) Now define the function g with rule 


x-1 
JS 5 


and find the value of g(3). 


Remember to use brackets appropriately in the rule for g. 


(d) Now try to use Maxima to find g(2). 


This gives an error. 


This is because when Maxima tries to evaluate g(2), it finds that this 


involves dividing by zero, so it returns an error message. The function 
g is undefined at x = 2. 


(e) Enter h(1); 


Maxima returns the expression h(1), unevaluated. 


Since the function h has not been defined, Maxima cannot calculate a 
value. It knows only that the value is h(1). 


(f) Enter functions; 


This displays the value of the system variable functions, which is a 


list of all the functions that you have defined in this Maxima session. 


Notice the difference between defining a function, which has an argument 
and which can be evaluated at different values of the argument, using :=, 
and assigning an expression to a variable (which has no argument) using :. 
As with variables, once you have defined a function in Maxima it is 
remembered by the system for the rest of your session, unless you remove 
it using the kill command that you met in Subsection 2.6 of this Guide. 


Working with functions 


Operation Command 


Define a function ‘= f(x) :=2*x+3; 


Evaluate a function at a value | function (_) f); 
List all user-defined functions | functions functions; 
Remove a defined function kill(function) | kill(f); 


4 Functions 
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You can plot the graph of a function using the wxplot2d command 
introduced in Subsection 3.3 of this Guide, as demonstrated in the next 
activity. 


Computer activity 28 Plotting functions 


(a) Plot the function f defined in Computer activity 27 over the range 
—2 < x < 2, by entering 


wxplot2d(f(x), [x,-2,2]); 


Note that the first argument must be f (x) rather than just f. 


(b) Plot the function g defined in Computer activity 27 over the range 
—5<a<5. 


Notice the vertical scale of the graph, and the large spike near x = 2. 


—1 
The function g has the rule g(x) = 2 57 50 it is not defined at x = 2. 
i bf — 


Maxima has evaluated it at values of x close to 2, obtained large 
values and attempted to join up the corresponding points. 


To obtain a better graph, you should restrict the range of the vertical 
axis. You are asked to do this next. 


(c) Plot g again, for values of x in the range —5 < x < 5, but restricting 
the range of the vertical axis to between —5 and 5. 


Remember, to restrict the vertical range to —5 < y < 5, include the 


argument [y,-5,5] in the wxplot2d command. 


If you began studying this section from Activity 11 of Unit 3, return to the 
unit to continue your study of functions. 


5 Exponential functions and 
logarithms 


You saw in Subsection 2.3 of this Guide that the mathematical constant e 
is represented in Maxima by %e. The % symbol indicates that this is the 
name of a special quantity built into Maxima. 


5 Exponential functions and logarithms 


Computer activity 29 Working with exponentials and 
logarithms 


(a) Find e° using Maxima, by entering 


hed ; 


Since Maxima works symbolically, it does not return a decimal 


approximation. 


(b) Find a decimal approximation to e° by entering 
float (%e^5); 
(c) An alternative way to find powers of e is to use the exp function. 
Find e° by entering 
exp (5) ; 


Again, this is expressed exactly, as a power of e. 


(d) The Maxima command for a natural logarithm is log(_). 
Use this to find a decimal approximation for In 8. 


(e) Use the solve command introduced in Subsection 3.2 of this Guide to 
solve the equation 


3°7 = 5. 


You may be surprised to see two solutions to this equation. 
log(5) 
2log(3) 
log(—v5) 
log(3) 


One solution, ‚is the sort of solution that you would expect. 


The other solution, , involves the logarithm of a negative 


number. 


You have seen that only positive numbers have logarithms, but that 
applies when you are working only with real numbers. If you are 
working with complex numbers, then negative numbers do have 
logarithms. 


log(—v5) 
log(3) 

clearly by obtaining a decimal approximation for it. It contains the 

imaginary number i, which Maxima represents as %i. You will learn 


The solution is a complex number, as you can see more 


about complex numbers in Unit 12. 


There is only one real solution to the equation. 


45 


Computer algebra guide 


46 


The Maxima commands for the exponential function and logarithms are 
summarised below. Maxima has no command 1n(_); remember that the 
command for a natural logarithm is log(_). 


Exponentials and logarithms 


Exponential, for example e? 


Natural logarithm, In 


When you are working with logarithms in Maxima, the algebraic 
simplification commands that you met in Section 3 of this Guide can be 
useful. You can also use the commands shown below. 


Simplifying expressions involving exponentials and logarithms 


radcan (log (x/2)) ; 


Operation 


Simplify something radcan( ) 
involving exponentials and 
logarithms 


Combine logarithms 


logcontract( ) | logcontract (log(a)+log(b)) ; 


The name of the radcan command arises from the fact that it converts 
expressions involving radicals, that is, roots, powers and logarithms, into a 
canonical (or standard) form. 


These commands are demonstrated in the following activity. 
Computer activity 30 Simplifying expressions involving 
logarithms 
Simplify the expressions 
ln(z? + 2x? + x) — ln(z) 
ln(x + 1) 
by following these steps. 


and In(a) + 2 In(b) — In(c) 


(a) Define u to be the expression 
ln(z? + 2x? + x) — ln(z) 
ln(x + 1) 


Remember that the Maxima command for a natural logarithm is log. 


(b) Simplify u using the following command. 


radcan (u); 


Can you simplify this expression and obtain the same answer by hand? 


(c) Define v to be the expression In(a) + 2ln(b) — In(c). 
(d) Simplify v by using the following command. 


logcontract (v); 


If you began studying this section from Activity 48 of Unit 3, return to the 
unit to continue your study of functions. 


6 Trigonometry 


The Maxima commands for the trigonometric functions are given below. 
Like many computer systems, Maxima assumes that all angles are 
measured in radians. Notice that the names of the inverse trigonometric 
functions sint, cos~! and tan~! in Maxima are asin, acos and atan, 
respectively. This notation is used by many computer systems; it is short 
for arcsin, arccos and arctan, the alternative names often used for sin” *, 


cos! and tan~!. 


Trigonometric functions 


Function | Syntax Example 


sin sin( ) sin(1); 

cos cos ($) cos (3*%pi/2); 
tan tan( ) tan(%pi/4) ; 

cosec csc (i) csc(2); 

sec sec( ) sec (%pi/3); 

cot cot ($) cot (%pi/4); 

sin * asin( ) | asin(sqrt(3)/2); 
cos”! acos( ) | acos(1/sqrt(2)); 
tan! atan( ) | atan(1/2); 


Computer activity 31 Using trigonometric functions 
Calculate the following using Maxima. 


(a) cos (=) 


Remember that to enter the constant m in Maxima, you type “pi. 


Maxima returns an exact answer, in the form of a surd. Remember 


that Maxima always uses exact values where possible, rather than 
decimal approximations. 


6 Trigonometry 
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(b) sin 45° 


T 
180° 


Remember to convert the angle to radians, by multiplying 45° by 


(c) sec (3.4) 


Here, since the angle is entered as a decimal number, Maxima returns 
a decimal approximation to the answer. 


(e) cosec (=) 


There is no simple exact answer here, so Maxima leaves this 
expression unevaluated. 


To find a decimal approximation, you can use the command 


float( csc(%pi/7) ); 


Again, Maxima leaves this expression unevaluated, but it simplifies it 
slightly, using the fact that tan(—@) = — tan 8. 


(g) cos”! (2.0) 


Here, since the argument of the function is a decimal number, 
Maxima returns a decimal approximation to the answer. 


You might be surprised that it returns an answer at all, since 2.0 is 


not in the domain of cos~!. 


Notice that the answer is a multiple of %i. In Maxima, %i represents 


the imaginary number i whose square is —1, which you will learn 
about in Unit 12. Although there is no real number y such that 
cosy = 2, there is an imaginary number with this property, and 
Maxima returns this value. 


If you are working only with real numbers, you should ignore any 
Maxima answers that contain the constant %i. 


(h) cos! (2) 


This expression is equivalent to the expression in part (g), but does 
not contain a decimal number, so Maxima tries to find an exact 
answer. Since there is no simple exact answer, Maxima leaves the 


expression unevaluated. 


This hides the fact that the answer is not a real number. You need to 
take care when using inverse trigonometric functions. 


You can solve equations involving trigonometric functions in Maxima, 
using the solve command that you met in Subsection 3.2 of this Guide. 
However, Maxima will return only one solution to such an equation, even 
though many may exist. It warns you of this, as shown in the following 
activity. 


Computer activity 32 Solving trigonometric equations 


(a) Solve the equation cos(x) — 0.6 = 0 by entering 


solve(cos(x)-0.6=0) ; 


Since the solve command always tries to find exact solutions, 
Maxima first converts the decimal number —0.6 to the exact 
fraction —2. It tells you that it has done this. 


It also displays the following warning to tell you that there might be 
other solutions. 


solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


Then it returns a single solution. 


(b) Find a decimal approximation to the solution obtained in part (a), by 
entering 


float (%) 


Remember, % represents the answer of the last evaluated calculation. 


(c) Solve the equation sin(3x + 4) = 4. 
(d) Solve the equation cos(2x + 4) — sin(x/3) = 0. 
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Maxima returns the original equation, slightly rearranged but 


unsolved. This is because the solve command could not find an exact 
solution to this equation. 


In circumstances such as Computer activity 32(d) when the solve 
command fails to find the exact solution of an equation, you can use 
another Maxima command, find_root, to try to find an approximate 
solution. 


Whereas the solve command manipulates the equation mathematically to 
find a solution, the find_root command performs a series of numerical 
calculations to find a decimal approximation to a solution. This is known 
as solving the equation numerically, or finding a numerical 
approximation to a solution. 


You know from Unit 2 that the solutions of an equation of the form 
f(x) = 0 are the x-coordinates of the points at which the graph of 
y = f(x) crosses the x-axis, as illustrated in Figure 10. 


Xy 


Figure 10 The graph of y = f(x) with the solutions of f(x) = 0 marked 


The find-root command uses the following property of functions. If f is a 
function and [a,b] is an interval included in its domain such that one of 
f(a) and f(b) is negative and the other is positive, then there is at least 
one solution of the equation f(x) = 0 in the interval (a,b). Figure 11 
illustrates this property for two different functions. This property holds as 
long as the graph of f is continuous on the interval [a, bl; informally, this 
means that you can draw the graph of f for input values in the interval 

(a, b] without taking your pen off the paper. 


The find_root command looks for a single solution of an equation of the 
form f(x) = 0 that lies within an interval that you specify. The command 
requires that the values of f at each endpoint of the interval have different 
signs. If this is not the case, then Maxima gives an error message. 
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Figure 11 There is at least one solution of f(x) = 0 in the interval (a,b) if 
f(a) and f(b) have different signs 


Solving equations numerically 


Solve an equation numerically find_root (expression , variable , find_root (2*x*3-1,x,0,1); 


interval start value, 
interval end value ) 


Note 1: the command finds a value of the variable at which the 
expression is zero. 


Note 2: the values that the expression takes at the endpoints of the 
interval must have different signs. 


Computer activity 33 Solving equations numerically 


(a) Find a numerical approximation to a solution of the equation 
cos(2z + 4) — sin (=) = 0, 


as follows. 


(i) First, define the function f in Maxima, with the rule 


f(x) = cos(2x + 4) — sin (5) ; 


Remember, to define a function use :=. 


(ii) Next, plot the graph of f(x) over the interval —2 < x < 2, say. 
This will help you to find an interval that contains a solution. 


Remember, you can plot graphs using the wxplot2d command. 


From the graph, you can see that f(0) < 0 and f(1) > 0, so there is a 
solution in the interval (0,1). 
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If there were no solution in the interval over which you plotted the 


function, then you could plot the function over a larger interval, and 
so on, until you find a suitable interval, if one exists. 


(iii) Find a numerical approximation to the solution of the equation 
in the interval (0,1) by entering 


find_root (f(x), x, 0, 1); 


Remember to type the underscore symbol (_) in the middle of the 
command name. 


Notice that the first argument of the find_root command is an 
expression (or a function defined by an expression) not an equation. 
Maxima finds an input value at which the expression takes the value 0. 


(b) Using similar methods, find a solution of the equation 
sin?(x +2) = 3a. 


Hint: consider the function g(x) = sin?(a + 2) — 3x and find a value of 
x for which g(x) = 0. 


Remember that sin?(x + 2) means (sin(x + 2))?. You can enter this 
expression in Maxima by typing either (sin(x+2) )%2 or sin(x+2)%2. 


Maxima can use trigonometric identities to rewrite expressions involving 
trigonometric functions in different forms. 


These behave like the algebraic simplification commands you met in 
Section 3 of this Guide. In particular, you need to guide Maxima in how to 
conduct such rearrangements. Commands for doing this are demonstrated 
in Computer activities 34 and 35. 


Computer activity 34 Manipulating trigonometric 
expressions 


(a) Enter 


trigreduce( sin(x)’2 ); 


The trigreduce command attempts to express powers of sines and 


cosines of x in terms of sines and cosines of multiples of x. It reduces 
the powers of the functions. 
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(b) Enter 
trigexpand(%) ; 


Remember that % represents the result of the previous calculation. 


The trigexpand command attempts to express trigonometric 


functions of sums, differences and multiples of angles in terms of 
trigonometric functions of the individual angles. 


(c) Enter 


trigsimp(%) ; 


The trigsimp command attempts to simplify expressions involving 
trigonometric functions. 


This returns you to the expression that you started with! 


Computer activity 35 Simplifying rational trigonometric 
expressions 


S 


n(4z) 


i ; 
by entering 
sin(x) 


Simplify 


trigrat( sin(4*x)/sin(x) ); 


The trigrat command attempts to simplify algebraic fractions 


involving trigonometric functions. 


Trigonometric manipulation commands 


Operation Command Example 
Expand trigonometric functions trigexpand( ) | trigexpand(sin(A+B)); 


of sums, differences and 
multiples of angles 
Express powers of sines and trigreduce( ) | trigreduce(sin(x)%2); 
cosines in terms of sines and 
cosines of multiple angles 


Simplify trigonometric expressions | trigsimp( ) trigsimp(sin(x)2+cos(x)%2) ; 
Simplify algebraic fractions trigrat(_ ) trigrat ((sin(2x))/cos(x)); 
containing trigonometric 
functions 


You can use these commands, and others that you learned previously, in 
the final activity of this section. 
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Computer activity 36 Finding a new trigonometric identity 


(a) Find an expression for cos(5x) in terms of products of sin(x) and 
cos(x). 

(b) Plot the graphs of cos(5x) and the expression you found in (a), to 
check that they look the same. 


If you began studying this section from Activity 37 of Unit 4, this completes 
your study of Unit 4. 


7 Plotting circles 


In this section you will learn how to use Maxima to plot circles. 


In Subsection 3.3 of this Guide, you learned how to use the wxplot2d 
command to plot curves. You can use this command to plot a curve only if 
the equation of the curve is in explicit form; that is, only if it is written in 
the form y = f(x), where f is a function. 


To plot a curve represented by an equation in implicit form, such as 
x? +y? = 1, you need to use the wximplicit_plot command. This is 
demonstrated in Computer activity 37. 


The wximplicit_plot command is not part of the basic Maxima collection 
of commands, so you have to add it by loading an additional Maxima 
package, the implicit_plot package, which was installed on your 
computer when you installed Maxima. The Maxima command to load an 
additional package is load, which is also demonstrated in Computer 
activity 37. Each package needs to be loaded only once per Maxima 
session. 


If you are using a Maxima interface other than wxMaxima, then you 
should use the implicit_plot command instead of the wximplicit_plot 
command. You load and use it in exactly the same way as 
wximplicit_plot. 


Computer activity 37 Plotting a circle 


Plot the circle represented by the equation x? + y? = 1, as follows. 
(a) Load the implicit_plot package by entering 


load(implicit_plot) ; 


Remember to type the underscore symbol (_) in the middle of the 


package name. 


This package provides the wximplicit_plot command. 


(b) Plot the circle by entering 


wximplicit_plot (x‘2+y*%2=1, [x,-2,2], [y,-2,2]); 


Again, remember to type the underscore symbol (_) in the middle of 
the command name. 


The first argument of the command is the equation to be plotted. 
Notice that in Maxima, equations are given using the usual equals 
sign. 

The second argument specifies the variable to be plotted on the 
horizontal axis, and the range of values of that variable to be plotted. 


The third argument specifies the variable (and its range) to be plotted 
on the vertical axis. You must include this third argument. 


The circle with equation x? + y? = 1 has centre (0,0) and radius 1. So 
the circle lies within —2 < x < 2 and —2 < y < 2 and hence these 
seem suitable ranges to use. 


The curve that Maxima plots, however, does not look very circular! 
This is discussed below. 


Note that the first argument of the wximplicit_plot command is an 
equation, whereas the first argument of the wxplot2d command is an 
expression (or a function or variable representing an expression). If an 
expression is given as the first argument, wximplicit_plot will plot the 
graph of the equation obtained by setting the expression equal to zero. 


The curve that you plotted in Computer activity 37 does not look circular 
because it has different scales on the axes. Whatever ranges of xz- and 
y-values you specify for a plot, by default Maxima makes the y-axis appear 
3 as long as the x-axis. 


To specify the use of equal scales, you can add the following additional 
argument to the wximplicit_plot command: 


same_xy 


Note that in some older releases of Maxima, the option needed to 
specify equal scalings on coordinate axes was 


[gnuplot_preamble, "set size ratio -1"] 


If the same_xy option does not change the scales of your graph, then 
use this option instead. 
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Computer activity 38 Plotting a circular circle! 


Plot the circle z? + y? = 1, using equal scales on the z- and y-axes by 
entering 


wximplicit_plot (x*2+y*%2=1, [x,-2,2], [y,-2,2], same_xy); 


Including the same_xy (or gnuplot_preamble) argument in the 
wximplicit_plot command can be tiresome if you are plotting lots of 
curves and need equally scaled z- and y-axes for each. An alternative is to 
change the default behaviour of Maxima so that it uses axes with equal 
scales for the rest of your session. You will see how to do this in the 
following activity. 


Computer activity 39 Plotting another circular circle! 


Plot the circle (x — 1)? + y? = 3, using equal scales on the z- and y-axes as 
follows. 


(a) Tell Maxima to always use equal scales on the z- and y-axes by 
entering 


set_plot_option(same_xy) $ 


This sets same_xy to be a default option for all plots during the rest of 
your Maxima session. The line ends with $ to prevent the display of 
the output of the command. Without this, the output consists of a 
long list of all the current default settings, which can safely be ignored. 


If you are using an older release of Maxima that uses the 
gnuplot_preamble option, use 


set_plot_option([gnuplot_preamble, "set size ratio -1"])$ 


instead. 


(b) Plot the circle by entering 
wximplicit_plot ((x-1)*2+y%2=3, [x,-1,3], [y,-2,2]); 


If you wish to stop plotting graphs with equal scales and restore the 
default Maxima plotting behaviour, then set the same_xy plotting option 
to be false using 


set_plot_option([same_xy, false])$ 


If you are using an older version of Maxima, the command needed to reset 
the default view is 


set_plot_option([gnuplot_preamble, ""])$ 


7 Plotting circles 


To plot more than one curve at once you can list their equations within 
square brackets as the first argument of wximplicit_plot, as you did with 
the wxplot2d command. The equations do not necessarily have to be in 
implicit form. This is demonstrated in the following activity. 


Computer activity 40 Plotting a line and a circle 


Plot the line y = 3x — 5 and the circle (x — 2)? + (y — 3)? = 2 as follows. 


(a) For simplicity, assign the equation of the line to the variable line, 
and the equation of the circle to the variable circle. 


Remember, to assign an expression to a variable, use the colon (:) 
operator. So to assign the equation of the line to the variable line, 


use 


line: y=3*x-5; 


(b) Plot the line and circle on axes with equal scales by entering 


wximplicit_plot([line, circle], [x,0,5], [y,0,5]); 


The first argument is a list of the equations to be plotted. Since you 
have assigned the equations to variables, you can simply list the 
names of the variables. 


Since the circle has centre (2,3) and radius V2 ~ 1.41, the ranges 
O0<a<5and0<y< 5 are large enough for the whole circle to be 
plotted. 


The commands introduced in this section are summarised below. 


Loading additional packages 


Load a package | load( package name) | load(implicit-_plot); 


Plotting curves represented by equations in implicit form 


Operation Command Example 

Plot a curve represented | wximplicit_plot( equation , | wximplicit_plot(x2ty*%2=1, 
by an equation in horizontal range, Bells pst, A03 
implicit form vertical range, ...) 

Plot several curves wximplicit_plot( wximplicit_plot( 
represented by list of equations , [x^2+2*y^2=4, y=x^3], 
equations in implicit horizontal range , [x,“1,1], [ys-1,1)3 
form vertical range, ...) 
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Note 1: the implicit_plot package needs to be loaded before this 
command can be used. 


Note 2: the first argument is an equation or a list of equations. 
Alternatively it can be an expression, or a list of expressions, in which 
case the equation formed by setting the expression equal to zero is plotted. 


Note 3: equations in explicit form can also be plotted using this command. 


Graph plotting options 


op i e 


Set axes to have equal scales same_xy same_xy 


(in older versions of Maxima) | [gnuplot_preamble, [gnuplot_preamble, 
"set size ratio -1"] "set size ratio -1"] 


Setting default plotting options 


Set the default plotting options | set_plot_option( option) | set_plot_option(same_xy) 


Circles are not the only type of curve that can be represented by equations 
in implicit form. In the next activity you will plot some other equations of 
this type and see the curves that they represent. 


Computer activity 41 Plotting equations given in implicit 
form 
Plot the curves represented by the following equations, using equal axis 


scales for each curve. Make sure you use suitable ranges of x and y. 


Remember, you can use the editing facilities of wxMaxima to avoid 


having to type the whole wximplicit_plot command each time. 


(a) a +4y?=1 (b) a? —4y?=1 (©) yr =a etl 
(d) grats (e) gba (tryin 


If you began studying this section from Activity 16 of Unit 5, return to the 
unit to continue your study. 
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8 Differentiation 


The Maxima command for differentiating an expression is diff, which is 
demonstrated in the following activity. This command finds derivatives by 
using the derivatives of standard functions together with the sum rule, 
constant multiple rule, product rule, quotient rule and chain rule. 


Sometimes Maxima may give the result in a different, but equivalent, form 
to the expression that you might obtain when differentiating by hand. If 
this happens, then you might like to check that the two answers are 
equivalent. You can do this by finding the difference between the two 
answers and checking that it simplifies to zero. A rougher check is to use 
Maxima to plot graphs of the two answers and check they appear to be the 
same. 


Computer activity 42 Demonstrating differentiation with 
Maxima 


(a) Find the derivative of f(x) = zt + 5x? + 15, by entering 
diff (x^4+5*x^2+15, x); 
The first argument of the diff command is the expression to be 


differentiated, and the second argument is the variable to differentiate 
with respect to. 


You must include the second argument. If you do not, then the 


answer returned by Maxima will be the derivative multiplied by 
del(x). This notation is outside the scope of MST124. 


As you would expect, the first argument can also be a variable that 
you have previously defined to be an expression, or a function you 
have previously defined. 


£2 
(b) Define the variable p to be equal to sind)’ and then find its 
sin 
derivative with respect to t by entering 
diff (p‚t); 


(Remember, to assign a value to a variable you use a colon (:).) 


Using diff in this way can help you check the expression to be 


differentiated has been entered correctly before you differentiate it. 


(c) Define the function g(u) = (u? +3) In(u?), and find its derivative using 
diff (glu) ,u); 


(Remember, to define a function you use :=, and the Maxima 
command for natural logarithms is log.) 
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Notice that the first argument given here is g(u), not just g. 


(d) Find the second derivative of the expression p defined above, by using 
diff (p,t,2); 


The third argument of diff specifies the number of times that the 


expression is to be differentiated. 


(e) Find the third derivative of g(u) defined above. 


The diff command is summarised below. You can practise using it in the 
next activity. 


Differentiation 
Operation Command Example 
Differentiate diff( , variable) diff (sin(x*2) ,x); 
Differentiate diff(_, variable, diff (log(x) ,x,3); 
multiple times positive integer ) 


Computer activity 43 Differentiating with Maxima 


Use Maxima to differentiate the following, simplifying the answer where 
possible. 


Remember, you can enter e” in Maxima either as %e^x or exp (x). 


(a) 20° + = (b) In(e?+2) (©) = 


You can use the diff command to assign the derivative of a function to 
another named function. This is useful when you want to evaluate the 
derivative at a particular input value, for example. 


The next activity shows you how to do this — unfortunately it’s not 
completely straightforward. In the activity, the original function is denoted 
by f and its derivative is denoted by df, rather than by f’. This is 
because you cannot use f’ as a function name in Maxima. 


Computer activity 44 Assigning a derivative to a function 


sin(x) 


(a) Define the function f(x) = : 
et +1 
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(b) Try to define df (x) to be the derivative of f(x), by entering 
df (x) :=diff (f(x) ,x); 


(c) Try to find the value of the derivative of f at x = 3 by evaluating 
df (x) at x = 3; that is, try entering 


df (3); 


This gives an error. 


The cause of this error lies in the previous command: 


df (x) :=diff (f(x) ,x); 


The problem is that when you enter this command, Maxima does not 
perform the differentiation. Instead, it just notes that df(x) means 
diff (f(x) ,x). Then, when you enter df (3), Maxima understands 


this to mean diff ( (3) ,3), which is meaningless, as you cannot 
differentiate with respect to a particular number. 


To force Maxima to perform the differentiation, you can use the 
syntax ''()) (which contains two single quotation marks, not a 
double quotation mark, and a pair of round brackets). This is 
demonstrated below. You need to use this syntax whenever you want 
to assign the result of a Maxima command to a function. 


(d) Define df(x) to be the result of differentiating f with respect to x, by 
entering 


df (x) :='' (diff (f (x) ,x)); 


Remember that '' is two separate quotation marks. Use the upright 
quotation mark on your keyboard ('), not the sloping one (`). Don’t 
forget to include the outer round brackets too. 


Notice that this time the output line shows that df (x) is defined to 
be the result of differentiating f(x). 


(e) Find the value of the derivative of f at x = 3 by entering 
df (3); 


If you want a decimal approximation to the result, use the float 


command. 


(f) Usea similar approach to find the derivative of sin()e®®) at x = 1, 
giving your answer as a decimal number to 3 significant figures. 
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The command used to force the evaluation of the derivative in Computer 
activity 44 is summarised below. 


Forcing the evaluation of a command 


Operation Command | Example 


Force the evaluation of a g(x) :='' (diff (x%4,x)); 
command 


An alternative approach to Computer activity 44 parts (d) and (e) is to 
assign the derivative of f to a variable rather than use the derivative to 
define a function. You can do this by entering, for example 


df : diff (f(x) ,x); 
This calculates an expression for the derivative of f at x and assigns it to 
the variable df (the ''( ) command is not required). 
You can then evaluate the derivative at 3 by substituting the value 3 for zx, 
which you can do by entering 


subst (3,x,df) ; 


If you began studying this section from Activity 22 of Unit 7, return to the 
unit to continue your study. 


9 Integration 


The Maxima command to integrate an expression is integrate. This is 
demonstrated in the following activity. 


As with differentiation, the results that you obtain from integrating using 
Maxima might be in a different algebraic form to the expression that you 
might obtain when integrating by hand. To show that the two answers are 
equivalent integrals, you could find the difference between the two answers 
and check that it simplifies to a constant. A rougher check is to use 
Maxima to plot graphs of the two answers (taking the arbitrary constant 
to be zero, for example) and check they seem to be vertical translations of 
each other. 


Computer activity 45 Demonstrating integration with 
Maxima 


(a) Evaluate i 2x? + e?” dz by entering 


integrate (2*x*2texp(3*x), x); 


The first argument of integrate is the expression to integrate, and 


the second is the variable to integrate with respect to. 


9 Integration 


The output returned by Maxima does not include an arbitrary 
constant. That is, Maxima gives an antiderivative of the expression, 


not its indefinite integral. When you write down an indefinite integral 
that you have found using Maxima, you need to remember to add an 
arbitrary constant yourself. 


3 
1 
(b) Evaluate the definite integral | — dx by entering 
2 T 


integrate(1/x, x, 2, 3); 


The third and fourth arguments of integrate are the lower and 
upper limits of integration, respectively. 


Notice that, as usual, Maxima returns an exact result. Use the float 
command to find a decimal approximation, if required. 


(c) Define the function p(x) = , and define q(x) to be an 


x 
xr? +1 
antiderivative of p. 
If you try to define q(x) using 
q(x) :=integrate(p(x),x); 


then Maxima does not actually do the integration, but just notes that 


q(x) means fræ dz. 


As with differentiation, to force Maxima to perform the integration 
you need to use the ''(__) command and enter 


q(x) :='' (integrate (p(x) ,x)); 


Integration 


Operation Command Example 
Integrate integrate( , variable) integrate (x/\2,x); 


(find an antiderivative) 


Evaluate a definite integral | integrate( , variable, integrate(sin(x) ,x,0,1); 
lower limit , 


upper limit ) 


You can practise using the integrate command in the following activity. 
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Computer activity 46 Using integration 


Use Maxima to find the following. 


bl [(e-3e- Dee (b) ks (c) f estar 


It is not always possible to express the result of integrating a function in 
terms of the mathematical functions you have previously met. The 
following activity illustrates some of the things that you may encounter 
when integrating functions using Maxima. 


Computer activity 47 Calculating harder integrals 


Try to calculate the following using Maxima. 


(a) | ene) de (b) | ae 10 I ar 


Maxima was unable to find an antiderivative for part (a), and so 
returned the integral uncalculated. In parts (b) and (c) answers were 
returned containing the function erf, which is known as the Error 
function, and the function gamma_incomplete, the incomplete gamma 


function. These are examples of what are known as special functions 
in mathematics. They cannot be simply expressed in terms of 
functions you have already met, but many can be defined in terms of 
antiderivatives of simpler functions. 


Sometimes you may want to integrate a function that contains a variable 
other than the one you are integrating with respect to. For example, the 


integral I x” dx contains the variable n as well as the variable 7. When 


you use Maxima to find such an integral, you sometimes need to provide 
information about the other variable, as demonstrated in the following 
activity. 


Computer activity 48 Giving Maxima more information 


(a) Try to use Maxima to find the integral J x” dz, as follows. 


(i) First, ensure that n has no predefined value, by entering 
kill(n); 
(ii) Then, enter 


integrate(x“n, x); 


9 Integration 


Maxima asks the question: Is n equal to -1? 


This is because the answer to J x” dx depends on the value of n. 


1 
If n = —1 (so n +1 = 0), the integral is J — dx, which is equal to 
£ 


In |æ| + c, whereas if n is any other value, then the integral is equal 
grt 


to Ji +c. So Maxima cannot calculate the integral without 
n 
knowing more about n. 


Enter the answer n; 
(which is short for no) 


You could also answer the question with y; to mean yes. 


Maxima displays the expected result. 


(b) An alternative approach is to give Maxima information about n before 
you use the integrate command. You can do this as follows. 


(i) Tell Maxima that n # —1 by entering 


assume (notequal (n,-1)); 


The assume command tells Maxima information about variables. For 
example, assume (a>0) ; tells Maxima that a is positive. 


Maxima has no simple symbol for ‘not equal to’, so you have to use 
the notequal command for this, as above. 


(ii) Now integrate x” using 


integrate(x“n, x); 


No question is asked this time, since Maxima knows all it needs to 
know about n. 


(iii) Enter 
facts(); 


to see all the additional facts about variables known to Maxima. 


This lists all such facts. The only one shown should be the one you 


gave earlier. 


(iv) Tell Maxima to forget about the assumption on n by entering 


forget (notequal(n,-1)); 
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(v) Type facts(); again to see all the facts now known. 


No facts are displayed, since Maxima has now forgotten the 


assumption about n. 


The commands introduced in the previous activity are summarised below, 
together with the equal command, which is used in a similar fashion to 
notequal. 


Facts about variables 


Operation Command Example 


Make an assumption about a assume(_ ) assume (a>0) ; 

variable 
State two things are equal equal( , ) assume (equal (n,3)) ; 
State two things are not equal | notequal( , ) assume (notequal(n,-1)); 
Forget a property forget (property) | forget(a>0); 
List all known facts facts() facts(); 
List all known facts about facts (variable ) facts(a); 

a particular variable 


Computer activity 49 Using facts 


Calculate the integral 


1 
|= dz, where a > 0. 
T+ a 


You saw in Computer activity 47 that Maxima cannot always find an 
antiderivative for a function. Even if Maxima cannot find an antiderivative 
of a function f, it can still find an approximate value for a definite integral 


b 
of the form f f(x)dx. 


In Unit 8, you saw how the value of the definite integral can be 
approximated by breaking the area between the graph of the function and 
the x-axis into a number of rectangles. Maxima uses a more sophisticated 
version of this basic method to calculate approximate values of definite 
integrals. 


You will learn the Maxima command for this in the following activity. 
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Computer activity 50 Approximate values of definite 
integrals 


1 
(a) Try to calculate f In(cos(a)) dx by entering 
0 


integrate(log(cos(x*2)), x, 0, 1); 


The integral is returned unevaluated, since Maxima cannot 


calculate it. 


(b) Find an approximate value of the integral by entering 


quad_qags(log(cos(x*2)), x, 0, 1); 


The quad_qags command calculates an approximate value for a 
definite integral. 


The ‘quad’ part of the name of this command arises from the fact the 
process of finding the value of a definite integral (or any area) is 
sometimes called quadrature. The letters q, a, g and s in the second 
half of the name specify the particular method Maxima uses to find 
the approximate value of the definite integral. This method gives 
good results for a wide variety of integrals. 


Notice that the output of the quad_gags command is a list. The first 
element of the list is the approximate value of the definite integral, the 
second element is an estimate of the accuracy of the approximation, 
the third element is the number of values at which the function to be 
integrated was evaluated during the calculation, and the final element 
is an error code. An error code of 0 means that no errors occurred. 


1 
(c) Find an approximate value for the definite integral Í e” dt. 
0 


Finding approximate values for definite integrals 


Find an approximate value | quad-gags(_ „variable, quad_qags (exp (-x%2) ,x,0,1); 
of a definite integral lower limit , upper limit ) 


If you began studying this section from Activity 47 of Unit 8, return to the 
unit to continue your study. 


Integration 
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10 Matrices 


In this section you will learn how to use Maxima for matrices. 


10.1 Matrix algebra 


There are several ways to input matrices in Maxima. Two of these are 
demonstrated in the following activity. 


Computer activity 51 Entering matrices 


2 —3 3 
(a) Assign the matrix | —2 1 2} to the variable P by entering 
1 -l1 4 


P:matrix( [2,-3,3], [-2,1,2], [1,-1,4] ); 


The matrix command allows a matrix to be defined. You give each 
row of the matrix as a list (that is, within square brackets with the 
elements separated by commas). 


By default, wxMaxima displays matrices using square brackets: 
2 3 3 
-2 1 2 
1 -1 4 


If you are using a different interface to Maxima, or have chosen 
different display settings, matrices may be displayed as follows. 


2 -3 3 ] 


[ 
[ 
[ 
[ 
[ 


(b) If you are using wxMaxima, you can also input a matrix using an 
on-screen form. Use this form to assign the matrix 


5 —3 
=i =2 
0 4 


to the variable Q as follows. (If you are not using wxMaxima, then 
enter this matrix using the matrix command, as in part (a).) 


(i) From the Algebra menu select Enter Matrix.... 
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Make sure you select Enter Matrix..., not Generate Matrix.... 


A window opens requesting properties of the matrix, as shown below. 


Matrix 


(ii) The matrix for this example has size 3 x 2, so enter 3 as the 
number of rows and 2 as the number of columns, replacing the 
default values shown. 


It is also possible to give information on the type of matrix, but for 
now leave this set to general. 


Finally, enter Q as the variable name to which the matrix will be 
assigned, and click ‘OK’. 


This opens a second window requesting the elements of the matrix. 


Enter matrix 


(iii) Complete entry of the matrix by filling in the elements in this 
second window. You may like to use the Tab key ((—)) to move 
the cursor between the different elements. When you have 
finished, click ‘OK’. 
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Completion of this on-screen form has the effect of including a matrix 
input line in your worksheet, followed by the appropriate output. 


(c) Assign the matrix 4 7 0) to the variable R using an 
-5 1 8 
appropriate method. 


The typed command for entering a matrix is summarised below. 


Entering matrices 


Note: matrices can also be entered using the Enter Matrix... option 
of the Algebra menu. 


Once a matrix has been input, you can obtain its size and individual 
elements by using the commands demonstrated in the following activity. 


Computer activity 52 Using matrices 


(a) Find the size of the matrix Q defined in Computer activity 51 by 
entering 


matrix_size(Q); 


The output is a list containing the number of rows and the number of 


columns of the matrix. 


(b) Extract the element in the third row and second column of R by 
entering 


R[3,2]; 


Alternatively, R[3] [2] gives the same result. 


Notice the use of square brackets to reference elements of a matrix is 


the same syntax as used in Section 3.2 of this Guide to obtain 
individual elements of a list. 


(c) Extract the third row of R by entering 
R[3]; 
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(d) Change the element in the third row and second column of R to be 0 
by entering 


R[3,21:0; 
then display the revised matrix by entering 


R; 


You can use commands of this type to correct any errors that you 


might make when entering a matrix. 


Using matrices 


Operation Command Example 
Size of a matrix matrix_size (matrix ) matrix_size(A) ; 
Row of a matrix matrix [row] A[2]; 
Element of a matrix | matrix [row , column] A[3,4]; 
or matrix [row] [column] | A[3] [4]; 


The syntax to use for matrix operations in Maxima is given below. You 
can practise using them in Computer activity 53. 


Matrix operations 


Operation Syntax | Example 
Matrix addition + A+B; 


where A and B are matrices of the same size 


Matrix subtraction = A-B; 

where A and B are matrices of the same size 
Scalar multiplication | * 2*A; 
Matrix multiplication | . A.B; 

where A and B are matrices of appropriate sizes 
Matrix powers OP AS 


where A is a square matrix 


Warning 


Remember to use . (not *) for matrix multiplication and ^^ (not ^) to 
find the power of a matrix. 


If you use the symbol * when multiplying two matrices or the symbol 
^ to find the power of a matrix, then Maxima may return an answer, 
but in general it will be the wrong answer. (It will be an answer found 
using a different matrix operation.) 


10 Matrices 
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Computer activity 53 Matrix operations 


(a) Define the following matrices in Maxima 


2 0 0 3 13 
p-(3 4). a-(5 4), re 9 
—4 6 1 3 


(b) Find P+Q. 


The matrices P and Q are the same size, so this calculation is 


possible. 


(c) Find P-R. 


The matrices P and R are not the same size, so this calculation is not 
possible. Maxima gives an error message: 


fullmap: arguments must have same formal structure. 
-- an error. To debug this try: debugmode(true) ; 


(d) Find 2P. 
This is a scalar multiple of a matrix, so use the * operator. 
(e) Find PR. 


This is matrix multiplication, so use the . operator. 


The sizes of these matrices means the calculation is possible. 


(f) Find PQ. 


Remember to use the . operator. 


This calculation is not possible due to the sizes of the matrices. 
Maxima gives an error message: 


MULTIPLYMATRICES: attempt to multiply nonconformable 
matrices. 
-- an error. To debug this try: debugmode(true) ; 
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(g) Find PS. 


The matrices P and S have sizes 3 x 2 and 2 x 1, respectively, so this 


is possible and the result is a 3 x 1 matrix, that is, a column vector. 


(h) Find PT. 


The matrices P and T have sizes 3 x 2 and 1 x 2, respectively, so this 
product should not be possible. However, Maxima does return an 
answer! 


Maxima does not distinguish between single row matrices (which are 
sometimes called row vectors) and single column matrices (column 
vectors). It considers them interchangeable as circumstances require. 
You should take care to multiply only appropriately sized matrices. 


(i) Find R?. 


Since R is a matrix use the matrix power operator ^^. 


(j) Enter 
R‘2; 


The result is not the same as in part (i). 


Here the ^ symbol was used and Maxima squared each individual 
element of R. 


(k) Enter 
P*Q; 


The result is not PQ, since this product does not exist! 


Here the * symbol was used and Maxima multiplied corresponding 
elements of P and Q together. 


If you began studying this subsection from Activity 13 of Unit 9, return to 
the unit to continue your study. 
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10.2 Determinants and inverses 


The Maxima commands for finding the determinant and inverse of a 
matrix are demonstrated in the following activity. 


Computer activity 54 Finding determinants and matrix 


inverses 
Consider the matrices 
1 2 0 1 2 0 
A=l|l-—4 4 2 and B= {5 6 2 
—-1 2 1 2.2 1 


(a) First, enter the matrices A and B into Maxima. 


(b) Find the determinant of A by entering 


determinant (A); 


The determinant command calculates the determinant of a square 


matrix. 


(c) Find the inverse of A, and assign the result to the variable C by 
entering 


C:invert (A); 


The invert command calculates the inverse of a square matrix. 


An alternative way to specify the inverse of the matrix A in Maxima is 
A^^(-1), which accords with the notation usually used to denote a 
matrix inverse. 


(d) Check that C really is the inverse of A by calculating the products AC 
and CA. 


Remember to use the symbol . for matrix multiplication. 


Both products should be the 3 x 3 identity matrix. 


(e) Try to find the inverse of B. 


This gives an error: 


expt: undefined: 0 to a negative exponent. 


-- an error. To debug this try: debugmode(true) ; 


This matrix does not have an inverse. You can check this by 
calculating its determinant. 
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Another command that is sometimes useful when you are working with 
matrices is ident, which specifies an identity matrix. For example, 


ident (2) gives | ih the 2 x 2 identity matrix. 


1 
O 1 
Matrix operations 


Operation Example 


Determinant (of a square matrix) determinant (matrix) | determinant (A); 


Matrix inverse (of a square matrix) | invert (matrix) invert (A); 
or matrix ^^(-1) A^^(-1); 
Identity matrix ident (size) ident (2) ; 


Computer activity 55 Finding determinants and inverses of 
larger matrices 


For each of the following matrices, check whether it is invertible. If it is, 
find the inverse. 


-1 0 3 -2 
2 —1 1 it op 6 
(a) [2 O1 (b) B dt oe 
4 21 T 
-1 2 1 -1 


If you began studying this subsection from Activity 24 of Unit 9, return to 
the unit to continue your study. 


11 Sequences and series 


In this section you will learn how to use Maxima to work with sequences 
and series. 


11.1 Plotting graphs of sequences 


In Subsection 3.3 of this Guide, you learned how to use the wxplot2d 
command to plot the graphs of algebraic expressions. Here, you will see 
how to use the same command to plot graphs of sequences. First, in 
Computer activities 56 and 57, you will learn how to define and work with 
sequences in Maxima. 


Computer activity 56 Working with a sequence specified by 
a closed form 


Consider the sequence (an) specified by the following closed form and 
range of values of n: 


üs=08 (mE dak 
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(a) Define this sequence in Maxima, by entering 


a[n] :=0.8%n; 


This command defines the general term an of the sequence using the 
same define operator (:=) as is used when defining functions. 


Notice that you type the subscript n of the general term within square 
brackets, but the Maxima output uses a properly formatted subscript. 


(b) Find the term a100 of the sequence by entering 


a[100]; 


Notice that when you defined the sequence in Computer activity 56 you 
did not need to enter the range of values of n. In fact, Maxima places no 
restriction on the values of n that you can use when calculating terms in a 
sequence defined by a closed form. For instance, for the sequence (an) 
defined in Computer activity 56 it is possible to calculate values for 
a[3.5] and a[-4] using Maxima. 

However, these numbers are not terms of the sequence (an), which are 
defined only for integer values of n greater than or equal to 1. 


Computer activity 57 Working with a sequence specified by 
a recurrence relation 


Consider the sequence (bn) specified by the recurrence relation 
mel, be = Dyer) (nd ddr 
(a) Define this sequence in Maxima by entering the two commands below. 


bli :1; 
b[n] :=2*b[n-1] +1; 


Notice that you assign the value of the first term, b[1], using the : 
operator that is used to assign values to variables, but you specify the 
recurrence relation using the := operator that is used to define 
functions. 


Unlike for a sequence defined using a closed form, Maxima can 
calculate a term of a sequence defined using a recurrence relation only 
when the term number is an integer greater than or equal to the 
starting value, which in this case is 1. 


(b) Find the 5th and 100th terms of this sequence. 
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(c) Change the value of bj to be 10, by entering 
b[1]:10; 


and try to find the revised value of bs. 


Notice that the value of b5 displayed is unchanged. This is because 
when Maxima uses sequences, it remembers all the values of the 
sequence that it has previously calculated, and bases subsequent 
calculations on them. This improves the efficiency of the calculations. 


However, it also means that sequences cannot be easily changed. 


To change the definition of a sequence, you first need to remove the 
sequence from Maxima’s memory using the kill command, and then 
redefine it. 


(d) Revise the sequence (bn) correctly, by first entering 
kill(b); 
and then entering the two commands 
b[1] :10; 
bin] :=2*b[n-1] +1; 
to specify the revised sequence. 


Now calculate the value of bs. 


The commands for defining and working with sequences are summarised 
below. 


Sequences 


Operation Command Example 
Define a sequence sequence [n]:= expression in n a[n] :=n^2; 


using a closed form 


Define a sequence sequence [first term number]: initial value | p[1]:1; 

using a recurrence system | Sequence [n] := expression in sequence[n-1] bin] :=2*b[n-1]; 
Calculate a term in sequence [term number ] b [50]; 

a sequence 


Note: to change the definition of a sequence, you first have to delete it 
using kill. 


In the following activity you will learn how to plot a graph of a sequence. 
Remember, a graph of a sequence (an) is a set of points (n, an) plotted for 
a suitable range of values of n. 
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Computer activity 58 Plotting the graph of a simple 
sequence 
Plot the graph of the sequence (cn) specified by 
cn (9 =1,2,3) 
as follows. 
(a) First, define the sequence in Maxima by entering 
c[n] :=n^2; 


(b) Create a list of the coordinates of the points to be plotted; that is, the 
points (1,c1), (2,c2) and (3,c3) and assign it to a variable, say, pts 
(for ‘points’), by entering 


pte: [ [1,c[1]], [2,c[2]], [3,c[3]] 1s 


Notice that you type the coordinates of each point as a list in square 
brackets. The variable pts is a list of lists. 


Do not use the variable name points for your list of points — this is a 
word used by Maxima to set plot styles, as you will see next. 


(c) Plot the graph of the sequence using the command 


wxplot2d([discrete, pts], [style, points], 
[xlabel, "n"], [ylabel, "c[n]"]); 


Here, instead of the first argument of the wxplot2d command being 
an expression, or a list of expressions to plot, it is a list containing two 
elements. The first element is the keyword discrete. This tells 
Maxima that a set of discrete points are to be plotted. The second 
element is the list of the coordinates of the points to be plotted. 


The second argument of the wxplot2d command is 
[style, points]. 
This tells Maxima to simply plot points, and not join them with lines. 


The final two arguments specify labels for the x- and y-axes. These 
arguments were introduced in Subsection 3.3 of this Guide. 


You can also use other standard arguments of wxplot2d that you met 
earlier. For example, including [x,0,4] changes the range of the 
horizontal axis to 0 < x < 4. 
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The use of wxplot2d for plotting sequences is summarised below. 


Plotting graphs of sequences 


Plot a sequence | wxplot2d([discrete, list of coordinate pairs], | wxplot2d([discrete, 


[style, points], ...) [{1,5], (2,10]]], 
[style, points]); 


Graph plotting options 


Plot points | [style, points] | [style, points] 


In Computer activity 58 the sequence plotted was short, so it was 
relatively easy to create the list of coordinates to plot. For a longer 
sequence this can be laborious, but Maxima can help by creating the list 
for you. This is demonstrated in the next activity. 


Computer activity 59 Plotting the graph of a longer 
sequence 


Plot the graph of the sequence given by the recurrence relation 
di=0; dn =}4(2—dn-1) (W=2,3,4,..2,15) 

as follows. 

(a) First define the sequence (dn) in Maxima. 


(b) To plot the sequence, we need to plot the points 
(Lea) (2, d2), (8de). ..., (15, di5). 


Create a list of these points, and assign the list to the variable pts, 
using the command 


pts:makelist([n,d[n]], n, 1, 15); 


The makelist command creates a list. The first argument of the 
command is the general form of the elements of the list. In this case, 
each element has the form [n,d[n]]; that is, it is itself a list 
containing two elements: the term number and the value of the 
corresponding term. These are the coordinates of a point to be 


plotted. Here, the general form is expressed using the dummy 
variable n. 


The second argument of the makelist command is the name of the 
dummy variable used, and the third and fourth arguments are the 
minimum and maximum values taken by the dummy variable. 
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Creating lists 


When you create a long list you may not want to see all the output! 
Remember that you can end a command with a dollar sign ($) rather 


than a semicolon (;) to suppress the display of the output. 


(c) Plot the points in the list pts by using the command 


wxplot2d([discrete, pts], [style, points], 
[xlabel, "n"], [ylabel, "d[n]"]); 


as seen earlier. 


The makelist command used in the previous activity is summarised below. 


Create a list | makelist( general term, dummy variable, | makelist(2*n,n,1,100); 
start value , end value) 
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As you know from your study of Unit 10, a sequence can have infinitely 
many terms; such a sequence is called an infinite sequence. Since 
Maxima cannot plot an infinite number of points, when plotting a graph of 
an infinite sequence, you need to restrict the number of points to a large, 
but finite, number. Try to choose a number that is large enough to reveal 
the long-term behaviour of the sequence. 


Computer activity 60 Plotting part of some infinite 
sequences 


Use Maxima to plot graphs of the sequences specified below, by restricting 
the range of values of n to be n = 1,2,3,...,50 in each case. 


How do you think each sequence will behave as more and more terms are 
considered? 


(a) r1=1, Tn = l1.5rn-1 — 0.1(rn-1)? (n= 2,3,4,...) 
(b) sn =n? -10 (n=1,2,3,...) 
(c) tn=1-—1.1” (n=1,2,3,...) 


If you began studying this subsection from Activity 17 of Unit 10, return to 
the unit to continue your study, and in particular learn about the long-term 
behaviour of sequences. 


11.2 Summing series 


You can calculate the sums of various series using the Maxima command 
sum which is demonstrated in the following activity. 


Computer activity 61 Summing series 


(a) Find the sum of the squares of the first 100 integers, that is 
100 


Dn, 
n=1 
by entering 
sum(n*2, n, 1, 100); 


The first argument of the sum command is the general term of the 
series to be summed, which in this case is expressed in terms of the 


index variable n. The second argument is the index variable itself and 
the third and fourth arguments are the lower and upper limits of the 
summation. 


(b) Write the sum 


6) +6) +6) 


in sigma notation, and then use the Maxima sum command to 
calculate it. 
A 1 
(c) Try to calculate the sum of the infinite series De T by entering 
n=1 


sum(1/(4%n), n, 1, inf); 


Maxima uses inf to represent oo. 


Notice that Maxima displays the sum (unevaluated) using sigma 
notation. By default, Maxima does not attempt to calculate the sums 
of infinite series. 


You can, however, tell Maxima to simplify the result of a summation, 
which can sometimes result in it finding a value for a sum of an 
infinite series, if it exists. To do this, you change the value of the 


system variable simpsum, as shown below. (simpsum is short for 


‘simplify sums’.) 


(d) Set the system variable simpsum to have the value true by entering 


simpsum: true; 


oo 
1 
Now try to calculate `> Je again. 


n=1 
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Maxima now returns the sum of the series. 


< 1 
n=1 


The commands introduced in Computer activity 61 are summarised below. 


Summing series 


Operation Command Example 


Sum a series sum( general term, index variable, | sum(n‘2,n,1,10); 
lower limit , upper limit ) 


Simplify summations | simpsum: true simpsum: true; 


Infinity 


Constant | Syntax | Example 


oo inf sum(1/n“2,n,1,inf); 


If you began studying this subsection from Activity 32 of Unit 10, return to 
the unit to continue your study. 


12 Taylor polynomials 


Taylor polynomials is an optional topic and you will not be assessed on this 
section. 


You could use Maxima to find a Taylor polynomial for a function by using 
the diff command to calculate the required derivatives and then 
combining these to form the required polynomial. However, Maxima also 
has a command for calculating Taylor polynomials more efficiently. This is 
demonstrated in the following activity. 


Computer activity 62 Finding Taylor polynomials 


(a) Find the cubic Taylor polynomial about 7/6 for the function 
f(x) = cosa, by entering 


taylor(cos(x), x, %pi/6, 3); 


The taylor command calculates a Taylor polynomial for the function 


that is its first argument. Its second argument specifies the input 
variable of the function. 
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The third argument of the command is the point about which the 
Taylor polynomial is to be found, which is also known as the centre of 
the polynomial. The final argument of the command is the degree of 
the Taylor polynomial required. 


Notice that the Maxima output ends with ..., which indicates that 
more terms could have been found. Also, the output line number is 
followed by /T/. This indicates that the output is a (truncated) 
Taylor series. 


(b) Define p to be the quintic Taylor polynomial about 0 for the function 
f(x) = tan z, by entering 


p(x):=''(taylor(tan(x), x, 0, 5)); 


Here, the ''( ) command is needed to force Maxima to calculate the 
Taylor polynomial which is then used to define p(x). 


Without the ''( ) command, p(x) would be given by the expression 
taylor(tan(x), x, 0, 5). So if you tried to find, for example, 
p(3), then Maxima would interpret this to mean the expression 
taylor(tan(3), 3, 0, 5); but this is meaningless and would 
therefore produce an error. 


You first met the ''( ) command in Section 8 of this Guide. 


(c) Express the Taylor polynomial p(x) found above in terms of the 
variable y by typing 


p(y); 


As soon as a Taylor polynomial is used in a further calculation or 
command, the trailing ... is dropped, and the Taylor polynomial is 
converted to a rational expression. This is indicated by the marker 
/R/ following the output line number. 


(d) Find the decimal value of p(x) when x = 0.1. Also find the difference 
between p(x) and tanz at x = 0.1. 


Maxima warns you it has replaced certain decimal numbers with 


rational approximations when calculating these values. 


83 


Computer algebra guide 


84 


The taylor command is summarised below. 


Taylor polynomials 


Find a Taylor | taylor (function , variable, | taylor(sin(x) ,x,0,3); 
polynomial centre , degree ) 


Computer activity 63 Plotting Taylor polynomials 


Consider the function f(x) = e%™(@), 


Plot a graph of the function and the graph of its quintic Taylor polynomial 
about x = 1 on the same axes over the range 0 < x < 2. 


(Remember, you learnt how to plot more than one function on the same 
graph in Subsection 3.3 of this Guide.) 


If you began studying this section from Activity 21 of Unit 11, return to the 
unit to continue your study. 


13 Complex numbers 


In this section you will learn how to use Maxima to work with complex 
numbers. 


13.1 Working with complex numbers 


The imaginary number 7 is represented in Maxima by %i. You saw this in 
Subsection 3.2 of this Guide, when you solved equations that have complex 
solutions as well as real solutions. The usual arithmetic operations work 
with complex numbers in Maxima, and there are commands for finding the 
real and imaginary parts of a complex number, its conjugate, its modulus 
and principal argument, as demonstrated in the following activity. 


Computer activity 64 Using complex numbers 


(a) Assign the complex number 3 + 4% to the variable a, by entering 


a:3+4x/i; 


When Maxima displays a complex number, it usually arranges it with 


the term in 7 as the first term, unless the imaginary part is negative 
and the real part is positive. 


(b) Similarly, assign the value 2 — 2i to b. 


(c) Calculate atb, a-b, a*b and a/b. 


The sum and difference are displayed in simplified form, but the 
product and quotient are not. 


To simplify the result of a calculation and display it in Cartesian form 
(or rectangular form, as Maxima calls it) use the command 


rectform. 


(d) Find a*b in Cartesian form, by entering 
rectform(a*b) ; 


Find a/b in Cartesian form similarly. 


(e) Find the real and imaginary parts of a by entering 
realpart (a); 
and 
imagpart (a) ; 
(f) Find the conjugate of a by entering 
conjugate(a) ; 
(g) Find the modulus and the principal argument of b by entering 
abs(b); 
and 
carg(b); 


respectively. 


Notice that the abs command is the same command that you used 
earlier to find the magnitude (also called the modulus, or the absolute 
value) of a real number. 


The name of the carg command is short for ‘complex argument’. The 
argument calculated is the principal argument, that is, the argument 
in the interval (—7, 7]. 


As you saw in Unit 12, any non-zero complex number can be expressed in 
polar form r(cos@+ isin 0), where r is the modulus and @ is one of its 
arguments. The next activity uses this form. 
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Computer activity 65 Working with the polar form of 
complex numbers 


(a) Assign the complex numbers 3 (cos (57/13) + isin (57/13)) and 
4 (cos (77/13) + isin (77/13)) to the variables c and d, respectively. 


(b) Use Maxima to find the moduli and principal arguments of c*d and 
c/d. 


Hint: use the trigreduce command to simplify the initial answers 
given. 


The commands introduced in the activities above are summarised below. 


Working with complex numbers 


Operation Command Example 

Real part of a complex number realpart( ) realpart (1+2*%i) ; 
Imaginary part of a complex number imagpart( ) imagpart (1+2*/i) ; 
Conjugate of a complex number conjugate( ) | conjugate(1+2*%i) ; 
Modulus of a complex number abs( ) abs (1+2*%i) ; 
Principal argument of a complex number carg (M) carg(1+2*%i) ; 
Express a complex number in Cartesian form | rectform(_) rectform(2/%i) ; 


If you began studying this subsection from Activity 27 of Unit 12, return to 
the unit to continue your study. 


13.2 Solving polynomial equations 


In Subsection 3.2 of this Guide you met the solve command that is used to 
find solutions of equations. This command tries to find all the solutions of 
the equation (if they can be found using algebraic manipulation), including 
solutions that are non-real complex numbers. You might remember that in 
Computer activity 20(e) you solved an equation with non-real solutions. 


Computer activity 66 Solving polynomial equations 


(a) Use the solve command to find the solutions of the equation 
+z2+1=0, 


and assign the result to the variable solns. 


Notice that the solutions are given exactly, although in a rather 


unwieldy form! 
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It is often more useful to find decimal approximations to solutions, 


expressed in Cartesian form, which you can do using the float and 


rectform commands. 


(b) Express the solutions found in part (a) in Cartesian form, with each 
real number, and the real and imaginary parts of each complex 
number, correct to two decimal places. 


(c) Find, in Cartesian form, all the solutions of 


24 — 2023 + 17127 — 626z + 962 = 0. 


You may have noticed that in Computer activity 66 all the non-real 
complex solutions occur in complex conjugate pairs. In other words, for 
each of the equations, if a particular non-real complex number z is a 
solution, then so is its complex conjugate Z. You saw in Unit 12 that this 
property holds for any quadratic equation that has real coefficients. In 
fact, it holds for any polynomial equation that has real coefficients. It 
doesn’t hold for polynomial equations that don’t have real coefficients, 
however. For example, the solutions of the polynomial equation 


2 —2=0 


are 1 +i and —1 — i, which are not complex conjugates of each other. 


Plotting solutions in the complex plane 


It is sometimes informative to plot the solutions of an equation in the 
complex plane. You can do this in Maxima using the methods that you 
learned in Subsection 11.1 of this Guide for plotting the discrete points 
that form the graph of a sequence. 


To plot the solutions of an equation, you first need to create a list of the 
points to be plotted, with the first coordinate of each point being the real 
part of a solution, and the second coordinate being the corresponding 
imaginary part. 


The following activity shows how you can do this. 
Computer activity 67 Plotting the solutions of 
Bot =) 
(a) Solve the equation 
z+42+5=0 


and assign the list of solutions to the variable s. 


Each element of the list s is an equation of the form z=value, where 


value is a solution of the equation. 
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(b) Create a second list containing just the solutions of the equation 
(without the z=) and assign it to the variable v, by entering 


v:makelist(rhs(s[k]), k, 1, length(s)); 


The kth element of the first list s is s[k], so rhs(s[k]) is the 
right-hand side of this element, which is the kth solution of the 
original equation. 


You met the makelist command in Subsection 11.1 of this Guide. 
The index variable k starts with the value 1 and ends with the value 
given by length(s), which is the number of elements in list s. This 
ensures that all of the the solutions of the equation are included in the 
list v. 


(c) Create a list of the coordinates of the points to be plotted and assign 
it to the variable pts, by entering 


pts:makelist([realpart(v[k]), imagpart(v[k])], 
k,1,length(s)); 


The kth solution of the original equation is v[k]. The point in the 
complex plane corresponding to this solution has coordinates 


(Re(v [k] ), Im(v[k])), 
which can be expressed in Maxima syntax as 
[realpart (v[k]), imagpart(v[k])]. 


The command above uses makelist to create a list of such points, one 
for each solution of the equation. 


(d) Finally, plot the points found by entering 


wxplot2d([discrete,pts], [style,points], 
[xlabel,""], [ylabel,""]); 


This command plots a set of discrete points. The first two arguments 
of this wxplot2d command are the options that you used to plot 
graphs of sequences in Subsection 11.1 of this Guide. The final two 
arguments specify that the axis labels are to be left blank. 


Notice that since the complex solutions of polynomial equations with 
real coefficients occur in conjugate pairs, the points plotted are 
symmetric about the horizontal axis. 


Computer activity 67 introduced the length command, which is 
summarised below. 
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Length of a list 


Operation Command Example 


Computer activity 68 Finding and plotting the seventh 
roots of unity 


(a) Find the seventh roots of unity by solving the equation 
z =1. 


Give the solutions in Cartesian form, with each real number, and the 
real and imaginary parts of each complex number, correct to two 
decimal places. 


(b) Plot the solutions in the complex plane. 


If you began studying this subsection from Activity 35 of Unit 12, return to 
the unit to continue your study. 


14 Maxima accessibility guide 


This section outlines issues associated with the use of Maxima if you 
interact with a computer using only a keyboard, use a screenreader or 
require different colour settings. This information should help you decide 
whether the wxMaxima interface is suitable for you, or whether you should 
consider using an alternative interface, such as the command-line interface. 
It also explains how to configure Maxima to make it easier for you to use 
and contains alternative versions of Subsection 1.2 and parts of Section 2 
describing how to use the command-line interface. 


Activities within this section are numbered to match corresponding 
activities in Sections 1 and 2. If you decide to use the command-line 
interface, the format of the solutions given at the end of this Guide may 
not exactly match the output you obtain. 


14.1 Choosing and configuring your 
interface 


This subsection describes how wxMaxima can be configured to help those 
with additional requirements. If you find these changes are not sufficient, 
you may wish to use the command-line interface to Maxima instead. 
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Using wxMaxima with a keyboard alone 


The majority of the wxMaxima interface is accessible using a keyboard 
alone. On a Windows computer you can access the wxMaxima menus by 
pressing and releasing the key, navigating through the menus using the 
keyboard arrow keys, and then pressing to make your selection. 
Alternatively, you can press and hold down the key, and then press the 
keyboard key that corresponds to the underlined letter in a menu entry to 
select it. In both cases, you can leave menu navigation without making a 
selection by pressing the key. 


On an Apple Mac you can access the wxMaxima menus by pressing and 
releasing the key while holding down the and keys (or just 
on some keyboards), navigating through the menus using the keyboard 
arrow keys, and then pressing to make your selection. You can 
leave menu navigation without making a selection by pressing the key. 


The wxMaxima toolbar cannot be accessed using the keyboard alone, but 
this is not used extensively in this module and alternative typed 
commands for accessing the functions are given as required. 


Changing wxMaxima colours and fonts 


You can change the appearance of the text used by wxMaxima as follows. 
First, click on the following toolbar icon. 


Alternatively, select Configure from the Edit menu (or, on an Apple Mac, 
select Preferences from the wxMaxima menu). 


In the configuration window that opens, select the following ‘Style’ icon. 


it 


You can then change the colour, style, and sometimes the font and size, of 
individual elements of the wxMaxima interface by first choosing the 
appropriate element from the ‘Styles’ list, and then selecting your 
preferred settings. 


Where a selected element has a greyed-out ‘Choose font’ button this is 
usually because it shares its font with other elements of the interface. For 
example, those elements that form part of the mathematical output of the 
interface share a font that you can change by first clicking on the Math 
font ‘Choose font’ button and then selecting your preferred font and size. 
Most of the remaining elements with a greyed-out font button share a font 
that you can change in a similar way by clicking on the Default font 
‘Choose font’ button. If you choose a large font size, some horizontal 
scrolling may be necessary to read all the text displayed in a worksheet. 
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Unfortunately, the Configuration window is not keyboard accessible, and 
so if you only use a keyboard you may need the assistance of a non-medical 
helper to initially configure the system. 


The Configuration window does not change the colours and fonts used in 
graphs produced by Maxima. Subsection 14.7 explains how to change the 
colours, fonts and line thickness used in graphs. You should read 
Subsection 14.7 in conjunction with Subsection 3.3, which discusses how to 
plot graphs. 


Using a screen magnifier or screenreader 


The text in wxMaxima can become pixelated when you use a screen 
magnifier. An alternative is to increase the size of the text font as 
described above and/or use the zooming facilities on the View menu, which 
you can also access using Alt-I (zoom in) and Alt-0 (zoom out). 


The wxMaxima interface is not accessible to a screenreader. If you use a 
screenreader you might wish to consider using the command-line interface 
to Maxima instead of the wxMaxima interface. You should ensure your 
screenreader is set to read all the punctuation. Details of using Maxima 
with the command-line interface are given below. The graphs produced by 
Maxima are not accessible to a screenreader. You might need assistance 
when using these. 


14.2 Starting, testing and configuring 
Maxima using the command-line 
interface 


This subsection describes how to start, test and configure Maxima when 
using the command-line interface. If you are using the command-line 
interface you should read this subsection instead of Subsection 1.2, which 
this replaces. 


Computer activity 2 Starting Maxima 


Start the command-line interface to Maxima on your computer and keep it 
open to work with as you study this section. 


If you cannot remember how to do this on your computer, see the 
OU Maxima website: learnl.open.ac.uk/site/maxima. 


The command-line interface to Maxima is illustrated in Figure 12. 


91 


Computer algebra guide 


92 


E Command line Maxima 


Maxima 5.30.0 http://maxima.sourceforge.net 

using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING. 
Dedicated to the memory of William Schelter. 

The function bug report() provides bug reporting information. 
(%i1) 


KANS 


Figure 12 The initial command-line interface window 


At the top of the window is some introductory text, followed by 
(hit) 
This is how Maxima labels lines. The i stands for ‘input’: this is input line 


number 1 and is where you can type your first command. These commands 
instruct Maxima to either perform calculations or change system settings. 


The percentage symbol (%) is used by Maxima to indicate objects built 
into the system. There is more about this later. 


The following activity shows you how to use Maxima to perform a simple 
calculation and checks that your system is working properly. Follow the 
steps using Maxima as you read the activity. 


If you encounter unexpected problems when working through this activity, 
check the Frequently asked questions (FAQs) section on the OU Maxima 
website: learn1.open.ac.uk/site/maxima. 


Computer activity 3 Your first Maxima calculation 


Follow these steps to use Maxima to calculate 2 + 3 x $. 


(a) Click anywhere on the Maxima window, to ensure that the text you 
type next is directed to the correct place. 


(b) Type the following exactly as it appears here: 
2+3*4/5; 


including the semicolon (;) at the end. 


The text you type will appear next to the input line number (%i1), as 


indicated by a flashing horizontal line known as the editing cursor. 


(c) Press : 
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If you forget to type the semicolon (;) at the end of the command, 
the editing cursor simply moves to the next line, waiting for you to 
type more. Maxima will continue to do this until you indicate you 

have reached the end of the command by typing a semicolon. So if 
you forgot the semicolon, then type it now and press again. 


Pressing instructs Maxima to calculate the expression. This 
is known as evaluating the expression. 


The result is displayed under the input line, preceded by the label 
(401). The o stands for ‘output’, so this label identifies the line as 
output line 1. It gives the output corresponding to input line 1. 
Throughout each Maxima session the line numbering starts at 1 and 
then increases by one for each new input (and output). 


Maxima then displays another input line number, ready for you to 
enter your next command. 


E Command line Maxima 


Maxima 5.30.0 http://maxima.sourceforge.net 

using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING. 
Dedicated to the memory of William Schelter. 

The function bug report() provides bug reporting information. 
(Sil) 2+3*4/5; ` 


(%01) aS 
5 
(%i2) 


If, in Computer activity 3, Maxima behaved as described, then you have a 
working system to use throughout MST124. 

The result of the calculation in Computer activity 3 was displayed as 
formatted mathematics over several lines, replicating how mathematics is 
usually written. This can, however, cause problems if you are using a 
screenreader. The display behaviour can be changed for every Maxima 
session as described in the following activity. 


Computer activity 4 Configuring the Maxima display of 
mathematics 


You should only do this activity if you wish to change how Maxima 
displays mathematics. 
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(a) First exit Maxima by entering the command 
quit(); 
(on some installations you may also have to close the window). 


(b) Turn off the two-dimensional display of mathematics by adding the 
Maxima command display2d:false; to your maxima-init.mac 
configuration file, as follows. 


(i) Follow the instructions in the Accessibility section of the OU 
Maxima website learn1.open.ac.uk/site/maxima to locate 
your maxima-init.mac configuration file, or create one in the 
correct location as appropriate. 

(ii) Open your maxima-init.mac file in a plain text editor (such as 
Notepad on a Windows computer, or TextEdit on an Apple 
Mac), and add the text 


display2d:false; 


on a new line. There are further details on how to do this in the 
Accessibility section of the OU Maxima website. 


(iii) Save your modified maxima-init.mac, and close your text editor. 


(c) Restart Maxima and enter the calculation from Computer activity 3 
again, that is, enter 


2+3*4/5; 


Notice that Maxima now displays the result on a single line: 


22/5 


Another useful configuration change is to make Maxima display the result 
of a command left-justified against the output line number, instead of the 
default behaviour of displaying it centred on the line. This can be achieved 
by entering the command 


leftjust:true; 


either in the Maxima window, or, if you wish to make the change 
permanent, in your maxima-init.mac file as described above. 


Closing wxMaxima 


When you have finished using Maxima, close it by using one of the 
following methods. 


e Enter quit 0; 


e On a Windows computer, click the usual small cross button at the top 
right-hand corner of the Maxima window. 


14 Maxima accessibility guide 


14.3 Troubleshooting Maxima using the 
command-line interface 


If you are using the command-line interface to Maxima and it does not 
seem to be responding, try the following suggestions. 


Troubleshooting Maxima 


e Check you have ended your previous command with either a 
semicolon (;) or a dollar ($). Maxima will not do anything without 
them! 


e Check that Maxima is not waiting for you to type something. Is a 
question being displayed that you have not yet answered? 


e If Maxima is taking too long to evaluate a command, you can 
interrupt it by pressing Ctrl-C. 


If you do this, you may see the following output, or something 
similar, which you can safely ignore. 


Maxima encountered a Lisp error: 
EXT:GC: User break 
Automatically continuing. 
To enable the Lisp debugger set *debugger-hook* to nil. 


Note that pressing Ctrl-C at other times may close your Maxima 
session. 


If you encounter other problems when using Maxima, check the 
Frequently asked questions (FAQs) section of the OU Maxima website: 
learn1.open.ac.uk/site/maxima. 


14.4 Reusing and editing commands in the 
command-line interface 


In Computer activity 8 you may have entered the expression twice: once to 
evaluate it exactly and once for a decimal answer. This is, however, not 
necessary. When using the command-line interface, you can edit and 
re-evaluate previously entered commands or calculations, as shown in the 
following activity. 


Computer activity 10 Editing a command 


(a) Calculate 2/3 + 5/27 by entering the following line into Maxima. 
2*sqrt (3)+5*sqrt (27) ; 
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(b) Edit your entry in part (a) to calculate 2/3 — 5/27 as follows. 


(i) Press the keyboard up arrow key once, so that the previously 
entered command appears against the current input line number 
(on an Apple Mac you may need to hold down the key while 
pressing the keyboard up arrow key once). 


You can access all the commands you have previously entered in your 


Maxima session using the up and down keyboard arrow keys. 


(ii) Edit the expression by moving the editing cursor using the left 
and right keyboard arrow keys, typing new text and using the 


(backspace) key to make deletions. 
(iii) Press ‘ 


Maxima also permits the result of one calculation to be used in another 
calculation. The symbol % represents the result of the last evaluated 
command. So, for example, if you enter sqrt (%) ; Maxima calculates the 
square root of the last evaluated output. You can also use the output (or 
input) of other lines, by typing the line label in a calculation. For example, 
if you enter %05^3; Maxima calculates the cube of the expression in output 
line 5. 


Computer activity 11 Using the results of previous 
calculations 


(a) Enter 6*2; into Maxima and evaluate it. 


(b) On the following line type 3+%; press . 


Maxima returns 15, the sum of 3 and the previous answer. 


(c) Enter 40! + 5; with | replaced by the output line number 
corresponding to the result of 6*2; entered earlier. 


This command adds 5 to the result of the line referenced. 


14.5 Saving and printing your work using 
the command-line interface 


After working with Maxima you will probably want to save your 
calculations so that you can reuse or read them at a later date. 
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When you are using the command-line interface you cannot save the 
session directly, as you do when saving other files. Instead there are two 
stages in saving your work: 


e If you want to keep a record of the contents of the command-line 
window, then you need to save a transcript of the session. You can then 
print this transcript if needed. 


e If you want to save the values of any variables and functions you have 
defined, so that you can continue working with them in a new Maxima 
session, then you will need to save the current state of the system. Note 
that this does not give you access to the commands used to define the 
variables or functions, so you cannot edit them in the new session. 


These two processes are illustrated in Computer activity 16. 


When saving your work, you need to tell Maxima exactly where to save 
the file containing the session state or your transcript by giving the 
complete pathname of the file, which includes not only the name of the file 
but also the hierarchy of folders in which the file is located. Unless you do 
this, Maxima will try to save the file in the location in which the system is 
installed on your computer. Not only is this not recommended, but as a 
computer user, you may not have sufficient permission to do this. 


The pathname of a file depends on its location, the type of computer you 
are using, and your computer username. In the following activity, we 
assume that a user with the computer username ‘bill’ wishes to save a 
transcript of his Maxima session in a file called transcript.txt within a 
folder called maxima which he created in Computer activity 1. The maxima 
folder is within Bill’s ‘Documents’ area on a recent Windows computer 
(that is, Windows Vista, 7, 8 or 10). In this case, the full pathname of his 
file is 


C:/Users/bill/Documents/maxima/transcript.txt 


(Notice that Maxima separates parts of the pathname with a forward slash 
(/) rather than the usual Windows backward slash (\).) 


If Bill were using an older Windows computer (such as Windows XP), the 
equivalent pathname would be 


C:/Documents and Settings/bill/Documents/maxima/transcript.txt 
Using an Apple Mac the equivalent is 


/Users/bill/Documents/maxima/transcript.txt 


Computer activity 16 Saving your work 


(a) Create a transcript of your current session by following these steps. 


(i) Enter 
writefile("C:/Users/bill/Documents/maxima/transcript.txt") ; 
replacing the pathname with one appropriate to you and your 
system. 
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The path of the file is given in double quote marks. 


The command tells Maxima to save a copy of all the subsequent input 
and outline lines to the file named transcript.txt in your chosen 
folder. 


The output 


#<OUTPUT BUFFERED FILE-STREAM CHARACTER 
C:\Users\bill\Documents\maxima\transcript.txt> 


done 


is given, where the pathname is the one you gave in the command. 


(ii) Enter 
playback () ; 


This command redisplays each of the input and output lines of your 
session, and when used in conjunction with the writefile command 
results in a copy of these lines being stored in your chosen file. 


To playback only a selection of the lines from your session, you can 
use a command of the form 


playback([ first line number , last line number |); 


which plays back only those lines within the given range. 


(iii) Enter 
closefile(); 


This stops the recording of the session into the file. 
The output 


#<CLOSED OUTPUT BUFFERED FILE-STREAM CHARACTER 
C:\Users\bill\Documents\maxima\transcript.txt> 


done 


is given. 


(iv) Find the file transcript.txt on your computer and open it 
using a text editor, for example, Microsoft Word. 


The file contains a copy of the input and output lines from your 
session. (When opened in some text editors the file may appear all on 


one long line.) 
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(b) Save your current session then reload it, by following these steps. 
(i) Enter 
save("C: /Users/bill/Documents/maxima/mysession.lisp", all); 


replacing the pathname with one appropriate to you and your 
system. 


This tells Maxima to save your current session to the file 
mysession.lisp in your chosen folder. (Lisp is the computer 
programming language used to create Maxima, and the save 


command saves your work using commands from this language.) 


The path of the file is output. 


(ii) Close Maxima by entering 
quit(); 
(iii) Restart Maxima in the usual way. 
(iv) Load your previously saved session by entering 
load("C:/Users/bill/Documents/maxima/mysession.lisp") ; 


using the same pathname you used in part (b)(i). 


This tells Maxima to restore the state of the system to that saved in 
the specified file. 


(v) Enter 
values; 


to display the names of all known variables. 


You should see that all the variables from your previous session, if any, 
are known. You might like to check their values are as expected too. 


There is more about variables and the values command in 
Subsection 2.6. 


An alternative way of using the writefile and closefile commands is to 
enter the first at the very start of your session and the second at the end, 
so that everything you do is recorded. 
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Saving your work 


Operation Command 


Stop writing a transcript 


Save the state of a session 
Reload a session 


100 


Start writing a transcript writefile("file path"); 


Playback all the input and output in a session | playback() ; 
Playback a range of input and output playback([ start line number, end line number ]) ; 


closefile(); 


save("file path", all); 
load(" file path"); 


You can also obtain an image of the Maxima window by taking a 
‘screenshot’. To do this on a Microsoft Windows computer, first click on 
the Maxima window, then press Alt-PrintScreen (that is, hold down the 
keyboard key while pressing the PrintScreen key, which may be 
labelled ( PrtScr ) , (Prt Sern) or something similar). This stores the image of 
the Maxima window, which you can then paste into a suitable application. 
If you are using Windows Vista, Windows 7, 8 or 10, you might like to use 
the ‘Snipping Tool’, available by typing ‘snip’ into the Start menu. 


To take a screenshot on an Apple Mac computer, press Cmd-Ctr1-Shift-4 
(that is, hold down the Command) , (Ctrl) and Lt) keyboard keys while 
pressing the key). Next press the space bar, move the camera pointer 
over the Maxima window and click. This stores the image of the Maxima 
window, which you can then paste into a suitable application. You may 
also like to use the Grab application from the Utilities folder. 


14.6 Plotting graphs without wxMaxima 


When using the command-line interface, you should use the plot2d 
command rather than the wxplot2d command described in Section 3.3. 
These commands work similarly, and have the same arguments. The main 
difference is that plot2d opens a new computer window containing the 
graph. 


The graph window can be closed by pressing the usual cross at the top 
right-hand corner, or by pressing Alt-F4. You will not be able to enter 
any commands into the Maxima command-line interface until the graphing 
window is closed. 


The graph can be saved as a .png, .pdf or .svg file using the save button. 
These files can be opened with many image viewing programs and can be 
inserted into word-processed documents. 


14.7 Changing properties of graphs 


You can change the thickness of curves plotted in graphs by including an 
additional argument in the wxplot2d command, of the form 


[style, [lines, 3]] 
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where the number indicates how thick the line should be. By default, lines 
have thickness 1. 


So, to plot a graph of the equation y = 2x with a thicker line, enter, for 
example, 


wxplot2d(2*x, [x, 0, 1], [style, [lines, 311); 


Typing this additional argument each time you want to plot a graph can 
be cumbersome, so you might like to assign this argument to a variable 
that you can reuse. For example, to assign your preferred style to the 
variable thick enter 


thick: [style, [lines, 3]]; 
Then you can plot graphs using commands such as 
wxplot2d(2*x, [x, 0, 1], thick); 


To change the colour of a curve, you can include an argument of the form 
[color, red] as described in Subsection 3.3. Alternatively, you can add a 
further element to the styles argument described above. For example, 
including 


[style, [lines, 3, 5]] 


produces curves with thickness 3 and colour number 5, which is black. The 
colour number codes used by the style argument are given in the table 
below. 


style argument colour numbers 


Number | Colour 


blue 

red 
green 
magenta 


black 


DEN 


cyan 


Note that the value of any variable you define to help set curve properties 
(such as thick above) will be lost when you close Maxima at the end of a 
session unless you save your session. The style argument also only 
changes properties of the curve plotted, not other aspects of the graph 
such as the axes and labels. 


Changing the axes and labels, and ensuring such changes are automatically 
used each time Maxima is started, is more complicated. A computer file is 
provided in the Accessibility section of the OU Maxima website 
learni.open.ac.uk/site/maxima which you can download and configure 
to set the graph properties to your needs. These settings will then be 
automatically used for each graph in every Maxima session. Details of how 
to download and use this are given in the Accessibility section of the OU 
Maxima website. 
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Computer methods for CAS Activities 
in Books A-D 


Computer solution to Unit 5, Activity 17 


First, load the implicit_plot package and set the default plotting options 
to use equal scales on the x- and y-axes. 


(4i1) load(implicit_plot) ; 


(%01) C:\maxima—5.38.1\share\maxima\5.38.1_5\share\contrib\ 
implicit_plot.lisp 
(%i2) set_plot_option(same_xy)$ 


This line is ended with $ to prevent the display of the output of the 
command. 


(a) Assign the equation representing the parabola and the equation 
representing the circle to the variables p and c, respectively. 


(%i3) p:y=2*x^2-10*x+14; 
(403) y= 2x" — 102414 
(hi4) c:x^2-6*x+y^2-8*y+22=0; 
(%04) y?—8y+x?—6x+22=0 


Now plot the two curves. Appropriate axis ranges are chosen to 
ensure the whole of the circle is shown. 


(%25) wximplicit_plot((p,c], [,0,7], [y,0,7]); 


y = 2°x^2-10°x+14 — 
y^2-8*y+x^2-6*x422 =0 ——— | 


(%t5) 3 b 


(%05) 

(b) To find the points of intersection of the parabola and circle, we solve 
the equations representing each simultaneously, for x and y. The 
equations have previously been assigned to the variables p and c, 
respectively. 


(4i6) solve([p,c], [x,y]l); 
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(%06) [[x = 1.304970075461879, y = 4.356193142057383], 
[x = 1.664051355206847, y = 2.897620365246265], 
[x = 3.121372031662269, y = 2.272207016670814], 
[x = 3.909606299212598, y = 5.473980309423348]] 


So, to 2 d.p., the four points of intersection are 
(1.30, 4.36), (1.66, 2.90), (3.12, 2.27) and (3.91, 5.47). 


The maximum number of intersections of a circle and a parabola is 
four. This fact can be justified by the geometric argument that each 
‘half’ of the parabola (each part on one side of the vertex) can meet 
the circle at most twice. 


Alternatively, if the equation of the parabola is used to substitute for 
y in the equation of the circle, then an equation for x containing 
terms in zf, z’, x?, x and a constant term is obtained, and such an 
equation can have at most four solutions (by a result that you will 


meet in Unit 12). 


Computer solution to Unit 7, Activity 23 


(a) 


(b) 


The function f is defined as follows. 
(Hit) £(x) :=(sqrt (1+x%2) )/8+(2-x)/16; 


IFE ee 


h f = ; 
(%01) f(x) 3 16 
The graph of f(x) is plotted below. 
(%i2) wxplot2d(f(x), [x,0,21); 


0.28 
© 0.275 
= 0.27 
Á (0.265 
+ 0.26 
Z 0.255 
Cht2) & 0.25 
% 0.245 
= 0.24 
a 0.235 
0.23 
0 0.5 1 1,5 2 
X 
(h02) 


The value of x that gives the minimum value of f(x) seems to be 
about 0.6. 


The derivative of f can be found and assigned to df (x) as follows. 
(413) df (x) s="" (diff (f (x);x)]; 


x Al 
%o3) af(x) := ——— — — 
ee) SUG) aaa i 
x 1 
So the derivative of f is f'(x) = == — —. 
fis f(z) 8Vx2+1 16 
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(e) 


(£) 


First, we try to use the solve command to find the solution of 
df (x)=0. 


(i4) solveldf(x)=0); 


z 
hoo) [x= YE 4 

2 
Maxima returns a rearranged version of the equation, but does not 
solve it. The solve command cannot find an exact solution of 
equation df (x)=0, so we look for an approximation to the solution 


instead. 

To do this we need to find an interval containing the solution. We 
know from the graph of f(a) that its gradient is negative to the left of 
the minimum, and in particular at x = 0, so df(0)< 0. Also, the 
gradient of f(x) is positive to the right of the minimum, in particular 
at x = 1, so df (1)> 0. So a solution of df(x)=0 must lie in the 
interval [0,1]. An alternative way to find this interval is to plot a 
graph of df (x). 

(415) find-root (df (x) ,x,0,1); 

(405) 0.5773502691896258 


The solution of the equation f'(x) = 0 is x = 0.577 (to 3 s.f.). 


The man should join the road 0.58km (to the nearest 10 metres) 
along from the point that is closest to his initial position. 


(Because the graph of f indicates that the minimum value of f(a) occurs 
at a local minimum rather than at one of the endpoints of the interval 
[0,2], in the solution above we don’t consider the values of f(0) and f(2).) 


Computer solution to Unit 8, Activity 48 


(a) 


The domain of f is the interval [—1, co). 


(b) The graph of f can be plotted as follows. 


(h11) £ (x) :=(x%2-10*x+15)/(1+sqrt (xt+1)); 
x?—10x+15° 


t+yxti ’ 


(%12) weplotod(f(s), [z,;=1;15]); 


(hoi) f{x):= 


am 30 T T T T T T T T 
+ 
= 25 7 
+ 
= 20 4 
8 15 4 
(442) © 49 | 
+ 
es 5 J 
oO 
a 0 i 
X -5 i L L L L fi fi L 
0 2 4 6 8 10 12 14 
X 
(%02) 


(c) 
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The z-intercepts of f can be found by solving f(x) = 0. 

First, solve the equation, assigning the result to the variable solns. 
(413) solns:solve(f(x)=0) ; 

(403) [x =5—V10,x = v10 + 5] 


Next, extract the two solutions from the list and assign them to the 
variables sol1 and sol2, respectively. 


(414) soli:rhs(solns[1]); 

(404) 5—V10 

(%15) sol2:rhs(solns[2]); 

(405) V10+5 

Finally, express the solutions as decimal approximations. 
CAi6): float(sol1); 

(%06) 1.83772233983162 

(417) float(s012); 

(%07) 8.16227766016838 

So the two z-intercepts of f are 1.84 and 8.16 (to 3 s.f.). 


The area between the graph of f and the x-axis, between the two 
x-intercepts, is 


a z2 — 10x +15 
= = GE: 
1.837... l+vz+1 


We shall concentrate on evaluating 


o r2 — 10x +15 
eee T, 
1:837.. l+yr+i 


that is, 
8.162... 
| Fod, 
1.837... 
remembering the required answer is the negative of this integral. 


Maxima can integrate this function exactly. 


(%18) integrate (f(x) ,x,sol1,so12) ; 


(%08) — (450 log (lv 6 + 1) — 450 log (ve-v + 1) 


— 15102 + /6 — V10 (s8vi0 +66) 
+ \/V10+6 (8v0 - 66) ) J15 


(%i9) float(%); 
(%09) —12.37959499949233 
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So 


8.162... 
f f(x) dz = —12.4 (to 3 sf). 
1.837... 


Remembering the required area is the negative of this, the required 
area is 12.4 (to 3 s.f.). 


Computer solution to Unit 8, Activity 49 


The area of the cross-section, in square centimetres, is 
1.5 
f $V x? +1dz. 
-1 


The integrate command cannot evaluate this definite integral, so we find 
an approximate value instead. First, assign the integrand to the variable y. 


(hal) yi (1/3) eaqrt (x*St1); 
xs 1 
3 
Then find an approximate value for the definite integral. 


(hoi) 


(%1i2) quad_qags(y,x,-1,1.5); 
(%02) [0.939242718890935, 7.102141097448111 10 *°, 231, 0] 


The value of the definite integral is 0.939 (to 3 s.f.). So the volume of 
plastic required to make one metre of the edging is approximately 93.9 cm. 


Computer solution to Unit 8, Activity 50 
The change in the displacement of the object from time t = 0 to t = 10 is 


10 2 
f sin | — | dt. 
0 150 


The integrate command calculates this definite integral in terms of 
values of the Error function, erf. Instead, we shall find an approximation 
to its value using the quad_qags command. 


(%il) quad_qags(sin(t%2/150) ,t,0,10) ; 
(%01) [2.15266541148388, 2.389938704143591 10 **, 21, 0] 
So the value of the definite integral is 2.15 (to 2 d.p.). 


The displacement of the object at time t = 10 is (6 + 2.15) m = 8.15 m, to 
the nearest centimetre. 
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Computer solution to Unit 9, Activity 28 


(a) The matrix form of the equations is 


22 3\ /z 2 
23 Al {yl=l| 1 
4 7 10/ \z 3 


Assign the coefficient matrix to the variable A, and the vector formed 
by the right-hand sides of the equations to b. 


(hill) A:matrix([2,2,3],[2,3,4],[4,7,10]); 


22 3 
(%01) |2 3 4 
4 7 10 


(%12) b:matrix([2],[-1],([3]); 


2 
(%02) |-1 
3 
The inverse of A can be found as follows. 


(113) invert(A); 


1 1 _1 
2 2 
(%03) |-2 4 —1 
1 —3 i 
So we have 
=l 1 1 
22 3 1 2-2 
2.3 A = | -2 4 —1 
4 7 10 1 —3 1 
Therefore 
1 1 
yļ=ļ-2 4 ak 
z 1 —3 1 3 


This can be calculated using 
(4i4) invert(A).b; 


0 
(404) |—11 
8 


The solution is z = 0, y= —11, z = 8. 


Note that an alternative way to solve these equations using Maxima is 
to use the solve command, as follows. 


(4i5) solve([2*x+2*y+3*z=2, 2*x+3*y+4*z=-1, 4*x+7*y+10*z=3], [x,y,z]); 
(%05) [[x = 0,y = —11,z = 8]] 
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(b) The matrix form of the equations is 


3 2 1\ fa 1 
43 1){y|=[2 
7 5 1/ \z 1 


Assign the coefficient matrix to the variable A, and the vector formed 
by the right-hand sides of the equations to b. 


(4i6) A:matrix([3,2,1],[4,3,1],[7,5,1]); 


3 2 1 
(%06) |4 3 1 
T5 f 
(%4i7) b:matrix([1],[2],[1]); 
1 


(%07) |2 
1 


The inverse of A can be found as follows. 
(4i8) invert(A); 
2 —3 1 
(%08) —3 4 —i1 
1 1 —1 
So we have 


3 2 1 2 —3 1 
4 3 1 =|-3 4 =l 
7 5 1 1 1 -l 


Therefore 


This can be calculated using 
(%i9) invert (A) .b; 
=3 


(%09) 4 
2 


The solution is x = —3, y = 4, z = 2. 


Note that an alternative way to solve these equations using Maxima is 


to use the solve command, as follows. 


(%i10) solve([3*x+2*y+z=1, 4*x+3*y+z=2, 7*x+5*y+z=1], [x,y,z]) 


(4010) [[x —3,y 4,2 2]] 
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Computer solution to Unit 9, Activity 29 


(a) The matrix form of the equations is 


i, 2 3 x 5 
23 1] lyl= {6 
4 7 7 z 1 
1 2 3 
The coefficient matrix is | 2 3 1 
4 7 7 


Its determinant can be found using Maxima as follows. 
(Lil) determinant (matrix([1,2,31,[2,3,11,[4,7,71)); 
(fol) 0 


The determinant of the coefficient matrix is zero, so the coefficient 
matrix is non-invertible and the equations do not have a unique 
solution. 


(b) The matrix form of the equations is 


2 2 3 x 5 
23 1] ly|=({6 
4 7 7 z 1 
2 2 3 
The coefficient matrix is | 2 3 1 
4 7 7 


Its determinant can be found using Maxima as follows. 
(%12) determinant(matrix([2,2,31,[2,3,11,[4,7,71)); 
(%02) 14 


The determinant of the coefficient matrix is 14. This is non-zero, so 
the coefficient matrix is invertible and the equations have a unique 
solution. 


(c) The matrix form of the equations is 


1 1 -l 1 £ 0 
1-1 i illiyi |2 
2 3 —l wl 11 
6 2 4 0 z 0 
1 1 —1 1 
: as 1 —1 1 1 
The coefficient matrix is 9 0 3 1 


6 2 4 O 
Its determinant can be found using Maxima as follows. 


(k13) determinant (matrix [1,45=1.4), [t.=1.4,.1) 4.123053 5=1], [6,254.0] >: 
(403) 0 
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The determinant of the coefficient matrix is zero, so the coefficient 
matrix is non-invertible and the equations do not have a unique 
solution. 


Computer solution to Unit 10, Activity 20 


(a) 


Tn = —4(1.2)} (n= 1,2,3,...). 


The first 30 points on the graph of the sequence (£n) can be plotted 
as follows. 


First, define the sequence. 
(4i1) x[n] :=-(1/3)*(1.2)%n; 


Cooly Bi (-3) 12> 


Then generate the points to be plotted and plot them. 
(⁄i2) xpts:makelist([n,x[n]],n,1,30)$ 


(%13) wxplot2d([discrete,xpts], [style,points], 
[xlabel,"n"], [ylabel,"x[n]"]); 


0 CeeS sooo, T T 
-10 + bd 4 
-20 L °. 4 
-30 } e 7 
-40 L e | 
(AES) -50 L e 
-60 L 4 
-70 L 4 
-80 j ji 1 f 1 ® 


x[n] 


(h03) 
The sequence (zn) is decreasing and £n — —oo as n — oo. 
i,=5(=-09)" Al) 


The first 30 points on the graph of the sequence (yn) can be plotted 
as follows. 


First, define the sequence. 

(414) yn] :=5*(-0.9)4n; 

(404) yo := 5(—0.9)" 

Then generate the points to be plotted and plot them. 
(%i5) ypts:makelist([n,y[n]] ,n,1,30)$ 


(%16) wxplot2d([discrete,ypts], [style,points], 
[xlabel,"n"], [ylabel,"y[n]"]); 
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O-]-NWAH 
T 
e 
e 
L 


yin] 


(%t6) dp get | 


(%06) 

The sequence (yn) alternates in sign and yn > 0 as n > oo. 

Zn = —7p(-2.5)” (n =1,2,3,...). 

The first 15 points on the graph of the sequence (zn) can be plotted 
as follows. 

First, define the sequence. 


(4i7) zn] :=-(1/10)*(-2.5)4n; 


Oot): 2. = (-3) (—2.5)” 


10 
Then generate the points to be plotted and plot them. 
(%18) zpts:makelist([n,z[n]],n,1,15)$ 
(4i9) wxplot2d([discrete,zpts], [style,points], 
[xlabel,"n"], [ylabel,"z[n]"]); 
100000 T T T T T T T 
80000 p 4 
60000 p J 
40000 - 7 


(%+9) 20000 + e 4 


Opeeeeeeeeee®, 4 


-20000 + J 
-40000 fi fi I L L L È 


z[n] 


(%o9) 


The sequence (zn) alternates in sign and is unbounded. 
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Computer solution to Unit 10, Activity 21 
(a) an =17-— U (n=1,2,3,...). 


The first 30 points on the graph of the sequence (an) can be plotted 
as follows. 


First, define the sequence. 


(%11) aln]:=17=(1/3)*(1:2)^n; 
1 
(hot) an := 17 — 51.2" 
Then generate the points to be plotted and plot them. 


(%i2) apts:makelist([n,a[n]],n,1,30)$ 


(4i3) wxplot2d([discrete,apts], [style,points], 
[xlabel,"n"], [ylabel,"a[n]"]); 


vevedeovedoooes, 


(463) 


(03) 
The sequence (an) is decreasing and an — —oo as n — oo. 
(b) bn =5(-0.9)" +45 (n=1,2,3,...). 


The first 30 points on the graph of the sequence (bn) can be plotted 
as follows. 


First, define the sequence. 

(hi4) bln] :=5*(-0.9)%n+45; 

(%04) by := 5(—0.9)" + 45 

Then generate the points to be plotted and plot them. 
(4i5) bpts:makelist({[n,b[n]],n,1,30)$ 


(4i6) wxplot2d([discrete,bpts], [style,points], 
[xlabel,"n"], [ylabel,"b[n]"]); 
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(%t6) 44 4 Pee ee 


b[n] 
& 


0 5 10 15 20 25 30 


(06) 
The terms of (bn) alternate either side of 45 and bn — 45 as n > oo. 
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Solutions to Computer activities 


Solution to Computer activity 3 
(hil) 2+3*4/5; 


, 22 
(%o1) 5 
hare. 
5 5 


Solution to Computer activity 5 
(a) AEL) 2771/3); 

(hol) 3 
(b). (eae) 27^1/3; 

(%02) 9 


Without any brackets, 27^1/3 is interpreted according to the 
BIDMAS rules, as (271)/3 which is equal to 9. 


Do not be surprised if your line numbering differs from ours. Each of our 
solutions is prepared as a separate Maxima session with the line numbering 
starting at 1 (unless it is a continuation of a previous activity); your 
numbering will continue to increase for the duration of your session. 


Solution to Computer activity 7 
(a) (%i1) 1309/3451; 
u 
29 
(b) (%i2) 1309.0/3451; 

(402) 0.3793103448275862 

So 1309.0/3451 ~ 0.379 310 344 827 5862. 
(c) (%13) float(1309/3451) ; 

(403) 0.3793103448275862 

So 1309/3451 ~ 0.379 310 344 827 5862. 


Your answers may differ in the last decimal place or so to those given here. 
This is due to small differences in the way floating-point numbers are 
computed by different systems. 


(hot) 


Solution to Computer activity 8 
(a) (%i1) sqrt (147); 

(%o1) 7V3 

(%i2) float (sqrt (147)) ; 

(%o2) 12.12435565298214 

So /147 = 7/3 ~ 12.124355 652 982 14. 
(b) (413) (8/5)*(2/3); 

(%03) — 

(494) float ((8/5)"C2/3)); 

(%04) 1.367980757341358 


ace 4 
So (2) = — x 1.367 980 757 341 358. 
5 52/3 


(c) (%i5) 2^500; 


(%05) 32733906078961418700 1318969682[9 1digit s] 
217256545885393053328527589376 
(%16) float(2^500}; 


(%06) 3.273390607896142 101°° 
So 2°00 ~ 3.273 390 607 896 142 x 101%. 


Solution to Computer activity 9 
(a) (%i1) %pit%pi/3+he%2; 
(hot) = + he? 
(%12) float (%pi+4pi/3+he!2) ; 
(402) 11.57784630371704 


4 
Sori 5 Les = 4+ e? ~ 11.577 846 303 717 04. 


Solution to Computer activity 10 
(a) (%ii) 2*sqrt (3)+5*sqrt (27) ; 

(401) 17V3 

So 2V3 + 5V27 = 17v3. 


(b) (%i2) 2*sqrt(3)-5*sqrt (27) ; 
(%o2) —13V3 
So 2/3 — 5/27 = —13V3. 


Solutions to Computer activities 
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Solution to Computer activity 11 


At the end of this activity, your worksheet should look like the following, 
although line numbers may differ. 


(hi4) 6*2; 
(%04) 12 
(%15) 3+%; 
(%05) 15 
(hi6) ho4+5; 
(406) 17 


Solution to Computer activity 12 
(a) (%i1) a:sqrt(8) ; 
Go 2 
(b) KAD a; 
(ho2) 25 
(c) (%i3) braxsart (2); 
(%03) 4 
(d) hi4) a:sqrt(7); 
(404) VT 
(e) (hid) 
(05) 
(£) (%16) values; 
(406) [a,b] 
(œ) (hi7) killa); 
(%07) done 
(h) i8) a; 
(%08) a 
G) AALO b; 
(%09) 4 
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Solution to Computer activity 13 
(a) (%i1) a:42; 
(%o1) 42 
(b) (%i2) b:a*3-3*a+4; 
(%02) 73966 
(c) (%i3) c:float(sqrt(b)) ; 
(403) 271.9669097519035 


Solution to Computer activity 14 
(a) (%i1) fpprintprec:8; 
(hoi) 8 


(This sets the number of significant figures to be shown in the display 
of decimal numbers.) 


(b) (%i2) float (%pi) ; 
(%02) 3.1415927 
So the value of m is 3.141 5927 (to 8 s.f.). 


Solution to Computer activity 15 


The circumference of a circle of radius 5 is 31.42 (to 2 d.p.) and the area is 
78.54 (to 2 d.p.). 


Solution to Computer activity 17 
(a) (%i1) ? float; 


-- Function: float (<expr>) 
Converts integers, rational numbers and bigfloats in 
<expr> to floating point numbers. It is also an 
‘evflag’, ‘float’ causes non-integral rational numbers 
and bigfloat numbers to be converted to floating point. 
There are also some inexact matches for ‘float’. 
Try ‘?? float’ to see them. 
(%01) true 
(b) (i2) ?? sqr; 
O: isqrt (Functions and Variables for Number Theory) 
1: sqrt (Root Exponential and Logarithmic Functions) 
2: sqrtdenest (Package sqdnst) 


3: sqrtdispflag (Functions and Variables for Display) 
zee: 


-- Function: sqrt (<x>) 
The square root of <x>. It is represented internally by 
‘<x>4(1/2)’?. See also ‘rootscontract’ and ’radexpand’. 
(%02) true 
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Solution to Computer activity 18 
(a) (%i1) expand( (2+x)*(3*#x+4)%3 ); 
(%o1) 27x* + 162x? + 360x? + 352x + 128 
So (2+ x)(3x + 4)’ = 27x4 + 162z? + 360x? + 352x + 128. 
(b) (%i2) factor( 2*x^3-7*x^2-10*x+24 ); 
(%02) (x — 4)(x + 2)(2x — 3) 
So 223 — 7x? — 10x + 24 = (x — 4)(x + 2)(2x — 3). 
kel KAAS pele sine: 
dl 
x+1 
(d) (%i4) fullratsimp(p) ; 
(ied) x?—x+1 
(e) (%i5) subst (3,x,p); 
(405) T 


(%03) 


Solution to Computer activity 19 
(a) (i1) r:x^2+2*x=4; 
(Yoi) x? +2x=4 
(b) (412) Iħs(r); 
(%02) x? +2x 
(c) (419) rhs(r); 
(%03) 4 


Solution to Computer activity 20 
(a) (i1) solve( x^2+2*x=1 ); 
(hoi) [x =-v2-1,x=v2-1] 
The solutions of the equation are z = — v2 — 1 and z = v2 — 1. 
(b) (%i2) solvel 4*x^2-12*x+9=0 ); 
(%02) [x = A 
3 
5- 
(c) (%i3) solvel 2*x^3+x^2-5*x+2=0 ); 


The solution of the equation is 7 = 


1 
(%03) [x me 1] 


The solutions of the equation are z = 5, xr = —2 and x = 1. 
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(d) (%i4) solve( x^9+2*x-4=0 ); 

(ies) [O0 = x? Aad 

Maxima was unable to solve this equation. 
(e) (%i5) solve( x^2-4*x+5=0 ); 

(%05) [x=2—%i,x = hit 2] 


The two complex solutions are x = 2 — i and x = 2 + i. There are no 
real solutions. 


Solution to Computer activity 21 
(a) TAL) At (2,855,711); 
(loi) |2,3,5, 7,14] 
(b) (%i2) A[4]; 
(%02) 7 
So the fourth element of A is 4. 
(c) (i) i3) solns:solve(x^2-x-1=0) ; 


5—1 5+1 
(%03) ea ee 
2 2 
— 1 il 
So the solutions are x = a and x = = 8 


(ii) (%i4) float(solns); 
(%04) [x = —0.61803398874989, x = 1.618033988749895| 
So decimal approximations to the solutions are 
x = —0.618 033 988 749 89 and x = 1.618 033 988 749 895. 
(iii) (415) rhs(solns[1]) ; 
V5-1 


2 


Ck05)- = 


Solution to Computer activity 22 
(a) (%il) eq:3*y+4=3*x^2+9*x; 
(%o1) 3y+4= 3x? + 9x; 
(b) (%i2) eq-4; 
(%02) 3y = 3x? + 9x — 4; 
(c) (%13) 7/3; 


3x? +9x-—4 


(h03) 3 
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Solution to Computer activity 23 
(a) (hit) eqn:c=(2*at+b+3*c)/(a-b) ; 

3c+b+ 2a 
Pm 
(b) (%i2) solve(eqn, a); 
(b+3)c+b 

C— 2 
(b+3)c+b 

GZ 


(%01) 


(%02) [a= ] 


So a = 


Solution to Computer activity 24 
(a) (i1) eq1:3*x^2+2*y=20; 

(foi) 2y + 3x? = 20 

(hi2) eq2:4*x-3*y=-4; 

(%02) 4x -— 3y = —4 
(b) (%i3) solvel [eqi, eq2], [x,y] ); 


26 68 
% 2, 4j, a Sn 
(403) [[x = 2,y = 4], [x Rl > all 
The two solutions are x = 2, y=4 and z = -%, y= — $3. 


Solution to Computer activity 25 
(a) (%i1) wxplot2d(4*x*2-8*x+2, [x,-1,2]); 


14 
12 + 4 
a 10 F 4 
x 8p - 
$ 6p 4 
ety = al l 
st 2L 
0 L A 
-2 ji Í L 1 1 
-1 0.5 0 0.5 1 1.5 2 
xX 
(%01) 


(b) (412) vzplot2d(4*x^2=-8*x+2; [x,=1,2], [y.-5,5]); 


plot2d: some values were clipped. 
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N 
+ 
x< 
čo 
A 
t2) & 
St 

4E 4 

L Í L L L 
-1 -0.5 0 0.5 1 1.5 2 
X 
(%02) 


(e) A23) wxplot2d(4*x*2-Sex+2, [x,-1,2], Ly,-5,5), 
[color, green]); 


plot2d: some values were clipped. 


4 | ` J 
y 2} N J 
z Ng Z 
a op E a 
t3) Ñ Rd 
& -2+ a 4 
-4 L 4 
| 1 1 L | 
-1 -0.5 0 0.5 1 1.5 2 
X 
(%03) 


Solution to Computer activity 20 
(a) (hi1) wxplot2d([15*t^3+10*t, 60*t], [t,0,2]); 


140 T T T 
15*t43+10*t — 
60*t 


(%t1) 60 


(%01) 
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(b) (%i2) wxplot2d([15*t%3+10*t, 60*t], [t,0,2], 
[color, green, magenta] ) ; 


140 
120 
100 


(%t2) 60 


(ho2) 


l 45*43410% —— 
60*t 


0.5 


(c) CAi3) wxplot2d( [15*t*3+10*t, 60*t], [t,0,2], 
[color, green, magenta], 
[legend, "si", 


140 
120 
100 

80 


(%t3) 60 


(403) 


isg] yJ; 


(d) (hi4) wxplot2d([15*t^3+10*t, 60*t], [t,0,2], 
[color, green, magenta], 
[legend, "si", "s2"], 
[xlabel, "time, t (h)"], 
[ylabel, "displacement, s (km)"]); 


140 


(%t4) 


displacement, s (km) 


(%04) 


122 


120 
100 


80 
60 
40 
20 

0 


s1 a 
i s2 


0.5 


Solution to Computer activity 27 


(a) 


(b) 


(f) 


(Ail) f(x) :=x%5-34x%2+4; 
(401) £(x):=x°- 3x? +4; 
(12) £12); 
(%02) 24 
Sof (2) = 24. 
(413) g(x):= (x-1)/(x-2); 
x—1 
(%03) E(X) := nt 
(hi4) g(3); 
(h04) 2 
So g(3) =2. 
(hi5) g(2); 
expt: undefined: 0 to a negative exponent. 
#0: g(x=2) 
-- an error. To debug this try: debugmode(true) ; 
(416) MEI; 
(%06) h(1); 
(LiT) functions; 


(%07) [£(x), g(x)] 


Solution to Computer activity 28 


(a) 


(%i8) wxplot2d(f(x), [x,-2,2]); 


x^5-3*x^2+4 
S o 
T T 


(%t8) 
-20 + 
-30 + 
-40 L L L i 
2 15 -1 -05 0 
X 
(%08) 
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(b) (%i9) wxplot2d(g(x), [x,-5,5]); 


800 T T T T T 
600 J 
400 4 
200 4 
0 ma =| 
-200 J 
-400 J 
-600 J 
-800 J 

-1000 L L i L fi 


(x-1)/(x-2) 


(%t9) 


(%09) 
(c) KALLO) wxpletad eG). [x;=5;5]; [y>=5;5]1); 


plot2d: some values were clipped. 


4 I= 
OP 
N 
x 
= OF 
(4t10) x 
-2 8 
4 L 
L if | L i 
4 2 0 2 4 
X 
(%010) 


Solution to Computer activity 29 
(a) (411) Hes; 
(hot) Ke” 
(b) (%i2) float (%e%5) ; 
(%o2) 148.4131591025766 
So e° = 148.41 (to 5 sf). 
(c) (A13) exp(5); 
(403) He” 
(d) (424) float( log(8) ); 
(o4) 2.079441541679836 
So log(8) = 2.079 (to 4 s.f.). 
(e) (%i5) solve( 3^(2*x)=5 ) 


log(5) 
2log(3) 


] 


The only real solution is 


Solution to Computer activity 30 


(a) 


(b) 


(411) url log (x*3+24#x“2+x)-log(x) )/log(x+1); 
log(x? + 2x? + x) — log(x) 


(hot) log(x + 1) 

(412) radcan(u); 

(%02) 2 

Notice, 

ln(z? + 2x? + x) — In(z) 
ln(x + 1) ne KET tla) 
In r? z? +2? +r +2 
7 (x De In(a +1) 
Te 2 need) (n(x? +2 +1) )) 

_ 2in(x +1) 
~ In(x +1) 


=2. 
(%13) v:log(a)+2*log(b)-log(c) ; 
(%03) —log(c) + 2log(b) + log(a) 


(%i4) logcontract(v) ; 


ab? 
(404) log (=) 


Solution to Computer activity 31 


(a) 


(%i1) cos(%pi/6) ; 


3 
(%01) v3 


T 
So cos (=) = TA 
(%12) sin(45*%pi/180) ; 


1 
(%02) Va 


1 
So sin 45° = 


v2 
(%i3) sec(3.4); 
(%03) —1.034342024711439 


So sec(3.4) ~ —1.034 342 024 711 439. 
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(d) (%i4) asin(sqrt(3)/2) ; 
T 
(h04) 3 


So sin”! (5) = 7 


(e) (4i5) esc{hpi/7); 
T 
tos) ese (Z) 

(f) (426) tan(-%pi/12); 


T 
(406) —tan (5) 
(€) (%i7) acos(2.0); 
(%07) 1.316957896924817%i 
So cos 12.0 ~ 1.316 957 896 924 817i. 
(h) (4i8) acos(2); 


(%08) acos(2); 


Solution to Computer activity 32 
(a) (%i1) solve(cos(x)-0.6=0) ; 


rat: replaced -0.6 by -3/5 = -0.6 
solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


3 
(Joi) [x = acos G) 


3 
So a solution is x = cos™ t! G) 


(b) (i2) float) 
(%02) [x = 0.9272952180016121] 
So a solution is x = 0.927 (to 3 s.f.). 
(c) (i3) solve(sin(3*x+4)=1/2) ; 


solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


wa ks 
== 
a 18 
a —24 
S tones , 
O a Solution 1S v% 18 


(d) (%i4) solve(cos(2*x+4)-sin(x/3)=0) ; 
. x a 
(%04) [sin (=) = cos(2x + 4)| 


Maxima cannot find an exact solution of this equation, so returns the 
equation slightly rearranged. 
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(a) 


(i) hi1) f(x) :=cos(2*x+4)-sin(x/3) ; 
te 

(%o1) f(x) := cos(2x + 4) — sin (=) 
(ii) €4i2) weplot2d(f(x), [x,-2,2]); 


2 T T T T T T T 
oe 15 F J 
X 
= 
o 
F 
(t2) x 
N 
D 
=] 
(S 
x 
(%02) 


Since f (0)< 0 and f (1) > 0, there must be a solution in the 
interval (0, 1). 

(ii) (fae) find root(i(a), x, 0; 1); 
(%03) 0.4274333882308137 


So a solution of cos(2x + 4) — sin(x/3) within (0,1) is 0.427 
(to 3 s.f.). 


(4i4) g(x):= (sin(x+2))%2-3*x; 
(%04) g(x) := sin(x + 2)? — 3x 
(415) wxplotad(e(x), [x,=2,2)); 


6 T T T T T T T 
4 4 
x 
ge Br - 
` 
ai | j 
WES) E al | 
7?) 
4 + 4 
6 1 1 1 1 1 1 
2 A5 -1 05 0 05 1 #15 2 
, x 
(%05) 


Since g(0)> 0 and g(1)< 0, there must be a solution in the 

interval (0, 1). 

(If you chose an interval other than (—2,2) then your graph will look 
different from the one given above. If, in particular, your graph does 
not cross the x-axis, then you will need to plot the graph of g again 
using a larger interval.) 


(4i6) find_root(g(x), x, 0, 1); 
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(406) 0.2135670173061208 
So a solution of sin?(x + 2) = 3x is x = 0.214 (to 3 s.f.). 


Solution to Computer activity 34 
(a) (i1) trigreduce( sin(x)’2 ); 
(ot) 1 — cos(2x) 
2 
(b) (%i2) trigexpand(%) ; 
sin(x)? — cos(x)? +1 
2 
(c) (%i3) trigsimp(%) ; 
(%03) sin(x)? 


(402) 


Solution to Computer activity 35 
(f4i1) trigrat( sin(4*x)/sin(x) ); 
(%01) 2cos(3x) + 2cos(x) 

sin(x) 


So = 2 cos(3x) + 2 cos(x). 


Solution to Computer activity 36 
(a) (411) f:cos(5+x); 
(%01) cos(5x) 
(hi2) g:trigexpand(f) ; 
(%02) 5cos(x)sin(x)* — 10cos(x)*sin(x)? + cos(x)® 
So cos(5z) = 5 cos(x) sin*(x) — 10 cos? (x) sin? (x) + cos°(x). 
(b) Choosing the range —27 < x < 2r to plot the graphs, we obtain 
(%13) wxplot2d([f,g],[x, -2*/%pi, 2*%pil); 


j 
T T T 
EEA E ei 

0.5 
OF J 
(4t3) 
-0.5 F J 
| | | | | | 
6 4 2 0 2 4 6 
(%03) ü 


The two graphs appear identical, which suggests the two expressions 


are equal. 


Solution to Computer activity 37 
(a) (%i1) load(implicit_plot) ; 
(%01) C:\maxima—5.38.1\share\maxima\5.38.1_5\share\contrib\ 
implicit_plot.lisp 
(The output from this command is the location of the package on 
your computer system, and may differ from that shown here.) 


(b) (X12) wximplacit plot(.*2+y"2=1, [x,-2,;2], Ly,-2,21); 


2 T T T T T T T 


155 4 


(%t2) 


(%02) 


Solution to Computer activity 38 


(It is assumed the implicit-plot package has already been loaded.) 


Solutions to Computer activities 


(hill) wximplicit_plot(x%2+y*%2=1, [x,-2,2], [y‚-2,21, same- zy); 


(%t1) 


(%o1) 


Solution to Computer activity 39 

(It is assumed the implicit plot package has already been loaded.) 

(a) (%i1) set_plot_option(same_xy)$ 

(b): (%12) . wximplicit plettla1) dy 23, [x;=1,3],. [y;-2,2)); 
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(he) 


-2 1 fi 1 fi 1 i fi 
105 0 05 1 15 2 25 3 


X 


(%02) 


Solution to Computer activity 40 


(It is assumed the implicit-plot package has already been loaded, and 
the plot settings changed to use equal scales for the z- and y-axes.) 


(a) (i1) line: y=3*x-5; 
(%01) y=3x—5 
(%i2) circle: (x-2)%2+(y-3)*2=2; 
Geo) p=) dla =2 
(b) (%i3) wximplicitplot([line, circle], [x,0,5], [y,0,5]); 


i y = 3/5 a 
Z x-2)4 — p 
4E J 
3 E = 
> 
CHESS) 2L J 


(%03) 


Solution to Computer activity 41 


(It is assumed the implicit_plot package has previously been loaded and 
the plotting options set for equal-scale axes.) 


(a) (4i1) wximplicit_plot (x*2+4*y*2=1, [x,-2,2], [y,-2,2]); 
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(c) 


0.5 - -= 
> OF as, = 
(%t1) a 


2 15 -1 05 0 05 1 1.5 2 
x 


(ho1) 

(The curve is an ellipse.) 

(4i2) wximplicit_plot (x^2-4*y^2=1, [x,-4,4], ly,-4,4]); 
4 T E. 


3 F sH 


Qk 4 


wD | Pi i | 


(h02) 
(The curve is a hyperbola.) 
(4i3) wximplicit_plot (y^2=x^3-2*x+1, [x,-4,4], [y,-4,4]); 


Solutions to Computer activities 
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a 
ces) | D l 


(%03) 
(d) (%i4) wximplicit_plot (y*2=x%3-2*x+5, [x,-4,4], [y,-4,4]); 


4 T T T T 
‘| v l 
ab ra. 3 
4 ls a 
(44) 
-1 F 4 
abo | 
ab TN | 
-4 | | | j I | l 
4 -3 -2 A O0 1 2 3 4 
; x 
(%04) 
(e) (Aid) wvzimplicit-plot(y^6=x"2=-x^6; [x,-2,21,. [y2 2]; 
2 
154 | + 
1 | = 
0.5 - 7 
> Or = 
(%t5) 
-0.5 - 
ab J 
1.5 F 7 
-2 | | | i l | | 
2 15 -1 05 0 05 1 15 2 
; x 
(%05) 
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(£) (%i6) wximplicit_plot ((x^2+y^2-1)^3=x^2*y^3, 


2 T T T T 


1.5 F 


ak 
T 


0.5 F 


(%06) 


Solution to Computer activity 42 
(a) (Kil) diff (x*4+5#x%2+15, x); 

(%o1) 4x°+10x 

So the derivative of xt + 52? + 15 is 4z? + 10z. 
(b) (%i2) p:t“2/sin(4+t); 


2 
(%02) sin(4t) 
(4i3) diff (p,t):; 
(403) 2t At? cos(4t) 
= sin(4t) sin(4t)? 
2 2 
So the derivative of a is _ 2t _ 4t cos(4t) 


sin(4t)  sin(4t) sin? (4t) ` 


(c) (hi4) g(u):=(u^2+3)*log(u^2); 
(404) glu) := (u? + 3)log(u?) 
(415) diff(g(u),u); 


2 2 
(%05) 4ulog(u) + Aa 


So the derivative of (u? + 3) In(u?) is 4u In(u) + 


(d) 4X16) ditt (pt 2); 


16t? 2 16tcos(4t) 32t?cos(4t)? 


taoti sin(4t) k sin(4t)  sin(4t)? 
2 


So the second derivative of ———~ is 
sin(4t) 


Solutions to Computer activities 


B22] al: 
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16¢? 2 16tcos(4t)  32¢?cos?(4t) 


sin(4t) Tanay sin? (4t) % sin? (4t) 
(e) (4i7) diff (glu) ,u,3); 


4(u?+3 
(oT) sw) 
A(u? +3 
So the third derivative of (u? + 3) In(u?) is a a 


Solution to Computer activity 43 
(a) 4i1) diff(2*x^6+5/x, x); 


5 
(h01) 12x°- = 
x 


5 


5 
The derivative of 27° + — is 127° — a: 
x x 


(b) (%i2) diff (log(exp(x)+x), x); 
her +1 
her +x 


(%02) 


ef +1 

em +n 

(c). (Rasy ditt tan(t)=t°2) /exp(t), t); 
(403) %e~*(sec(t)? — 2t) — %e*(tan(t) — t?) 


The derivative of ln(e” + x) is 


This can be simplified as follows. 
(4i4) fullratsimp(%) ; 
(%04) —%e*(tan(t) — sec(t)? — t? + 2t) 
(Remember, % in input line 4 refers to the result of the last evaluated 
command.) 
an(t) — t. 


de eh t _ 
The derivative of at is —e~' (tan(t) — sec(t)? — t? + 2t). 


Solution to Computer activity 44 

(a) i1) f(x) :=sin(x)/(exp(x)+1); 
sin(x) 

exp(x) +1 

(b) (%i2) df (x) :=diff (f(x) ,x); 


(%01) f(x) := 


d 
Mad di{x) = ie) 


(c) CALS) ares); 
diff: second argument must be a variable; found 3 
#0: df(x=3) 
-- an error. To debug this try: debugmode(true) ; 


Solutions to Computer activities 


(d) (hid) df (x) :='' (diff (f(x) ,x)); 
cos(x)  %e*sin(x) 


Choa) AB) = Tara werf 


cos(x) _ e” sin(z) 


So f'(x) = a 
(0) f (x) er +1 (er 4 iy 
(e) €415) af(3); 
cos(3) %e?sin(3) 


DO) Kerr (hee +1)? 
(416) float(%): 


(%06) —0.05332658917585964 


So the value of f’ at x = 3 is —0.053 (to 2 s.f.). 
(f) First, define the function. 
(4i7) £ (x) :=sin(x)*exp(cos(x)); 
(%07) £(x):= sin(x)exp(cos(x)); 
Then calculate the derivative, defining df (x) to be the result. 
Chis) df(x):=''(diff(f(x),x)); 
(%o8) df (x) := pe cos(x) — Kes) sin(x)? 
Finally evaluate the derivative at x = 1, giving the result as a decimal 
approximation. 
(419) float (df(1)); 
(%09) —0.2879834260858308 
So the derivative of sin(x)e°8) at x = 1 is —0.288 (to 3 s.f.). 


Solution to Computer activity 45 
(a) (%i1) integrate (2+*x’2t+exp(3*x), x); 
Yes 2x3 


3 °° 3) 


(hot) 


oa 3 
So f 22? + de = +S te 
(b) (4i2) integrate(1/x, x, 2, 3); 
(%02) log(3) — log(2) 
Chis) Float): 


(%03) 0.4054651081081645 
3 
So | de = In(3) — In(2) = 0.405 (to 3 s.f.). 
2 
(c) (%i4) px) :=x/(x%2+1) ; 


(%04) p(x) := 


x 
x? +1 


135 


Computer algebra guide 


(4i5) q(x):='' Cintegrate(p(x) ,x)); 


1 24 
(405) q(x) = reste EN 
2 
So frod = mosh +c. 


Solution to Computer activity 46 
(a) (Ail) integrate((x-3)*(x-1), x); 
x? —6x?49x 
3 
3 Bel 
So [(e-3)e-1) ar = a 


(b) (%i2) integrate(1/sqrt(9-t*2), t); 


(Hot) 


rt 
(%02) asin (=) 
So la = sin! (5) +c. 


(c) (%i3) integrate(exp(t)*sin(t), t, 0, %pi); 


u) +! 
G 2 72 


(%i4) float(%); 
(404) 12.07034631638963 


T T 1 
So | esintdt= — + = 12.1 (to 3 s.f.). 
g 2 2 


Solution to Computer activity 47 


(a) (i1) integrate(exp(sin(x)), x); 
(%01) | ein) dx 


Maxima was unable to calculate this integral. 
(b) (%i2) integrate(exp(-x%2), x); 


(%02) vaert) 


Maxima has calculated this integral in terms of a function called the 
Error function and denoted erf. 


(c) (%i3) integrate(exp(x)/x, x); 
(%03) —gamma_incomplete(0, —x) 


Maxima has calculated this integral in terms of a function called the 
incomplete Gamma function. 
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Solution to Computer activity 48 
(a) (i) (Ail) kill@); 

(%01) done 

(4i2) integrate(x*n, x); 


Is n equal to -1? 
xet 1 


n; 


(%02) 


nti 
(4i3) assume (notequal(n,-1)); 
(%03) 
(hi4) 


[notequal(n, —1)] 
integrate(x^n, x); 
xett 

(%04) 
(i5) 
(%05) 
So Maxima knows n # —1. 

(4i6) forget (notequal(n,-1)); 
(406) 
CAF) facts(}); 


n+1 
facts(); 
[notequal(n, —1)| 


[notequal(n, —1)| 


No facts are now known. 


Solution to Computer activity 49 
First, tell Maxima what we know about a. 
(%il) assume(a>0) ; 

(%01) 
Then do the integration. 


[a > 0] 


(%i2) integrate(1/(x*2ta) ,x); 


atan (=) 
a= 
So 


1 
face- 
ve +a 


where a > 0. 


(402) 


Solution to Computer activity 50 
(a) 


(%11) integrate(log(cos(x*2)), x, 0, 1); 
1 

Cloi) f log(cos(x”))dx 
0 


Maxima cannot evaluate this integral exactly. It returns the integral 
unevaluated. 
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b) (%12) quad_gqags(log(cos(x*2)), x, 0, 1); 
qags (log 
(%02) [—0.1115121024361476, 3.331001582652842 10 **, 21,0] 


(Remember 3.331001582652842 10°11 means 
3.331 001 582 652 842 x 10711.) 


1 
So | In(cos(a?)) dz = —0.112 (to 3 s.f.). (Note that due to the 


0 
inherent differences in how different computer systems compute 
floating-point arithmetic, there may be differences in the final few 
decimal places between your result and that shown here.) 


(c) (%i3) quad_qags(exp(-t%2), t, 0, 1); 
(403) [0.7468241328124271, 8.291413475940727 10715, 21, 0] 


I 
So f et dt = 0.747 (to 3 s.f.). 
0 


Solution to Computer activity 51 
(a) (411). Pematrix( [2-85], (2.1.2), [1,-1.4] 3; 


2 —3 3 
(%01) |-2 1 2 
1 -1 4 
(b) (412) Q:matrix( [5,-3], [-1,-2], [0,4] ); 
5 —3 
(%02) |—-1 -2 
0 4 
(c) Gas? Rematriz( [iy-2,-1],. 4.7.0), [5,43] J; 
1 -=2 =f 
(%03) 4 T 0 
—5 1 3 


Solution to Computer activity 52 
(a) (%i4) matrix_size(Q) ; 
(%04) [3,2] 
So the matrix Q has size 3 x 2. 
(b) 415) RIS;21; 
Cos). d 
The element in the third row and second column of R is 1. 
(c) (hi6) RIS: 
(%06) [—5,1,3] 
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(d) 


Chit) RL3,21 20: 
(%07) O 
(418) R; 
1 —2 -1 
(%08) 4 7 0 
—5 0 3 


Solution to Computer activity 53 


(a) 


(Lil) P:matrix([2,0],[3,-1],[-4,6]); 


2 (0) 
(hol) 3 -1 
—4 6 


(4i2) Q:matrix([0,3],(5,-4],[1,3]); 


0 3 
(%02) [5 —4 
1 3 


(413): Rematrix( [1,3]; (2,41): 


, 1 3 
(403) : al 


Cia) Ssmatrix([1]; (21): 


404) 5) 


(415) - Timatrix¢[3;4))% 
(%05) [3 4] 
(%16) P+Q; 


2 3 
(%06) 8 —5 
—3 9 


2 3 
So P+Q= 8 —5 |. 
—3 9 


Cait) PR; 
fullmap: arguments must have same formal structure. 
-- an error. To debug this try: debugmode(true) ; 


P — R cannot be calculated, since the matrices are of different size. 


An error is given. 
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(d) (%i8) 2*P; 


4 (0) 
(%08) 6 -=2 
—8 12 
4 0 
So 2P = 6 —2 
—8 12 
(e) (%i9) P.R; 
2 6 
(%09) |1 5 
8 12 
2 6 
SoPR=Į|1 5 
8 12 


(f) (Aino) Pq; 


MULTIPLYMATRICES: attempt to multiply nonconformable 
matrices. 
-- an error. To debug this try: debugmode(true) ; 


PQ cannot be calculated, since the matrices are not of the 
appropriate sizes. An error is given. 


(g) i11)P.S; 


2 
(4o11) |1 
8 


2 
So PS= |1 
8 


(h) (%i12) P.T; 
6 


(fot): |S 
12 


Maxima interprets the row vector T as a column vector and obtains | 5 
12 


(i) (i13) R™2; 
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(%i14) R'2; 


1 9 
(4014) 5 En 


1 

y whose elements are the squares 
4 16 
of corresponding elements of R. 


(%115) P*Q; 


Maxima returns the matrix ( 


0 (0) 
(%015) |15 4 
—4 18 
0 0 
Maxima returns the matrix 15 4] whose elements are the 
—4 18 


products of corresponding elements of P and Q. 


Solution to Computer activity 54 


(a) 


(4i1) A:matrix([1,2,0],[-4,4,2], [-1,2,1]); 


1 2 0 
(fol) |—4 4 2 
—1 2 1 


(12) Bematrix( [1.2.0] ;(5,6, 21, [2,2,41): 


1 2 0 
(%02) |5 6 2 
2 2 1 


(%i3) determinant (A) ; 


(403) 4 
So det A = 4. 
(4i4) C:invert (A) ; 
o -i 1 
hoa) | a 73 
=1 -1i 3 
0 —4 1 
soat=[ a 4 4 
— -1 3 
(415) AC; 


1 0 0 
(%05) jO 1 O 
O O 1 
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(4i6) C.A; 


Oo ©: 


1 (0) 
(%06) 10 (0) 

0 1 
So AC = CA =I, and C is the inverse of A. 
(4i7) invert (B); 


expt: undefined: 0 to a negative exponent. 
-- an error. To debug this try: debugmode(true) ; 


B does not have an inverse. 


Solution to Computer activity 55 


(a) 
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First, enter the matrix into Maxima, and assign it to a variable. 


(Ali) Aimatrix([2,-1,1], [2,0,11]; [4;2;1]); 


2 —1 1 
(%01) |2 0 1 
4 2 1 


Then, check if the matrix is invertible by finding its determinant. 
(%i2) determinant (A) ; 
(hoa) —2 
We have 
2 -1 1 
det | 2 0 1] =-2, 
4 2 1 
so this matrix is invertible. 
Now find the inverse. 


(413) invert(A); 


1 —3 1 
2 2 
Mos) |-1 1 (0) 
—2 4 —1 
1 3 1 
2 -1 1 1 =3 5 
So | 2 0 1 = |-l1 1 0 
4 2 1 —2 4 —1 


Again, first enter the matrix and calculate the determinant. 


(%i4) B: matrix( [-1,0,3,-2], [4, 1/2, -1 ,0], 
[2,0,1,-1], [-1,2,1,-1] ); 


—1 3 —2 
—1 0 


1 —-1 
1 =i 


N ONIF O 


(hoa) |5 
= 


(%i5) determinant (B); 


(%05) : 
oO 5 
We have 
—1 0 3 —2 
1 
5 —l 0 1 
det 2 = —, 
0 1 —1 2 
—] 2 1 —1 


so this matrix is invertible. 
Now find the inverse. 


(416) invert (B); 


—4 —4 7 1 
: 6 -6 10 2 
(06) 19 —20 33 5 
—27 —28 46 7 
zi 0 3 -A ! A A 71 
del Bah 0 _ | =o, <6 1-2 
“lo 9 a tl ~ | -19 0 33 5 
-1 2 1 -1 -27 —28 46 7 


Solution to Computer activity 56 
(a) (%i1) a[n] :=0.8^n; 
(7,61). an:= 0.8" 
(b) (%12) 21100); 
(%02) 2.037035976334508 10 +° 
(Remember, this means 2.037 035 976 334508 x 10 1.) 
So a100 = 2.037 x 1071? (to 4 s.f.). 


Solution to Computer activity 57 
(a) La DEAL: 
(fo1) 1 
(4i2) bln] :=2*b[n-1]+1; 
(%02) By :=2b,-14+1; 
(b) (%i3) b[5]; 
(%03) 31 
(%i4) b[100]; 
(%04) 1267650600228229401496703205375 
So bs = 31 and b100 = 1 267 650 600 228 229 401 496 703 205 375 
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(c). (id) 
(%05) 
(hi6) 
(%06) 

(d) i7) 
(%07) 
(%18) 
(%08) 
(hi9) 
(%09) 


BLA +10; 

10 

DISI; 

31 

kill(b) ; 

done 

BELI LO; 

10 

b[n] :=2*b[n-1]+1; 
bn := 2 bn—1 + 1; 


(4110) b [5]; 
(4010) 175 
So bs is now 175. 


Solution 
(a) (hii) 
(%o1) 
(b) hi2) 
(02) 
(c) (%i3) 


(%t3) 


(%03) 


Solution 
(a) i1) 
(hot) 
(%12) 
(%02) 
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to Computer activity 58 

e[nl.s=n°2: 

=n 

ptsal Lem]; [2sel211; Bsc] 1; 

(lt, 1], [2,4], [3, 91) 

wxplot2d([discrete, pts], [style, points], 
[xlabel, "n"], [ylabel, "cln]"]); 


de) 


c[n] 

NN WO Ff ODN OO 
T 
ri 


to Computer activity 59 
ali]:0; 
0 
d[n] :=1/2*(2-d[n-1]); 
1 


dy = 5 (2— eee 


Solutions to Computer activities 


(b) (fag) pts:makelist([n,d[n]], n, 1, 15); 
1 3 5 11 21 43 85 


(%03) [[1, 0], [2, 1], [3, alt [4, al [ j 3l | ’ T | , T. [8, TL [9, o! 
171 341 683 1365 2731 5461 
[10, sl (fa, ——]; [12, = [13, sl [14, ~ [15, ——|] 
256 512 1024 2048 4096 8192 
(c) (%i4) wxplot2d([discrete, pts], [style, points], 
[xlabel, "n"], [ylabel, "d[n]"]); 
1 ® T T T T T T 
0.8 F ê 7 
SC eeeeseeee 
= 0.6 + e zi 
= ° 
ta) ° 04+ J 
0.2 + J 
0 2 j | | | L it | 
0 2 4 6 8 10 12 14 16 
n 
(ho4) 
Solution to Computer activity 60 
(a) First, define the sequence. 
Chad) fili; 
(hol) 1 
(412) xr[n] :=1.5*r[n-1]-0.1*(r[n-1] )%2; 
(J62) Ta:= 15744 — Gee, 
Next, create a list of points to plot, suppressing the output. 
(%i3) rpts:makelist([n,r[n]],n,1,50)$ 
Finally, plot the graph. 
(4i4) wxplot2d([discrete, rpts], [style, points], 
[xlabel, "n"], [ylabel, "r[n]"]); 
5 EER a | 
45 F e J 
4+ ®° J 
35+ ° J 
= 3p ° - 
(44) 25+ e 
ae J 
1.5 e 4 
1 le L L | | | | | L | 
0 5 10 15 20 25 30 35 40 45 50 
n 
(ho4) 
The terms of the sequence seem to be getting closer and closer to the 
value 5. 
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(b) First, define the sequence. 


(%15) s[n] :=n%*3-10; 

(%05) Sn := n° — 10 

Next, create a list of 50 points to plot, suppressing the output. 
(%16) spts:makelist([n‚s[nll‚n,1,50)$ 

Finally, plot the graph. 


(4i7) wxplot2d([discrete, spts], [style, points], 
[xlabel, "n"], [ylabel, "s[n]"]); 


140000 _ rrr 
120000 
100000 
80000 
60000 
Cht7) 40000 
20000 

0 


-20000 1 1 1 1 1 1 1 1 1 
0 5 10 15 20 25 30 35 40 45 50 


n 


s[n] 


oT) 

The terms of the sequence seem to be continually increasing. 
First, define the sequence. 

(%i8) t[n]:=1-1.1^n;} 

(408) ta:= 1— 1.1” 

Next, create a list of 50 points to plot, suppressing the output. 
(%i9) tpts:makelist([n,t[n]],n,1,50)$ 

Finally, plot the graph. 


(%i10) wxplot2d([discrete, tpts], [style, points], 
[xlabel, "n"], [ylabel, "t[n]"]); 


(%t10) 


-120 | | | | | | | | | ® 
0 5 10 15 20 25 30 35 40 45 50 


n 


(%010) 


The terms of the sequence seem to be continually decreasing. 
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Solution to Computer activity 61 


(a) 


(%4i1) sum(n’2, n, 1, 100); 
(701) 338350 
100 
So $ n” = 338 350. 
n=l 
20 i 
The sum can be written `> (5) and so can be calculated as follows. 
n=1 
(%12) sum((1/3)*n, n, 1, 20); 
1743392200 
(ho2) =z 
3486784401 
20 PA 
1 1 743 392 200 
EG) En 
= 3 3 486 784 401 


Notice that the sum is a geometric series with first term a = $, 
common ratio r = i and with 20 terms. So the sum can be calculated 
using 


a- HH) 


l-r 1— 


Del = wamne 


wo 
Do 
© 


— 1 
2 x 320 
3 486 784 400 


2 x 3486 784 401 
1 743 392 200 


~ 3486 784401 
~ 0.499 999 999 8566 


CALS -sumCl/ (4n); md inf); 


(%03) 5 
n=1 


(%i4) simpsum: true; 


|= 


n 


EN 


(%o4) true 


(%i5) sum(i/(4%n), n, 1, inf); 
(105) = 
1 1 : : : 
So 2. ae You could check this result using the expression for 


n=1 
the sum of an infinite geometric series. 
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(e) (%i6) sum(1/(n*2), n, 1, inf); 
2 


(406) £ 
oO B 


Solution to Computer activity 62 
(a) G1) tayler(cos (a), x, Ypi/e, 3); 


ag SD ed 


(401) /T/ EE z ' T 


So the cubic Taylor polynomial about 7/6 for f(x) = cos(x) is 


ast les) E-a 


z 2 4° @ 
(b) AAL plat Goaylor Cian), x; 0; BI): 
3B 955 
(402)/T/ p(x) := x + = + Ie: + 
So the quintic Taylor polynomial about 0 for f(x) = tan(x) is 
x 2x 
LL 3 + T5 


(e) (iis) py); 

2y° +5y° +15y 
—— 
(d) (%i4) float(p(0.1)); 


rat: replaced 0.10033466666667 by 75251/750000 = 
0.10033466666667 


(404) 0.10033466666667 
So p(0.1)= 0.10033 (to 5 s.f.). 
(%i5) float(tan(0.1)-p(0.1)); 


rat: replaced 0.1003346666666667 by 75251/750000 = 
0.1003346666666667 


(403) /R/ 


rat: replaced 0.1003346720854505 by 9334601/93034649 = 


0.1003346720854507 
(%05) 5.418783991614422 10 ° 


So the difference between the function and its quintic Taylor 
polynomial at x = 0.1 is 5.42 x 107° (to 3 s.f.). 


Solutions to Computer activities 


Solution to Computer activity 63 
First, define the function. 

(il) £(x) :=exp(sin(x’%2)); 

(%o1) f(x) := exp(sin(x?)) 


Next, calculate the required Taylor polynomial. The command is ended 
with $ to suppress displaying the result. 


(432) p(x):='" (taylor (f(x) ,x,1,5))$ 
Finally, plot the two curves over the required range. 
(413) weplot2d( (f(x) ,p(z)], [k,0,21, (legend, "£G@)", "“ptx)"]); 


4 T T j 


(%t3) 


(403) 


Solution to Computer activity 64 
(a) (Kil) a:3+4*%i; 
(hoi) 4%i1+3 
(b) (4a2) b:2=2*/i; 
(%02) 2—2%i 
(c) C%i3) atb; 
(%03) 2%i+5 
(%i4) a-b; 
(404) 6%i+1 
So a+b= 5 + 2i and a-b= 1 + 61. 
(%i5) a*b; 
o5) (2— 2%i)(4%i +3) 
(%16) a/b; 
4%i+3 
2—2%i 


(406) 
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(d) (%i7) rectform(a*b) ; 
(%07) 2%i +14 
(%18) rectform(a/b) ; 
Thi 1 


4 4 
So a*b= 14+ 27 and a/b= 4 + Ti. 


(%08) 


(e) (%i9) realpart (a); 

(%09) 3 

(%i10) imagpart (a) ; 

(%010) 4 

So the real and imaginary parts of a are 3 and 4 as expected. 
(f£) i11) conjugate (a) ; 

(%o11) 3— 4%i 

So the conjugate of a is 3 — 4i as expected. 
(g) (%i12) abs(b); 

(4012) 23/2 

(%i13) carg(b) ; 


(013) -Z 
(te) 4 


So |b| = 2/2 = (/2)° and the principal argument of b is 7 


Solution to Computer activity 65 
(a) i1) c:3*(cos(5*%pi/13)+,i*sin(5*,pi/13)) ; 


hon 3(tisin( =) + cos ( 2 
Re) ‘isin | 55 cos | 73 


(%i2) d:4*(cos(7*/,pi/13)+hi*sin(7*/pi/13)) ; 


(402) 4(%isin( Z = 
cO hl sin 13 + cos 13 


(b) The modulus of c*d can be found as follows. 


(%13) abs(cxd) ; 


(03) 12 . Br \ ° 4 5 \? . Tx \? + Tr\? 
A sın = COS en sin — COS — 
as 13 13 13 13 


Simplify this using the trigreduce command. Remember, % refers to 
the result of the previous calculation. 


(4i4) trigreduce(%) ; 
(%04) 12 
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The principal argument of c*d can be found as follows. Again the 
initial result requires simplification. 


(%15) carg(cxd) ; 
(=) 
sin | — 
13 
TT ST 
cos | — cos | — 
13 13 


(%i6) trigreduce(%); 


(%05) atan 


: 127 
(%o6) 3 


12 
So |c*d| = 12 and arg(c*d) = ce 
Similar methods can be used to find the modulus and principal 
argument of c/d. 


(Zit) absc/d): 


(%07) 


(%18) trigreduce(%); 
(%08) : 
00 4 


(%i9) carg(c/d) ; 
(409) —atan 


(%110) trigreduce (%) ; 

o 247 

(4010) T3 
Note that although carg gives a principal argument, that is an 
argument in the interval (—r, r], the trigreduce command simplifies 
the expression to an equivalent argument which is outside this 


interval. 
The principal Be 
e principal argument is -73 > 
2 
So |c/d| = ł and arg(c/d) = -= 
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Solution to Computer activity 66 


(a) (%i1) solns:solve(z*3+z+1=0) ; 


val ay veut 4 oS 
wo (ED (ES) 
23 2 2 2 Wer i 
233/2 2 
1/3 o _V8h%i 1 
= V31 1 V3%i 1 5 5 
LT (2372 2 23 32) 7 = P 
v31 1 
233/2 2 
Z r 1 
233/2 2 JA 1 1/3 
233/2 2 


(b) (%i2) float(rectform(solns)); 


(%02) [z = 0.34116390191401 — 1.161541399997252 Yi, 
z =1.161541399997252%i + 0.34116390191401, 
z = —0.68232780382802] 


So, to two decimal places, the solutions are 
z = 0.34 — 1.16i, z=0.344+1.16;, and z= —0.68. 
(c) (%i3) solve(z^4-20*z^3+171*z^2-626*z+962=0) ; 
(403) [z =3-—2%i,z = 2%i+3,z =7-—5%i,z=5%i+7] 


So the solutions are 3 — 2i, 3 + 2i, 7 — 5i and 7 + 5i. 


Solution to Computer activity 67 
(a) (i1) s:solve(z^4+4*z+5=0) ; 


vV4— 4i V4—4%hi 


(401) zi, 
de V&hi+4 Vit ,. 
Zhi En H i] 


(b) (%i2) v:makelist(rhs(s[k]), k, 1, length(s)); 
V4—4%i va Ali 


(%02) [—*i- 5 a — Ji, 
vi V4hi+4 VEE a 
0 EB oe 0 


(c) (4i3) pts:makelist([realpart(v[k]), imagpart (v[k])], 


k, 1, length(s)); 
28/2744 Vala Vga Vala 


(%03) [[ 93/2’ 9372 ‚| 23/2 ? 23/2 I 


y 25/2 +4 l= y 25/2 — 4 CEEE Vo= i 


[- 93/2” 93/3 l, 93/2 ? 93/2 
d) (414) wxplot2d([discrete,pts], [style,points], 
P P yte,p 
[xlabel,""], [ylabel,""]); 


1.5 T T T T IK 
gl | 
0.5 L e 4 
(iy Or 7 
-0.5 F e 4 
Pal | 
1.5 Le 


(104) 


Solution to Computer activity 68 
(Ail) sisolve(2*7=1): 


Qhint 4 him 6 hin 

(iol) [z=%e 7 zhe T zhe T7 3 
6 hin 4 hin 2 hir 
z=je 7 z= he. 7 z= he. a. 


z= Í] 
The solutions are given in terms of exponential functions. You will learn 
about the exponential form of complex numbers in Section 4 of Unit 12. 
The solutions can be written in Cartesian form as follows. 
(⁄12) float(rectform(s)); 


(%02) [z = 0.78183148246803 %i + 0.62348980185873, 
z = 0.97492791218182 %i — 0.22252093395631, 
z = 0.43388373911756 %i — 0.90096886790242, 
z = —0.43388373911756 hi — 0.90096886790242, 
z = —0.97492791218182 %i — 0.22252093395631, 
z = 0.62348980185873 — 0.78183148246803 %i, 
z = 1.0] 


So, to two decimal places, the seven seventh roots of unity are 
z = 0.62 + 0.78i, z = —0.22 + 0.97i, z = —0.90 + 0.43i, 
z = —0.90 — 0.43i, z = —0.22 — 0.97i, z = 0.62 — 0.78i 
and z = 1.00. 


Solutions to Computer activities 
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To plot the solutions in the complex plane, first create a list containing 
just the solutions of the equation (without the z=). The command is ended 
with $ to suppress the output. 


(4i3) v:makelist(rhs(s[k]), k, 1, length(s))$ 
Then create the list of points to be plotted. 


(hi4) pts:makelist([realpart(v[k]), imagpart(v[k])], 
k, 1, length(s))$ 


Finally, plot the points. 


(%i5) wxplot2d([discrete,pts], [style,points], 
[xlabel,""], [ylabel,""]); 


NE il: cilia 4 


-1 1 L 1 [| fi Í L L fi 
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 06 08 1 


(%05) 


Notice that the complex solutions occur in conjugate pairs, and hence lie 
symmetrically about the horizontal axis. The solutions also lie on the unit 
circle (which does not seem circular here due to the different axis scales 
used). 
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Mathematical operations and functions 


Operation or function Syntax Example 
Addition + 4+3; 


Subtraction 
Multiplication 
Division 

Brackets 

Powers, for example 23 


Square root, for example v5 
Exponential, for example e° 
Natural logarithm, ln 


Magnitude (absolute value) 
of a real number, 
or modulus of a complex number 


Equality (within an equation) 


Matrix addition 

Matrix subtraction 

Scalar multiple of a matrix 
Matrix multiplication 

Matrix powers 

Determinant (of a square matrix) 


Matrix inverse (of a square matrix) 


Identity matrix 


AA 
determinant (matrix ) 


invert (matrix ) 
or matrix ^^(-1) 


ident (size ) 


4-3; 
4*3; 
4/3; 

2* (3+4); 
2^3: 
2**3; 
sqrt(5); 


Ke 3; 
exp (3) ; 
log (8); 


abs (-3) ; 
abs (1+2*%1) ; 


sin(1); 

cos (3*/pi/2) ; 
tan (%pi/4); 

esc (2) ; 
sec(%pi/3) ; 

cot (%pi/4) ; 

asin (sqrt (3) /2) ; 
acos(1/sqrt (2); 
atan(1/2) ; 


Pig; 
determinant (P); 
invert (P); 

PAN (-1) ; 

ident (2); 


Maxima reference guide 
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Mathematical operations and functions (continued) 


Operation or function Syntax Example 
Real part of a complex number realpart( ) realpart (1+2*/i) ; 


Imaginary part of a complex number | imagpart(_) imagpart (1+2*%i) ; 

Conjugate of a complex number conjugate(_) conjugate (1+2*%i) ; 

Principal argument of a carg ($) carg(1+2*%i); 
complex number 

Express a complex number in rectform( ) rectform(2/%i) ; 


Cartesian form 


Mathematical commands 


Operation Command Example 
Get help on a command ? command ? float; 


command 


Load a package 


exponentials and logarithms 
Combine logarithms 


156 


load( package name ) 


logcontract( ) 


or describe ( command ) describe(float) ; 
Find help information whose ?? ?? flo; 
title contains the given text | or describe( , inexact) describe(flo, inexact) ; 
Force the evaluation of a ED g(x) :='' (diff (x4,x)); 


load(implicit_plot); 


Assign a value to a variable variable ; value a:23; 
Display the value of a variable | variable a; 
Define a function := f (x) :=2*x+3; 
Evaluate a function at a value | function ( ) f(1); 
Remove an assigned variable | kill( variable ) kill(a); 

or function kill(function) kill(£); 
Remove all assigned variables | kill (all) kill(all) ; 

and functions 
Reset all system variables reset () reset (); 
Convert to a decimal number | float ) float (sqrt (2)) ; 
Expand brackets expand( _) expand ((x+1) 2) ; 
Factorise factor ($) factor (2*x+4*x^2); 
Simplify fullratsimp( ) fullratsimp ((2*x+4*x%2) /x) ; 
Simplify something involving | radcan(_) radcan(log(x*2*y)) ; 


logcontract (log (a) +log (b)) ; 


Mathematical commands (continued) 


Operation Command Example 


functions of sums, 
differences and multiples of 
angles 

Express powers of sines and 
cosines in terms of sines and 
cosines of multiple angles 

Simplify trigonometric 
expressions 

Simplify algebraic fractions 
containing trigonometric 
functions 

Substitute 


Left-hand side of an equation 
Right-hand side of an equation 
Solve an equation exactly 
Solve an equation for a 
variable, i.e. change the 
subject 

Solve simultaneous equations 


Solve an equation numerically 


Form a list 
Assign a list to a variable 
Create a list 


Extract an element of a list 
Find the length of a list 


Specify a matrix 
Size of a matrix 
Row of a matrix 
Element of a matrix 


Expand trigonometric trigexpand( ) trigexpand(sin(A+B)) ; 


trigreduce( ) 


trigsimp( ) 


trigrat( ) 


subst (value , variable , expression ) 


lhs ( equation ) 
rhs ( equation ) 
solve( equation ) 


solve( equation, variable) 


solve( list of equations , 
list of variables ) 


find_root (expression , variable , 
interval start value, 
interval end value) 


IA M 
vyemisiole ii i iz 
makelist( general term , 
dummy variable , start value, 
end value) 
list [index ] 


length( list ) 


matrix(row,row,...) 
matrix_size (matrix ) 
matrix [row] 

matrix [row , column] 

or matrix [row ] [ column] 
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trigreduce(sin(x)%2) ; 


trigsimp(sin(x)2+cos(x)%2) ; 


trigrat ((sin(2*x) )/cos(x)); 


subst (4,x,x2+1); 
which substitutes 4 for x 
in x2+1 


lhs (4*x+1=2*x-2) ; 
rhs (4*x+1=2*x-2) ; 
solve (2*x’*2-1=0) ; 
solve(2*a*b-3*b=0, a); 


solve([2*x+y=4, x-2*y=1], 
[x,y]); 


find root (2*x*3-1,x,0,1); 


[a,b,c]; 
A: [a,b,c]; 
makelist(2*n,n,1,100); 


A[2]; 
length(A); 


matrix([1,2],[3,4]); 
matrix_size(A); 
A[2]; 

A[3,4]; 

A[3] [4]; 
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Operation Command Example 
Define a sequence using a sequence [n]:= expression in n a[n] :=n^2; 


closed form 
Define a sequence using a 
using a recurrence system 


Calculate a term in a sequence 
Sum a series 


Differentiate 
Differentiate multiple times 


Integrate 
(find an antiderivative) 
Evaluate a definite integral 


Find an approximate value 
of a definite integral 


Find a Taylor polynomial 


Make an assumption about a 
variable 

State two things are equal 

State two things are not 
equal 

Forget a property 

List all known facts 

List all known facts about 
a particular variable 


Mathematical commands (continued) 


sequence [first term number]: 
initial value 

sequence [n] := 
expression in sequence [n-1] 

sequence [term number ] 

sum( general term, index variable , 
lower limit , upper limit ) 


diff ( 
diff ( 


, variable ) 
, variable , positive integer ) 


integrate( , variable) 


integrate( , variable , 
lower limit , upper limit ) 
quad_qags( , variable , 
lower limit , upper limit ) 


taylor( function , variable , 
centre ,|degree ) 


assume() 


equal( , ) 
notequal( , ) 


forget ( property ) 
facts() 
facts (variable ) 


biti: i; 
b [n] :=2*b[n-1] ; 


b[50] ; 
sum(n*2,n,1,10); 


diff (sin(x*2) ,x); 
diff (log(x),x,3); 

to differentiate three times 
integrate (x*2,x) ; 


integrate(sin(x),x,0,1); 


quad_qags(exp(-x*2) ,x,0,1); 


taylor (sin(x) ,x,0,3); 


assume (a>0) ; 


assume (equal (n,3)) ; 
assume (notequal (n,-1)) ; 


forget (a>0) ; 
facts); 
facts(a); 
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Mathematical commands (continued) 


Operation Command Example 


Plot several graphs 


Plot a sequence 


Plot a curve represented by 
an equation in implicit 
form 


Plot several curves 
represented by equations 
in implicit form 


Set the default plotting 
options 


Plot a graph wxplot2d( expression , wxplot2d(x*2, [x,0,1]); 


horizontal range ,...) 


wxplot2d( list of expressions, 
horizontal range ,...) 


wxplot2d([discrete, 
list of coordinate pairs], 
[style, points],...) 


(use plot2d similarly if not 
using wxMaxima) 


wximplicit_plot (equation , 
horizontal range , 
vorcical Teneo peN) 


wximplicit_plot( 
list of equations , 
horizontal range , 
vertical range, ...) 


(use implicit_plot similarly if not 
using wxMaxima) 


Note: the wximplicit_plot and 
implicit_plot commands require 
the implicit_plot package to be 
loaded. 


set_plot_option( option ) 
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which plots the graph of z? 
for x between 0 and 1 
wxplot2d (x2, [x,0,1],[y,0,2]); 
which plots the graph of z? 
for x between 0 and 1, and with 
vertical axis between 0 and 2 
wxplot2d([x*2,2*x],[x,0,1]); 
which plots graphs of x? and 2x 
for x between 0 and 1 


wxplot2d([discrete, 


[{1,5], [2,10]]], 
[style, points]); 


wximplicit_plot (x*2+2*y*2=1, 


[x,-1,1], [y,-1,1]); 


wximplicit_plot( 
[x%2+2*y*2=1, y=x^3], 
[x,-1,1], [y,-1,1]); 


set_plot_option(same_xy) ; 
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Other mathematical symbols 


Constant | Syntax | Example 

e he he^2; 

T Lpi Q*hpi; 

i hi 2+3*%i; 

oo inf sum(1/n’*2,n,1,inf); 
Maxima system variables 

Variable Description Example 
fpprintprec | The number of significant figures of | fpprintprec:3; 
a decimal number to display 

functions The list of all user-defined functions | functions; 
values The list of all user-assigned variables | values; 
simpsum Whether to simplify summations simpsum: true; 


Maxima graph plotting options 


The following options can be added as additional arguments to plotting 


commands. 


Option Argument Example 


Curve colour 

Legend text 

Turn legend off 

Position legend at the bottom of 
the graph 

Position legend to the right of 
the graph 

Horizontal axis label 

Vertical axis label 

Set axes to have equal scales 


Change the size of the graph 
Plot points 
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Vertical range [y. E. Hp [y, 0, 5] 


[color, aad 
(legend, " ", ...] 
[legend, false] 
[gnuplot_preamble, 
"set key bottom"] 
[gnuplot_preamble, 
"set key outside"] 
[xlabel, " "] 
[ylabel, " "J 


[same_xy, true] 


wxplot.size:[ , ] 


[style, points] 


Maxima packages 


[color, red] 
[legend, "Car A"] 
[legend, false] 
[gnuplot_preamble, 
"set key bottom"] 
[gnuplot_preamble, 
"set key outside"] 
[xlabel, "t (h)"] 
[ylabel, "s (km)"] 


[same_xy, true] 


wxplot_size: [800,600] ; 
[style, points] 


Maxima package | Functions added 


implicit_plot Plotting curves given by equations in implicit form 


wxMaxima keyboard sequences (Microsoft Windows) 


Function 


Zoom out Alt-0 

Copy Ctrl-C 
Interrupt a calculation Ctrl-G 
Re-evaluate the worksheet | Ctr1-R 
Close wxMaxima Ctr1-Q 
Paste Ctrl-V 


Key sequence 


Zoom in Alt-I 


wxMaxima menu commands (Microsoft Windows) 


Close wxMaxima File Exit 

Save your work File Save or Save As 

Print the worksheet File Print 

Open a worksheet File Open 

Configure wxMaxima Edit Configure 

Copy Edit Copy 

Paste Edit Paste 

Zoom in View Zoom In 

Zoom out View Zoom Out 

Toggle display of main toolbar | View Main Toolbar 
Re-evaluate worksheet Cell Evaluate All Visible Cells 
Insert text cell Cell Insert Text cell 

Insert title cell Cell Insert Title cell 
Insert section cell Cell Insert Section cell 
Insert subsection cell Cell Insert Subsection cell 
Interrupt a calculation Maxima Interrupt 

Restart Maxima Maxima Restart Maxima 

Enter a matrix Algebra Enter Matrix... 
Maxima help window Help Maxima Help 


Maxima command-line interface commands 


Operation 


Interrupt a calculation 

Turn off 2D formatting of maths 

Display output left justified 

Start writing a transcript 

Stop writing a transcript 

Playback all the input and output in 
a session 

Playback a range of input and output 


Save the state of a session 
Reload a session 


Close Maxima quit(); 


Command 


Ctrl-C 

display2d:false; 
leftjust:true; 
writefile(" file path") ; 
closefile(); 

playback () ; 


playback([ start line number , 
end line number ]); 

save("file path", all); 

load("file path"); 
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abort calculation 
using the command-line interface 95 
using wxMaxima 19 
abs 14 
for complex numbers 85 
absolute value 14 
accessibility 89 
changing wxMaxima colours and fonts 90 
using a screen magnifier 91 
using a screenreader 91 
using wxMaxima with keyboard alone 90 
acos 47 
adding text to a wxMaxima worksheet 26 
addition 14 
of matrices 71 
annotating a wxMaxima worksheet 26 
approximation, numerical 50 
argument of a command 14 
argument of a complex number 85 
asin 47 
assign 
a derivative to function 60 
an integral to function 63 
functions 42 
variable 22 
assume 65, 66 


assumption 
about a variable 65, 66 
equals 66 
forget 65, 66 
list all 66 
not equals 65, 66 
atan 47 


axis scales, on a graph 55 
brackets 14 


carg 85 
Cartesian form of a complex number 85 
CAS see computer algebra system 
cell, in wxMaxima 16 
compress 16 
delete 17 
marker 11, 16 
re-evaluate 20 
section 26 


Index 


subsection 26 
text 26 
title 26 


changing a sequence 77 

changing the size of a graph, in wxMaxima 41 
circles, plotting 54 

closed form of a sequence 75 

closefile 98 

closing 


command-line interface 94 
wxMaxima 13 


colour of graph 40 
command 


argument of 14 

editing 
using the command-line interface 95 
using wxMaxima 20 

force evaluation of 61 

reusing 
using the command-line interface 95 
using wxMaxima 20 


command-line interface 7, 91 


abort calculation 95 
closing 94 

creating transcript 97 
cursor, editing 92 
editing commands 95 
editing cursor 92 
graphs 100 

input line number 92 
interrupt calculation 95 
output line number 93 
plotting graphs 100 
printing 97 

reusing commands 95 
saving work 96 
transcript 97 


complex numbers 35, 84 


argument (principal) 85 
Cartesian form 85 
conjugate 85 
imaginary part 85 
modulus 85 

polar form 85 

real part 85 
rectangular form 85 
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compress cell, in wxMaxima 16 
computer algebra system 5 
configure wxMaxima 12, 16 
conjugate 85 

conjugate of a complex number 85 
copy, using wxMaxima 21 


cos A7 
cosec 47 
cot A7 


create a list using makelist 80 
csc 47 
cube roots 15 
cursor 
editing 
in command-line interface 92 
in wxMaxima 20 
horizontal, in wxMaxima 20 


decimal approximation of a number 17 
define 
functions 42 
variables 22 
definite integral 62 
approximate value 66 
delete 
cell, in wxMaxima 17 
function 43 
variable 24 
derivative 59 
assign to function 61 
second, third 60 
describe 31 
determinant 74 
determinant of a matrix 74 
diff 59 
differentiate 59 
division 14 


editing commands 
using the command-line interface 95 
using wxMaxima 20 
editing cursor 
using the command-line interface 92 
using wxMaxima 20 
element 
index 36 
of alist 36 
of a matrix 71 
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equal 66 
equals 34 
assumption 66 
equations 34 
changing the subject 37 
left-hand side 34 
manipulating 37 
plotting 39 
rearranging 37 
right-hand side 34 
simultaneous 38 
solving 34 
solving numerically 50 
trigonometric 49 
erf 64 
evaluate 
an expression 11 
an expression using the command-line 
interface 93 
force command to 61, 62 
re-evaluate a cell using wxMaxima 20 
re-evaluate worksheet using wxMaxima 
exp 46 
expand 33 
expanding brackets 33 
exponentials 44 
expressions, evaluating 11 


factor 33 

factorising 33 

facts see assumption 

facts 66 

file types, wxMaxima 27 

find_root 50 

float 17 

floating-point number 17 
displayed precision 25 

forget 65, 66 

forget assumption 65, 66 

fpprintprec 25 

front-ends 7 

fullratsimp 33 

functions 42 


assign 42 
define 42 
delete 43 


trigonometric 47 
trigonometric, inverse 47 
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functions 43 


gamma_incomplete 64 

graphs 
changing default behaviour 56 
changing line thickness 100 
changing properties 100 
colour 40 
equal scales 55 
from equations in implicit form 54 
legend 40 
plotting 39 
size of, in wxMaxima 41 
using command-line interface 100 
vertical range 39 


help 30 


horizontal cursor, in wxMaxima 20 


ident 75 

identity matrix 75 

imaginary numbers 35 

imaginary part of a complex number 85 
imagpart 85 


implicit form of an equation, plotting 54 


implicit_plot 
command 54 
package 54 
indefinite integral 62 
index (of a list element) 36 
inf 82 
infinite sequence 80 
infinity 82 
input line number 11 
input prompt, wxMaxima 11 
inputting a matrix 68 
installing Maxima 8 
integral 62 
assign to a function 63 
definite 62 
definite, approximate value 66 
indefinite 62 
integrate 62 
integrate 62 
interfaces 7 
command-line 7, 91 
wxMaxima 7,9 
interrupt calculation 
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